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Abstract. The description of strongly interacting quantum fields is based on two- 
particle irreducible (2PI) approaches that allow for a consistent treatment of quan- 
tum systems out-of-equilibrium as well as in thermal equilibrium. As a theoret- 
ical test case the quantum time evolution of <? 4 -field theory in 2+1 space-time 
dimensions is investigated numerically for out-of-equilibrium initial conditions 
on the basis of the Kadanoff-Baym-equations including the tadpole and sunset 
self energies. In particular we address the dynamics of the spectral ('off-shell') 
distributions of the excited quantum modes and the different phases in the ap- 
proach to equilibrium described by Kubo-Martin-Schwinger relations for thermal 
equilibrium states. A detailed comparison of the full quantum dynamics to ap- 
proximate schemes like that of a standard kinetic (on-shell) Boltzmann equation 
is performed. We find that far off-shell 1 <-> 3 processes are responsible for chem- 
ical equilibration, which is not included in the Boltzmann limit. Furthermore, 
we derive generalized transport equations for the same theory in a first order 
gradient expansion in phase space thus explicitly retaining the off-shell dynam- 
ics as inherent in the time-dependent spectral functions. The solutions of these 
equations compare very well with the exact solutions of the full Kadanoff-Baym 
equations with respect to the occupation numbers of the individual modes, the 
spectral evolution as well as the chemical equilibration process. Furthermore, the 
proper equilibrium off-shell distribution is reached for large times contrary to the 
quasiparticle Boltzmann limit. We additionally present a direct comparison of 
the solution of the generalized transport equations in the Kadanoff-Baym (KB) 
and Botermans-Malfliet (BM) form; both solutions are found to agree very well 
with each other. The off-shell transport equation in the BM scheme allows for an 
explicit solution within an extended dynamical quasiparticle Ansatz. This leads 
to generalized equations of motion for dynamical quasiparticles that exceed the 
classical Hamilton equations and allow for the description of dynamical spectral 
functions. The dynamical quasiparticle model (DQPM) is used to extract par- 
tonic spectral functions from lattice QCD in the temperature range 0.8 < T/T c < 
10. By consideration of time-like and space-like sectors of 'observables' such as 
number densities, energy densities etc. mean-field potentials as well as effective in- 
teractions are extracted at different temperature T and quark chemical potentials 
\i q . The latter determine the off-shell dynamics in the Parton-Hadron-String Dy- 
namics (PHSD) transport approach. Illustrative examples for off-shell dynamics 
in the hadron and parton sector are presented for e + e~ production from nucleus- 
nucleus collisions from SIS to RHIC energies. The generalized transport approach, 
furthermore, qualifies for the description of hadronization without leading to the 
problem of entropy reduction as in conventional coalescence models. 
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Introduction 

Non-equilibrium many-body theory or quantum field theory has become a major topic of re- 
search for transport processes in nuclear physics, in cosmological particle physics as well as 
condensed matter physics. The multidisciplinary aspect arises due to a common interest to 
understand the various relaxation phenomena of quantum dissipative systems. Recent progress 
in cosmological observations has also intensified the research on quantum fields out of equi- 
librium. Important questions in high-energy nuclear or particle phyics at the highest energy 
densities are: i) how do noncquilibrium systems in extreme environments evolve, ii) how do they 
eventually thermalize, iii) how phase transitions do occur in real time with possibly nonequilib- 
rium remnants, and iv) how do such systems evolve for unprecedented short and nonadiabatic 
timescales? 

The very early history of the universe provides important scenarios, where non-equilibrium 
effects might have played an important role, like in the (post-) inflationary epoque (see e.g. [1- 
4]), for the understanding of baryogenesis (see e.g. [1,5]) and also for the general phenomena of 
cosmological decoherence [6] . Referring to modern nuclear physics the dynammics of heavy- ion 
collisions at various bombarding energies has always been a major motivation for research on 
nonequilbrium quantum many-body physics and relativistic quantum-field theories, since the 
initial state of a collision resembles an extreme nonequilibrium situation while the final state 
might even exhibit a certain degree of thermalization. Indeed, at the presently highest energy 
heavy- ion collider experiments at RHIC, where one expects to create experimentally a transient 
deconfined state of matter denoted as quark-gluon plasma (QGP) [7], there are experimental 
indications - like the build up of collective flow - for an early thermalization accompanied with 
the build up of a very large pressure. These examples demonstrate that one needs an ab-initio 
understanding of the dynamics of out-of-equilibrium quantum-field theory. 

Especially the powerful method of the 'Schwingcr-Keldysh' [8-12] or 'closed time path' 
(CTP) (non-equilibrium) real-time Greens functions has been shown to provide an appropri- 
ate basis for the formulation of the complex problems in the various areas of nonequilibrium 
quantum many-body physics. Within this framework one can derive valid approximations - 
depending, of course, on the problem under consideration - by preserving overall consistency 
relations. Originally, the resulting causal Dyson-Schwinger equation of motion for the one- 
particle Greens functions (or two-point functions), i.e. the Kadanoff-Baym (KB) equations [13], 
have served as the underlying scheme for deriving various transport phenomena and generalized 
transport equations. For review articles on the Kadanoff-Baym equations in the various areas 
of noncquilibrium quantum physics we refer the reader to Refs. [14-19]. We note in passing, 
that also the 'influence functional formalism' has been shown to be directly related to the KB 
equations [20]. Such a relation allows to address inherent stochastic aspects of the latter and 
also to provide a rather intuitive interpretation of the various self-energy parts that enter the 
KB equations. 

Furthermore, kinetic transport theory is a convenient tool to study many-body nonequilib- 
rium systems, nonrelativistic or relativistic. Kinetic equations, which do play the central role 
in more or less all practical simulations, can be derived by means of appropriate KB equa- 
tions within suitable approximations. Hence, a major impetus in the past has been to derive 
semi-classical Boltzmann-like transport equations within the standard quasi-particle approxi- 
mation. Additionally, off-shell extensions by means of a gradient expansion in the space-time 
inhomogenities - as already introduced by Kadanoff and Baym [13] - have been formulated: for 
an relativistic electron-photon plasma [21], for transport of electrons in a metal with external 
electrical field [22,23], for transport of nucleons at intermediate heavy-ion reactions [24], for 
transport of particles in ^ 4 -theory [18,25], for transport of electrons in semiconductors [19,26], 
for transport of partons or fields in high-energy heavy-ion reactions [30-33], or for a trapped 
Bose system described by effective Hartree-Fock-Bogolyubov kinetic equations [34]. We recall 
that on the formal level of the KB-equations the various forms assumed for the self-energy have 
to fulfill consistency relations in order to preserve symmetries of the fundamental Lagrangian 
[13,35,36]. This allows also for a unified treatment of stable and unstable (resonance) particles. 

In non-equilibrium quantum-field theory typically the nonperturbative description of (second- 
order) phase transitions has been in the foreground of interest by means of mean-field (Hartree) 
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descriptions [2,8,37-39], with applications for the evolution of disoriented chiral condensates 
or the decay of the (oscillating) inflaton in the early reheating era of the universe. 'Effective' 
mean-field dissipation (and decoherence) - solving the so called 'backreaction' problem - was 
incorporated by particle production through order parameters explicitly varying in time. How- 
ever, it had been then realized that such a dissipation mechanism, i.e. transferring collective 
energy from the time-dependent order parameter to particle degrees of freedom, can not lead to 
true dissipation and thermalization. Such a conclusion has already been known for quite some 
time within the effective description of heavy-ion collisions at low energy. Full time-dependent 
Hartrec or Hartree-Fock descriptions [40] were insufficient to describe the reactions with in- 
creasing collision energy; additional Boltzmann-like collision terms had to be incorporated in 
order to provide a more adequate description of the collision processes. 

The incorporation of true collisions then has been formulated also for various quantum-field 
theories [41-46]. Here, a systematic 1/N expansion of the "2PI effective action" is conventionally 
invoked [18,37,46] serving as a nonperturbative expansion parameter. Of course, only for large 
TV this might be a controlled expansion. In any case, the understanding and the influence of 
dissipation with the chance for true thermalization - by incorporating collisions - has become a 
major focus of recent investigations. The resulting equations of motion always do resemble the 
KB equations; in their general form (beyond the mean field or Hartree(-Fock) approximation) 
they do break time invariance and thus lead to irreversibility. This macroscopic irreversibility 
arises from the truncations of the full theory to obtain the self-energy operators in a specific 
limit. As an example we mention the truncation of the (exact) Martin-Schwinger hierarchy in 
the derivation of the collisional operator in Ref. [24] or the truncation of the (exact) BBGKY 
hierarchy in terms of n-point functions [24,41,42,47]. 

In principle, the non-equilibrium quantum dynamics is nonperturbative in nature. Unphys- 
ical singularities only appear in a limited truncation scheme, e.g. ill-defined pinch singularities 
[48] , which do arise at higher order in a perturbative expansion in out of equilibrium quantum- 
field theory, are regularized by a consistent nonperturbative description (of Schwinger-Dyson 
type) of the non-equilibrium evolution, since the resummed propagators obtain a finite width 
[49]. Such a regularization is also observed by other resummation schemes like the dynamical 
renormalization-group technique [45,50]. 

Although the analogy of KB-type equations to a Boltzmann-like process is quite obvious, this 
analogy is far from being trivial. The full quantum formulation contains much more information 
than a semi-classical (generally) on-shell Boltzmann equation. The dynamics of the spectral 
(i.e. 'off-shell') information is fully incorporated in the quantum dynamics while it is missing 
in the Bolzmann limit. A full answer to the question of quantum equilibration can thus only 
be obtained by detailed numerical solutions of the quantum description itself. 

We briefly address previous works that have investigated numerically approximate or full 
solutions of KB-type equations. A seminal work has been carried out by Danielewicz [51], who 
investigated for the first time the full KB equations for a spatially homogenous system with a 
deformed Fermi sphere in momentum space for the initial distribution of occupied momentum 
states in order to model the initial condition of a heavy-ion collision in the nonrelativistic do- 
main. In comparison to a standard on-shell semi-classical Boltzmann equation the full quantum 
Boltzmann equation showed quantitative differences, i.e. a larger collective relaxation time for 
complete equilibration of the momentum distribution /(p, t). This slowing down of the quan- 
tum dynamics was attributed to quantum interference and off-shell effects. Similar quantum 
modifications in the equilibration and momentum relaxation have been found for a relativis- 
tic situation in Ref. [52]. In the following, full and more detailed solutions of nonrelativistic 
KB equations have been performed by Kohler [53,54] with special emphasis on the build up 
of initial many-body correlations on short time scales. The role of memory effects has been 
clearly shown experimentally by femtosecond laser spectroscopy in semiconductors [55] in the 
relaxation of excitons. Solutions of quantum transport equations for semiconductors [19,56] - 
to explore relaxation phenomena on short time and distance scales - has become also a very 
active field of research [22,57]. 

In the last decade the numerical treatment of general out-of-equilibrium quantum dynamics, 
as described within 2PI effective approaches at higher order, have become more frequent. The 
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subsequent equations of motion are very similar to the KB equations, although more involved 
expressions for a non-local vertex function - as obtained by the 2PI scheme - are incorporated. 
The numerical solutions are also obtained only for homogenous systems, so far (cf. [58-63]). 

Apart from complete thcrmalization of all single-particle degrees of freedom the quantum 
dynamics of the spectral function is also a lively discussed issue in the microscopic modeling of 
hadronic resonances with a broad mass distribution. This is of particular relevance for simula- 
tions of heavy-ion reactions, where e.g. the Z\-resonance or the p-meson already show a large 
decay witdh in vacuum [64-69]. Especially the p vector meson represents the most promising 
hadronic particle for showing possible in-medium modifications in hot and compressed nuclear 
matter since the leptonic decay products are of only weakly interacting electromagnetic na- 
ture. Hence a consistent formulation for the transport of extremely short-lived particles beyond 
the standard quasi-particle approximation is needed. On the one side, there exist purely for- 
mal developments starting from a first-order gradient expansion of the underlying KB equations 
[67-69] , while on the other side already practical realisations for various questions have emerged 
[64-66,70]. The general idea is to obtain a description for the propagation of dynamical spectral 
functions, i.e. a propagation in the off-shell mass squared M 2 . 

The aim of these lectures is to provide a short survey on (nonperturbative) nonequilibrium 
quantum field theory with actual applications (and numerical illustrations) for the model case 
of scalar ^ 4 -theory. This includes the derivation of off-shell as well as on-shcll (Boltzmann) 
transport equations and their comparison to the exact solution for a variety of out-of-equilibrium 
initial conditions. Furthermore, a generalized testpartiele Ansatz is introduced that allows for a 
convenient solution of off-shell transport equations which involve only the complex (retarded) 
selfenergies for the relevant degrees of freedom. In case of hadrons these selfenergies may be 
determined from effective Lagrangian models whereas in case of partons suitable selfenergies 
can be fixed by the dynamical quasiparticlc model in comparison to lattice QCD at finite 
temperature T. Applications for dilepton production in nucleus-nucleus collisions provide a 
practical case for the spectral off-shell dynamics and bridge the gap to current spectra taken 
by various experimcnal collaborations [71]. On the partonic side we present an illustrative 
example for the build up of collective flow, strangeness equilibration as well as the dynamics of 
hadronization. 



1 Quantum field dynamics and thermalization 

1.1 Non-equilibrium dynamics 

1.1.1 The Kadanoff-Baym equations 

As mentioned in the Introduction, a natural starting point for non-equilibrium theory is pro- 
vided by the closed-time-path (CTP) method. Here all quantities are given on a special real-time 
contour with the time argument running from — oo to oo on the chronological branch (+) and 
returning from oo to — oo on the antichronological branch (— ). In cases of systems prepared at 
a time to this value is (instead of — oo) the start and end point of the real-time contour (cf. 
Fig. 1). In particular the path ordered Green functions (in case of real scalar fields <f>(x)) are 
defined as 

G(x,y) = (T p {^x)^y)}) (1) 

-e P (xo-2/o)(^(x)0(y)) + e p (yo-xo)(<f>(y)4>(x)) 

where the operator T p orders the field operators according to the position of their arguments 
on the real-time path as accomplished by the path step- functions O p . The expectation value in 
(2) is taken with respect to some given density matrix po, which is constant in time, while the 
operators in the Heisenberg picture contain the entire time dependence of the non-equilibrium 
system, i.e. 0(f) = exp(iH(t — to)) O cxp(—iH(t — to))- 
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Fig. 1. The closed time contour in the Schwinger-Keldysh formalism 



Self-consistent equations of motion for these Green functions can be obtained with help of 
the two-particle irreducible (2PI) effective action r[G]. It is given by the Legendre transform 
of the generating functional of the connected Green functions W as 

r[G] = r° + l - [ln(l-Q p G ® p E) + QpGQpS] + <P[G] (2) 

in case of vanishing vacuum expectation value <O|0(x)|O>= 0. In (2) _T° depends only on 
free Green functions and is treated as a constant, while the symbols P represent convolution 
integrals over the closed time path in Fig. 1. The functional <l> is the sum of all closed 2PI 
diagrams built up by full propagators G; it determines the self-energies by functional variation 
as 

^- 2 'SG^)- (3) 

From the effective action r(G) the equations of motion for the Green function are determined 
by the stationarity condition 

5r/SG = 0. (4) 



1.1.2 <2> 4 -theory 

The scalar (/> 4 -theory is an example for a fully relativistic field theory of interacting scalar 
particles that allows to test theoretical approximations [41,42] without coming to the problems 
of gauge-invariant truncation schemes as encountered for gauge (vector) fields [41,58-60]. Its 
Lagrangian density is given by (x = (i,x)) 

C(x) = l -d,4>(x)d^{x) - \m 2 4> 2 {x) - ±4>\x) , (5) 

where m denotes the 'bare' mass and A is the coupling constant determining the interaction 
strength of the scalar fields. 

In the present calculation we take into account contributions up to the 3-loop order for the 
^-functional (cf. Fig. 2) which reads explicitly 

i$ = jj^x G{x,x) 2 - ^ J^d d+1 x J^d d+1 y G{x,y)\ (6) 

where d denotes the spatial dimension of the problem. 

This approximation corresponds to a weak coupling expansion such that we consider con- 
tributions up to the second superficial order in the coupling constant A (cf. Fig. 3). For the 
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Fig. 2. Contributions to the ^-functional for the Kadanoff-Baym equation: two-loop contribution 
(l.h.s.) giving the tadpole self-energy and three-loop contribution (r.h.s.) generating the sunset self- 
energy. The ^-functional is built-up by full Green functions (double lines) while open dots symbolize 
the integration over the inner coordinates. 



superficial coupling constant order we count the explicit coupling factors A associated with the 
visible vertices. The hidden dependence on the coupling strength - which is implicitly incorpo- 
rated in the self-consistent Green-functions that built up the ^-functional and the self-energies 
- is ignored on that level. For our present purpose this approximation is sufficient since we 
include the leading mean-field effects as well as the leading order scattering processes that pave 
the way to thermalization. 

For the actual calculation it is advantageous to change to a single-time representation for 
the Green functions and self-energies defined on the closed-time-path. In line with the position 
of the coordinates on the contour there exist four different two-point functions 

i G c (x, y) = i G ++ (x, y) = < T c { cj>{x) cj>{y) } > , (7) 

iG<{x,y)=iG + -(x,y) = < { 4>{y) cj>{x) } > , 

iG > (x,y) = iG-+(x,y) = < { c/>(x) cf>(y) } > , 

iG a (x,y) = iG—(x,y) = < T a { cj>(x) ^y) } > . 

Here T c (T a ) represent the (anti-)time-ordering operators in case of both arguments lying on the 
(anti-)chronological branch of the real-time contour. These four functions are not independent of 
each other. In particular the non-continuous functions G c and G a are built up by the Wightman 
functions G > and G < and the usual ©-functions in the time coordinates. Since for the real 
boson theory (5) the relation G > (x, y) — G K (y, x) holds, the knowledge of the Green functions 
G < (x, y) for all x, y characterizes the system completely. Nevertheless, we will give the equations 
for G < and G > explicitly since this is the familiar representation for general field theories [25]. 

By using the stationarity condition for the action (4) and resolving the time structure of 
the path ordered quantities we obtain the Kadanoff-Baym equations for the time evolution of 
the Wightman functions [58,72]: 

- [d^+m 2 ] G*{x,y) - E 5 (x) G^(x,y) (8) 

+ ^dz Q Jd d z [S>(x,z)- E < (x,z)] G^(z,y) 

dz d d z S>(x,z) [G > (z,y)-G < (z,y)] 1 

t J 
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[dld^+m 2 ] G*(x,y) = E s {y) G*(x,y) 



+ / dz d d z [G > (x,z)-G < (x,z)] S>(z,y) 



dz d d z G > (x,z) [E > (z,y)-S < (z,y)], 





Fig. 3. Self-energies of the Kadanoff-Baym equation: tadpole self-energy (l.h.s.) and sunset self-energy 
(r.h.s.). Since the lines represent full Green functions the self-energies are self-consistent (see text) with 
the external coordinates indicated by full dots. 



Within the 3- loop approximation for the 2PI effective action (i.e. the ^-functional (6)) we get 
two different self-energies: In leading order of the coupling constant only the tadpole diagram 
(l.h.s. of Fig. 3) contributes and leads to the generation of an effective mass for the field quanta. 
This self-energy (in coordinate space) is given by 



S\x) = ~iG<(x,x) , 



(9) 



and is local in space and time. In next order in the coupling constant (i.e. A 2 ) the non-local 
sunset self-energy (r.h.s. of Fig. 3) enters the time evolution as 



E*{x,y) = ~G*(x,y)G*(x,y)G*(y,x) 
6 



E*(x,y) 



\2 r < -i 3 

y G>(x,y) 



(10) 
(11) 



Thus the Kadanoff-Baym equation (8) in our case includes the influence of a mean-field on the 
particle propagation - generated by the tadpole diagram - as well as scattering processes as 
inherent in the sunset diagram. 

The Kadanoff-Baym equation describes the full quantum nonequilibrium time evolution on 
the two-point level for a system prepared at an initial time i.e. when higher order correlations 
are discarded. The causal structure of this initial value problem is obvious since the time 
integrations are performed over the past up to the actual time xq (or yo, respectively) and do 
not extend to the future. 

Furthermore, also linear combinations of the Green functions in single-time representation 
are of interest and will be exploited for the spectral properties of the system later on. The 
retarded Green function G R and the advanced Green function G A are given as 



G R (x u x 2 ) = ©(ti-ta) [G > (x 1 ,x 2 )-G < (x u x 2 )] 



(12) 
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= G c (x 1 ,x 2 ) - G < (xi,x 2 ) = G > (xi,x 2 ) - G a (x u x 2 ) , 

G A (x 1 ,x 2 ) = -G(t 2 -t 1 ) [G > (x 1 ,x 2 )-G < (x u x 2 )] (13) 
= -G(t 2 -t 1 )([<j>(x 1 ),<f>(x 2 )]_) 

= G c ( Xl ,x 2 ) - G>( Xl ,x 2 ) = G < (x 1 ,x 2 ) - G a ( Xl ,x 2 ). 

These Green functions contain exclusively spectral, but no statistical information of the system. 
Their time evolution is determined by Dyson-Schwinger equations and given by (cf. Ref. [72]) 

-ffid£ 1 + m 2 + E s (x 1 )]G R {x u x 2 ) = (14) 



S (d+1 Hx!-x 2 ) + Jd d+1 z E R {x u z) G R (z,x 2 



-[d^d^ + m 2 + E s ( Xl )] G A (xi,x 2 ) = (15) 



5 {d+1 \ Xl -x 2 ) + jd d+1 z S A {x u z) G A {z,x 2 ) 



where the retarded and advanced self-energies S R , S A are defined via S > , S < similar to the 
Green functions. Thus the retarded (advanced) Green functions are determined by retarded 
(advanced) quantities, only. 



1.1.3 Homogeneous systems in space 

In the following we will consider homogeneous systems in space. To obtain a numerical solution 
the Kadanoff-Baym equation (8) is transformed to momentum space: 

d 2 1 G < (p,t 1 ,t 2 ) = -{p 2 +m 2 + £ s (t 1 )]G < (p,t u t 2 ) (16) 



- [df [E>(p,t u t')-E<(p,t 1 ,t')] G<(p,t',t 2 ) 

Jt 

+ [ t2 df E<(p,h,t') [G>(p,t',t 2 )-G<(p,t',t 2 )] 

Jta 



'to 

= -[p 2 + m 2 + Z s (t 1 )]G < (p,t u t 2 ) + I<(p,h,t 2 ) 

where both memory integrals are summarized in the function 1^. The equation of motion in 
the second time direction t 2 is given analogously. In two-time and momentum space (p, t, t') 
representation the self-energies read 



E>(p,t,t') = -- J-^J-^ G>(q,t,«') G>(r,«,f) G<(q+r-p,f,t) . 
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1.2 Numerical studies on equilibration 

1.2.1 Numerical implementation 

For the solution of the Kadanoff-Baym equations (16) a flexible and accurate algorithm has 
been developed: Instead of solving the second order differential equation (16) one can generate 
a set of first order differential equations for the Green functions in the Heisenberg picture, 

iG^ (t> (x 1 ,x 2 ) = < <p{x 2 )(t>{xi) > = iG < (x 1 ,x 2 ) , (18) 
iG^{x 1 ,x 2 ) = < 4>{x2)n(x!) > = 9 tl iG^(n,a;2) , 
iG^ n (x 1 ,x 2 ) = < n(x 2 )<t>(xi) > = 9 t2 iG^(xi,x 2 ) , 

iG< n (x!,X2) = < tt(x 2 )tt(x 1 ) > = d tl d t2 i G^ (j> (x 1 , x 2 ) , 

with the canonical field momentum ir(x) = d Xo <j>(x). The first index ir or <j> is always related to 
the first space-time argument. Exploiting the time-reflection symmetry of the Green functions 
some of the differential equations are redundant. The required equations of motion are given as 

d tl G<^(p,t 1 ,t 2 ) = G< <p (p,t 1 ,t 2 ), (19) 
di G&(p, i, t) = 2ilm{ G< (p, i, t) } , 

d tl G<^v,hM) = -fl 2 {ti) G^(p,ti,t 2 ) + IHvMM) , 
dt 2 G<^p,h,t 2 ) = G< v (p, t u t 2 ), 

dtG< 4> (p,i,t) = -n 2 (i)G< 4> (p,t,i) + G<„(p,i,t) + iHp,t,t), 

d tl G< v {p,t 1 ,t 2 ) = -n 2 ^) g^(pmm) + i^ 2 (p,t u t 2 ) , 

dtG< w (p,t,t) = - f2 2 {t) 2i Im{G<^p,i,t}} + 2iIm{I< 2 (p,i,F)} , 

where t — (t\ + t 2 )/2 is the mean time variable. Thus one explicitly considers the propagation 
in the time diagonal (cf. Ref. [53]). In the equations of motion (19) the current (renormalized) 
effective energy including the time dependent tadpole contribution enters, 

Q 2 (t) = p 2 + m 2 + Sm 2 tad + Sm 2 sun + Z 5 (t), (20) 

with Sm 2 ad and 5m 2 un specified in the next Subsection. The evolution in the t 2 direction has 
not be taken into account for G^ and G<^ since the Green functions beyond the time diagonal 
{t 2 > t\) are determined via the time reflection symmetry G < :(p,t\,t 2 ) = — [G < :(p, t 2 , t\) ]* 
from the known values for the lower time triangle in both cases. Since there is no time reflection 
symmetry for the G^^ functions, they have to be calculated (and stored) in the whole t\, 
t 2 range. However, we can ignore the evolution of G^ since it is obtained by the relation 
G^{p 1 h, t 2 ) = — [G^(p, t 2 , ti) }*. The correlation integrals in (19) are given by 

I<(p,t u t 2 ) = -f^dt' [E>(p,t 1 ,t')-E<(p,t 1 ,t')] G< (p,t',t 2 ) (21) 
Jo 

+f 2 dt' S<(p,t u t') [G< (-p,t 2 ,t')-G< ( P ,t',i 2 ) 

Jo 1 
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/i< 2 (p,*i,*2) = d t2 If{p,t u t 2 ) (22) 

V [S > (p,t 1 ,t')-E<(p,t u t')] G< n (p,t',t 2 ) 

G<^-p,t 2 ,t')-G<„(p,t>,t 2 ) . 

In (19) and (22) one can replace G^ n (p,ti,t 2 ) = — [G^(p, t 2 , t\) ]* such that the set of equa- 
tions is closed in the Green functions G^, G^ and G< v . 

The disadvantage, to integrate more Green functions in time in this first-order scheme, is 
compensated by its good accuracy. As mentioned before, we especially take into account the 
propagation along the time diagonal which leads to an improved numerical precision. The set 
of differential equations (19) is solved by means of a 4th order Runge-Kutta algorithm. For 
the calculation of the self energies a Fourier method similar to Refs. [15,53] is applied. The 
self energies (17), furthermore, are calculated in coordinate space where they are products of 
coordinate-space Green functions (that are available by Fourier transformation) and finally 
transformed to momentum space. 




+/ dt' E < (p,t 1 ,t') 



1.2.2 Renormalization of </> 4 -theory in 2+1 dimensions 

In 2+1 space-time dimensions both self-energies (9), (10) incorporated in the present case are 
ultraviolet divergent. Since the particles have a finite mass no problems arise from the infrared 
momentum regime. The ultraviolet regime, however, has to be treated explicitly. 

For the renormalization of the divergences we only assume that the time-dependent nonequi- 
librium distribution functions are decreasing for large momenta comparable to the equilibrium 
distribution functions, i.e exponentially. Thus one can apply the conventional finite temperature 
renormalization scheme. By separating the real time (equilibrium) Green functions into vacuum 
(T = 0) and thermal parts it becomes apparent, that only the pure vacuum contributions of 
the self-energies are divergent. For the linear divergent tadpole diagram a mass counter term 
(at the renormalized mass m) can be introduced as 

/d^p 1 
(2^2^' u; p = Vp 2 +m 2 , (23) 

which cancels the contribution from the momentum integration of the vacuum part of the Green 
function. 

In case of the sunset diagram only the logarithmically divergent pure vacuum part requires 
a renormalization, while it remains finite as long as at least one temperature line is involved. 
Contrary to the case of 3+1 dimensions it is not necessary to employ the involved techniques 
developed for the renormalization of self-consistent theories (in equilibrium) in Refs. [73] be- 
cause in 2+1 dimensions <P 4 theory is superrenormalizable. Since the divergence only appears 
(in energy-momentum space) in the real part of the Fcynman self-energy S c at T = (and 
equivalently in the real part of the retarded/advanced self-energies Jj ret / adv ^ it can be absorbed 
by another mass counterterm 

dm 2 sun = -Re S C T=0 (P 2 ) = -Re S r T % adv (p 2 ) (24) 

_ A 2 f d\ f cPr 1 w q +w r +w q+r _p 

6 7(27r) 2 J (2ir) 2 4 w q w r w q+r _ p [ u q +w r +cj q+r _ p ] 2 - pi 

at given four-momentum p = (p , p) and renormalized mass to. 

In Fig. 4 we demonstrate the applicability of the renormalization prescription. To this aim 
we display two momentum modes | p |/m = 0.0 (upper plots) and | p |/m = 2.0 (lower plots) of 
the equal-time Green function zG < (|p|, t, t) for various momentum cut-offs Pmax/m = 6 (2) 20 
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with renormalization 



without renormalization 



a 

V 



a 

V 




time t m 



time t m 



Fig. 4. Time evolution of two momentum modes | p |/m = 0.0, | p |/m = 2.0 of the equal-time Green 
function starting from the initial distribution D2 (as specified in Section 1.2.3) with coupling constant 
A/m = 14. One observes, that with the renormalization of the sunset diagram (left plots) a proper limit 
is obtained when increasing the momentum cut-off p ma x/m = 6(2)24, while without renormalization 
(right plots) the curves tend to infinity when increasing the ultraviolet cut-off 



with (left plots) and without (right plots) renormalization of the sunset self-energy. For both 
cases the renormalization of the tadpole diagram has been used. We mention, that a non- 
renormalization of the tadpole self-energy has even more drastic consequences in accordance 
with the linear degree of divergence. For the non-renormalized calculations - with respect to 
the sunset diagram - we observe that both momentum modes do not converge with increasing 
momentum space cut-off. In fact, all lines tend to infinity when the maximum momentum is 
enlarged (since the gridsize of the momentum grid is kept constant). Although the divergence 
as a function of the momentum cut-off is rather weak - in accordance with the logarithmic 
divergence of the sunset self-energy in 2+1 space-time dimensions - a proper ultraviolet limit 
is not obtained. 

This problem is cured by the sunset mass counter term (25) as seen on the left side of Figure 
4. For the momentum mode |p|/m = 2.0 the calculations converge to a limiting curve with 
increasing momentum cut-off. Even for the more selective case of the |p|/to = 0.0 mode of 
the equal-time Green function the convergence is established. We point out that this limit is 
obtained for the unequal-time Green functions as well (not shown here explicitly). In fact, it 
turns out that the equal-time functions provide the most crucial test for the applicability of 
the renormalization prescription, since the divergent behaviour appears to be less pronounced 
for the propagation along a single time direction t\ or t%. Thus we can conclude that the 
renormalization scheme introduced above, i.e. including mass counter terms for the divergent 
tadpole and sunset self-energies, leads to ultraviolet stable results. 
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In the following Subsections we will use the 'bare' mass m = 1, which implies that times 
are given in units of the inverse mass or t ■ m is dimcnsionlcss. Accordingly, the bare coupling 
A in (7) is given in units of the mass m such that X/m is dimensionless, too. 

1.2.3 Initial conditions 

In order to investigate equilibration phenomena on the basis of the Kadanoff-Baym equations for 
the 2+1 dimensional problem, one first has to specify the initial conditions for the time integra- 
tion. To this aim we consider four different initial distributions that are all characterized by the 
same energy density. Consequently, for large times (— > oo) all initial value problems should lead 
to the same equilibrium final state. The initial equal-time Green functions iG < (p, t = 0, t = 0) 
adopted are displayed in Fig. 5 (l.h.s.) as a function of the momentum p z . We first concentrate 
on polar symmetric configurations due to the large numerical expense for this first investiga- 
tion 1 . Since the equal-time Green functions G < (p, t, t, ) are purely imaginary we show only the 
real part of i G < in Fig. 5. 




-12 -10 -8 -6 -4 -2 2 4 6 8 10 12 -12 -10 -8 -6 -4 -2 2 4 6 8 10 12 

momentum p /m momentum p / m 



Fig. 5. Initial Green functions iG < (|p|,t = 0, t = 0) (l.h.s.) and corresponding initial distribution 
functions n(| p |, t = 0) (r.h.s.) for the distributions Dl, D2, D3 and DT in momentum space (for a cut 
of the polar symmetric distribution in p x direction for p y = 0). 

Furthermore, the corresponding initial distribution functions in the occupation density 
n(p, t = 0), related to iG < (p, t = 0, t = 0) via 

2u p iG < (p, t = 0, t = 0) = 2n(p, t = 0) + 1, (25) 

are shown in Fig. 5 on the r.h.s. While the initial distributions Dl, D2, D3 have the shape 
of (polar symmetric) 'tsunami' waves with maxima at different momenta in p z , the initial 
distribution DT corresponds to a free Bose gas at a given initial temperature To that is fixed by 
the initial energy density. According to (25) the difference between the Green functions and the 
distribution functions is basically given by the vacuum contribution, which has its maximum at 
small momenta. Thus even for the distributions Dl, D2, D3 the corresponding Green functions 
are non- vanishing for p as 0. 

Since we consider a finite volume V = a 2 we work in a basis of momentum modes charac- 
terized by the number of nodes in each direction. The number of momentum modes is typically 
in the order of 40 which is found to be sufficient for numerically stable results. For times t < 
we consider the systems to be noninteracting and switch on the interaction (~ A) for t=0 to 
explore the quantum dynamics of the interacting system for t > 0. 



1 In Section 1.4 we will present also calculations for non-symmetric systems. 
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1.2.4 Equilibration in momentum space 

The time evolution of various (selected) momentum modes of the equal-time Green function for 
the different initial states Dl, D2, D3 and DT is shown in Fig. 6, where the dimensionless time 
t ■ m is displayed on a logarithmic scale. We observe that starting from very different initial 
conditions - as introduced above - the single momentum modes converge to the same respective 
occupation numbers for large times as characteristic for a system in equilibrium. As noted above, 
the initial energy density is the same for all distributions and energy conservation is fulfilled 
strictly in the time integration of the Kadanoff-Baym equations. The different momentum modes 
in Fig. 6 typically show a three-phase structure: For small times (t ■ m < 10) one finds damped 
oscillations that can be identified with a typical switching-on effect at t=0, where the system 
is excited by a sudden increase of the coupling constant to A/to = 18. The damping of the 
initial oscillations depends on the coupling strength A/m and is more pronounced for strongly 
coupled systems. We mention that for very strong interactions A/m > 30 the initial oscillations 
are even hard to recognize. 




time t m 



Fig. 6. Time evolution of selected momentum modes of the equal-time Green function |p|/m — 
0.0,0.8, 1.6,2.4,3.2 (from top to bottom) for four different initial configurations Dl, D2, D3, and DT 
(characterized by the different line types) with the same energy density. For the rather strong coupling 
constant A/m = 18 the initial oscillations - from switching on the interaction at t=0 - are damped 
rapidly and disappear for t ■ m > 20. Finally, all momentum modes assume the same respective equi- 
librium value for long times t ■ m > 500) independent of the initial state. 

For 'intermediate' time scales (10 < t ■ to < 500) one observes a strong change of all 
momentum modes in the direction of the final stationary state. We address this phase to 'kinetic' 
equilibration and point out, that - depending on the initial conditions and the coupling strength 

- the momentum modes can temporarily even exceed their respective equilibrium value. This 
can be seen explicitly for the lowest momentum modes (p/m — or = 0.8) of the distribution 
Dl in Fig. 6, which possesses initially maxima at small momentum. Especially the momentum 
mode |p| = of the equal-time Green function G < , which starts at around 0.52, is rising 
to a value of ~ 1 .95 before decreasing again to its equilibrium value of ~ 1.29. Thus the time 
evolution towards the final equilibrium value is - after an initial phase with damped oscillations 

- not necessarily monotonic. For different initial conditions this behaviour may be weakened 
significantly as seen for example in case of the initial distribution D2 in Figure 6. Coincidently, 
both calculations Dl and D2 show approximately the same equal-time Green function values 
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for times t • m > 80. Note, that for the initial distribution D3 the non-monotonic behaviour is 
not seen any more. 

In general, one observes that only initial distributions (of the well type) show this feature 
during their time evolution, if the maximum is located at sufficiently small momenta. Initial 
configurations like the distribution DT - where the system at t = is given by a free gas 
of particles at a temperature Tq - do not show this property. Although the DT distribution 
is not the equilibrium state of the interacting theory, the actual numbers are much closer to 
the equilibrium state of the interacting system than the initial distributions Dl, D2 and D3. 
Therefore, the evolution for DT proceeds less violently. We point out, that in contrast to the 
calculations performed for </> 4 -theory in 1+1 space-time dimensions [60] we find no power law 
behaviour for intermediate time scales. 

The third phase, i.e. the late time evolution (t ■ m > 500) is characterized by a smooth 
approach of the single momentum modes to their respective equilibrium values. As we will see 
in Section 1.4 this phase is adequately characterized by chemical equilibration processes. 

1.3 The different phases of quantum equilibration 

1.3.1 Build-up of initial correlations 
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Fig. 7. Change of the different contributions to the total energy density in time. The sunset energy 
density e sun decreases rapidly in time; this contribution is approximately compensated by an increase 
of the kinetic energy density Skin- Together with the smaller tadpole contribution etad the total energy 
Etot is conserved. 



The time evolution of the interacting system within the Kadanoff-Baym equations is charac- 
terized by the build-up of early correlations. This can be seen from Fig. 7 where all contributions 
to the energy-density are displayed separately as a function of time and the total energy density 
£tot(t = 0) is subtracted. The kinetic energy density s^in is represented by all parts of etot that 
are independent of the coupling constant (oc A ). All terms proportional to A 1 are summarized 
by the tadpole energy density etad including the actual tadpole term as well as the correspond- 
ing tadpole mass counterterm (cf. Section 1.2.2). The contributions from the sunset diagram 
(oc A 2 ) are represented by the sunset energy density e sun . 



£tot(t) = ekin(t) + e ta d{t) + e sun (t) , 



(26) 
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£kin{t)= \ j-j^p (p 2 +m 2 o) %G^(p,t,t) + \ j-j^p iG< w (p,t,t), 

1 f d d p 1 f d d p 

etad(t) = - J -^-^ Stad{t) iG^(p,t,t) + - J j^-y Sm 2 tad iG^(p,t,t) , 

1 f d d p 1 f d d p 

£sun{t) = — A Jj^yi Hf(p,t,t) + - Jj-^ 5m sun »G< (p,t,t). 

The calculation in Fig. 7 has been performed for the initial distribution DT (which represents 
a free gas of Bose particles at temperature To w 1.736 m) with a coupling constant of X/m = 18. 
This state is stationary in the well-known Boltzmann limit (cf. Section 1.4), but it is not for 
the Kadanoff-Baym equation. In the full quantum calculations the system evolves from the 
uncorrelated initial state and the correlation energy density s sun decreases rapidly with time. 
The decrease of the correlation energy e sun which is - with exception of the sunset mass 
counterterm contribution - initially zero is approximately compensated by an increase of the 
kinetic energy density £ki n - Since the kinetic energy increases in the initial phase, the final 
temperature Tf is slightly higher than the initial 'temperature' To. The remaining difference is 
compensated by the tadpole energy density etad such that the total energy density is conserved. 

While the sunset energy density and the kinetic energy density always show a time evolution 
comparable to Fig. 7 the change of the tadpole energy density depends on the initial config- 
uration and may be positive as well. Since the self energies are obtained within a <?-derivable 
scheme the fundamental conservation laws, as e.g. energy conservation, are respected to all 
orders in the coupling constant. When neglecting the cx A 2 sunset contributions and starting 
with a non-static initial state of identical energy density one observes the same compensating 
behaviour between the kinetic and the tadpole terms. 



1.3.2 Time evolution of the spectral function 

Within the Kadanoff-Baym calculations the full quantum information of the two-point functions 
is retained. Consequently, one has access to the spectral properties of the nonequilibrium system 
during its time evolution. The spectral function A(x, y) in our case is given by 

A(x,y) = < [<f>(x), <P(y)}- >= i [G>(x,y) - G<(x,y)]. (27) 

From the dynamical calculations the spectral function in Wigner-space for each system time 
T = (trfi2)/2 is obtained via Fourier transformation with respect to the relative time coordinate 
At = h-t 2 : 

/oo 
dAt e iAt P0 A(p,t 1 =T+At/2,t 2 =T-At/2). (28) 
-oo 

Note that a damping of the function A(p,ti,t2) in relative time At corresponds to a finite 
width r of the spectral function in Wigner-space. This width in turn can be interpreted as the 
inverse life time of the interacting scalar particle. We recall, that the spectral function - for all 
times T = t and for all momenta p - obeys the normalization 

J°° ^Po A(p,po,T) = 1 Vp,T (29) 

which is nothing but a reformulation of the equal-time commutation relation (or quantization) 
for the fields. 

In Fig. 8 we display the time evolution of the spectral function for the initial distributions 
Dl, D2 and DT for two different momentum modes | p |/m = 0.0 and | p |/m = 2.0. Since the 
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Fig. 8. Time evolution of the spectral function A(p,po,T) for the initial distributions Dl, D2, D3 and 
DT (from top to bottom) for the two momenta |p|/m = 0.0 (l.h.s.) and |p|/m = 2.0 (r.h.s.). The 
spectral function is shown for several times t ■ m = 15, 30, 60, 120, 240, 360, 480 as indicated by the 
different line types. 
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spectral functions are antisymmetric in energy for the momentum symmetric configurations 
considered, i.e. A{p, — po,T) = —A(p,pQ,T), we only show the positive energy part. For our 
initial value problem in two-times and space the Fourier transformation (28) is restricted for 
system times T to an interval At € [— 2T, 2T\. Thus in the very early phase the spectral 
function assumes a finite width already due to the limited support of the Fourier transform in 
the interval At 6 [— 2T, 2T] and a Wigner representation is not very meaningful. We, therefore, 
present the spectral functions for various system times T = t starting from t • m = 15 up to 
t • m = 480. 

For the free thermal initialization DT the evolution of the spectral function is very smooth 
and comparable to the smooth evolution of the equal-time Green function as discussed in Section 
1.2.4. In this case the spectral function is already close to the equilibrium shape at small times 
being initially only slightly broader than for late times. The maximum of the spectral function 
(for all momenta) is higher than the (bare) on-shell value and nearly keeps its position during 
the whole time evolution. This results from a positive tadpole mass shift, which is only partly 
compensated by a downward shift originating from the sunset diagram. The time evolution 
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Fig. 9. Time evolution of the on-shell energies e(p,t) of the momentum modes |p|/m = 0.0 and 
p |/to = 2.0 for the different initializations Dl, D2, D3 and DT. The on-shell self energies are extracted 
from the maxima of the time-dependent spectral functions. 



for the initial distributions Dl, D2 and D3 has a richer structure. For the distribution Dl the 
spectral function is broad for small system times (see the line for t-m — 15) and becomes a little 
sharper in the course of the time evolution (as presented for the momentum mode | p \/m = 0.0 
as well as for | p |/m = 2.0). In line with the decrease in width the height of the spectral function 
is increasing (as demanded by the normalization property (29)). This is indicated by the small 
arrow close to the peak position. Furthermore, the maximum of the spectral function (which is 
approximately the on-shell energy) is shifted slightly upwards for the zero mode and downwards 
for the mode with higher momentum. Although the real part of the (retarded) sunset self-energy 
leads (in general) to a lowering of the effective mass, the on-shell energy of the momentum modes 
is still higher than the one for the initial mass m (indicated by the 'on-shell' arrow) due to the 
positive mass shift from the tadpole contribution. For the initial distribution D3 we find the 
opposite behaviour. Here the spectral function is quite narrow for early times and increasing its 
width during the time evolution. Correspondingly, the height of the spectral function decreases 
with time. This behaviour is observed for the zero momentum mode | p |/m = 0.0 as well as for 
the finite momentum mode | p |/m = 2.0. Especially in the latter case the width for early times 
is so small that the spectral function shows oscillations originating from the finite range of the 
Fourier-transformation from relative time to energy. Although we have already increased the 
system time for the first curve to t • m = 21 (for t • m = 15 the oscillations are much stronger) 
the spectral function is not fully resolved, i.e. it is not sufficiently damped in relative time At 
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Fig. 10. Time evolution of the on-shell widths —2 Im E (p, s(p, t),t)/2 e(p, t) of the momentum modes 
| p |/m = 0.0 and | p |/m = 2.0 for the different initializations Dl, D2, D3 and DT. 



in the interval available for the Fourier transform. For later times the oscillations vanish and 
the spectral function tends to the common equilibrium shape. 

The time evolution of the spectral function for the initial distribution D2 is somehow in 
between the last two cases. Here the spectral function develops (at intermediate times) a slightly 
higher width than in the beginning before it is approaching the narrower static shape again. 
The corresponding evolution of the maximum is again indicated by the (bent) arrow. Finally, 
all spectral functions show the (same) equilibrium form represented by the solid gray line. 

As already observed in Section 1.2.4 for the equal-time Green functions, we emphazise that 
there is no universal time evolution for the nonequilibrium systems. In fact, the evolution of 
the system during the equilibration process depends on the initial conditions. On the other 
hand, the time dependence of the spectral function is only moderate such that one might also 
work with some time-averaged or even the equilibrium spectral function. In order to investigate 
this issue in more quantitative detail we concentrate on the maxima and widths of the spectral 
functions in the following. 

Since the solution of the Kadanoff-Baym equation provides the full spectral information for 
all system times the evolution of the on-shell energies can be studied as well as the spectral 
widths. In Fig. 9 we display the time dependence of the on-shell energies e(p, t) of the momentum 
modes | p |/m = 0.0 (l.h.s.) and | p |/m = 2.0 (r.h.s.) for the four initial distributions Dl, D2, 
D3 and DT. We see that the on-shell energy for the zero momentum mode increases with 
time for the initial distribution Dl and to a certain extent for the free thermal distribution 
DT (as can be also extracted form Fig. 8). The on-shell energy of distribution D3 shows a 
monotonic decrease during the evolution while it passes through a minimum for distribution 
D2 before joining the line for the initialization Dl. For momentum | p |/m = 2.0 an opposite 
behaviour is observed. Here the on-shell energy for distribution Dl (and less pronounced for 
the distribution DT) are reduced in time whereas it is increased in the case of D3. The result 
for the initialization D2 is monotonous for this mode and matches the one for Dl already for 
moderate times. Thus we find, that the time evolution of the on-shell energies does not only 
depend on the initial conditions, but might also be different for various momentum modes. It 
turns out - for the initial distributions investigated - that the above described characteristics 
change around | p |/m = 1.5 and are retained for larger momenta (not presented here). 

Furthermore, we present in Fig. 10 the time evolution of the on-shell width for the usual 
momentum modes for the different initial distributions. The on-shell width r is given by the 
imaginary part of the retarded sunset self-energy at the on-shell energy of each respective 
momentum mode as 

r = -2Im S R (p, e(p, t),t) 1 2 e(p, t). (30) 

As already discussed in connection with Fig. 8 we observe for both momentum modes a strong 
decrease of the on-shell width for the initial distribution Dl associated with a narrowing of 
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the spectral function. In contrast, the on-shell widths of distribution D3 increase with time 
such that the corresponding spectral functions broaden towards the common static shape. For 
the initialization D2 we observe an non-monotonic evolution of the on-shell widths connected 
with a broadening of the spectral function at intermediate times. Similar to the case of the 
on-shell energies we find, that the results for the on-shell widths of the distributions Dl and 
D2 coincide well above a certain system time. As expected from the lower plots of Fig. 8 the 
on-shell width for the free thermal distribution DT exhibits only a weak time dependence with 
a slight decrease in the initial phase of the time evolution. 

In summarizing this Subsection we point out, that there is no universal time evolution of 
the spectral functions for the initial distributions considered. Peak positions and widths depend 
on the initial configuration and evolve differently in time. However, we find only effects in the 
order of <10% for the on-shell energies in the initial phase of the system evolution and initial 
variations of <50% for the widths of the dominant momentum modes. Thus, depending on 
the physics problem of interest, one might discard an explicit time-dependence of the spectral 
functions and adopt the equilibrium shape. 

1.3.3 The equilibrium state 

In Section 1.2 we have seen that arbitrary initial momentum configurations of the same energy 
density approach a stationary limit for t — > oo, which is the same for all initial distributions. 
In this Subsection we will investigate whether this stationary state is the proper thermal state 
for interacting Bose particles. 

As shown before, in the present calculations within the three-loop approximation of the 2PI 
effective action we describe kinetic equilibration via the sunset self-energies and also obtain 
a finite width for the particle spectral function, since the Green functions have a non-zero 
imaginary part (in energy space). It is not obvious, however, if the stationary state obtained 
for t — > oo corresponds to the proper equilibrium state. 

In order to clarify the nature of the asymptotic stationary state of our calculations we 
first change into Wigner space. The Green function and the spectral function in energy pg are 
obtained by Fourier-transformation with respect to the relative time At = t\ — t 2 at every 
system time i = (t\ + t 2 )/2, 

G > (p,p Q ,t) = / dAt e ipoAt G > (p,t 1 = t + At/2,t 2 = t - At/2) , (31) 

J — OO 



A(p,p ,F) = / dAt e ipaAt A(p,h=t + At/2, t 2 =t- At/2). (32) 

J — oc 

We recall, that the spectral function can also be obtained directly from the Green functions 
(cf. Section 1.1) by 

A(p, Po ,t) = i [G > (p,p ,i) - G<(p,po,t)] . (33) 

Now we introduce the energy and momentum dependent distribution function N at any system 
time t by the definition 

iG < (p,p ,t) = A(p,po,t) N(p,p ,t), 

iG > (p,p ,t) = A(p,p () ,t) [JV(p,po,t) + l] (34) 

In equilibrium (at temperature T) the Green functions obey the Kubo-Martin-Schwinger re- 
lation (KMS) for all momenta p [74-76] 

G>(p,po) = e^/ T G<(p,p ) V p. (35) 
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Fig. 11. Spectral function A for various momentum modes [p[/m = 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 as a 
function of energy for late times t ■ m = 540 (lower part). Corresponding distribution function N at 
the same time for the same momentum modes (upper part). All momentum modes can be fitted with 
a single Bose function of temperature T eq /m = 1.836 and a chemical potential close to zero. 



If there exists a conserved quantum number in the theory we have, furthermore, a contribution 
of the corresponding chemical potential in the exponential function which leads to a shift of 
arguments: Po/T — > (po — fJ-)/T. In the present case, however, there is no conserved quantum 
number and thus the equilibrium state has /i = 0. 

From the KMS condition of the Green functions (35) we obtain the equilibrium form of the 
distribution function (34) at temperature T as 

N eq (p, Po ) = N eq { Po ) = epo/ l _ 1 = N bose ( Po /T), (36) 

which is the well-known Bose distribution. As is obvious from (36) the equilibrium distribution 
can only be a function of energy and not of the momentum variable in addition. 

In Figure 11 (lower part) we present the spectral function A(p,po) for the initial distribution 
D2 at late times t-m — 540 for various momentum modes |p|/m = 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 as a 
function of the energy pq. - We note, that for all other initial distributions - with equal energy 
density - the spectral function looks very similar at this time since the systems proceed to the 
same stationary state (cf. Section 1.3). - We recognize that the spectral function is quite broad, 
especially for the low momentum modes, while for the higher momentum modes its width is 
slightly lower. 

The distribution function N(po) as extracted from (34) is displayed in Fig. 11 (upper part) 
for the same momentum modes as a function of the energy P q, We find that N(po) for all 
momentum modes can be fitted by a single Bose-function with temperature T/m = 1.836. 
Thus the distribution function emerging from the Kadanoff-Baym time evolution for t — > oo 
approaches a Bose function in the energy that is independent of the momentum as demanded 
by the equilibrium form (36). 

Fig. 11 (upper part) demonstrates, furthermore, that the KMS-condition is fulfilled not only 
for on-shell energies, but for all po- We, therefore, have obtained the full off-shell equilibrium 
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state by integrating the Kadanoff-Baym equations in time. In addition, the limiting stationary 
state is the correct equilibrium state for all energies po, i.e. also away from the quasi-particle 
energies. 

We note in closing this Subsection, that the chemical potential fi - used as a second fit 
parameter - is already close to zero for these late times as expected for the correct equilibrium 
state of the neutral (/> 4 -theory that is characterized by a vanishing chemical potential fi in 
equilibrium which is actually achieved in the calculations for t —* oo. This, at first sight, seems 
trivial but is a consequence of our 'exact' treatment. In contrast, the Boltzmann equation (cf. 
Section 1.4) in general leads to a stationary state for t — > oo with a finite chemical potential. 
We will attribute this failure of the Boltzmann approach to the absence of particle number 
nonconserving processes in the quasi-particle limit (see below) . 



1.4 Full versus approximate dynamics 

The Kadanoff-Baym equations studied in the previous Sections represent the full quantum- 
field theoretical equations on the single-particle level. However, its numerical solution is quite 
involved and it is of strong interest to investigate, in how far approximate schemes deviate 
from the full calculation. Nowadays, transport models are widely used in the description of 
quantum systems out of equilibrium (cf. Introduction). Most of these models work in the 'quasi- 
particle' picture, where all particles obey a fixed energy-momentum relation and the energy is no 
independent degree of freedom anymore; it is determined by the momentum and the (effective) 
mass of the particle. Accordingly, these particles are treated with their <5-function spectral 
shape as infinitely long living, i.e. stable objects. This assumption is very questionable e.g. for 
high-energy heavy ion reactions, where the particles achieve a large width due to the frequent 
collisions with other particles in the high-density and/or high-energy regime. Furthermore, this 
is doubtful for particles that are unstable even in the vacuum. The question, in how far the 
quasiparticle approximation influences the dynamics in comparison to the full Kadanoff-Baym 
calculation, is of widespread interest [53]. 



1.4.1 Derivation of the Boltzmann approximation 

In the following we will give a short derivation of the Boltzmann equation starting directly 
from the Kadanoff-Baym dynamics in the two-time and momentum-space representation. This 
derivation is briefly reviewed since we want i) to emphasize the link of the full Kadanoff- 
Baym equation with its approximated version and ii) to clarify the assumptions that enter the 
Boltzmann equation. 

Since the Boltzmann equation describes the time evolution of distribution functions for 
quasi-particles we first consider the quasi-particle Green-functions in two-time representation 
for homogeneous systems 

gL(p.^') ^ ^(TP) e ± -p(*-*') + [N qp (±p) + 1] e=F*-p(*-f) } (3?) 

GJ, 9P (P,M') = \{ TAWTP) e±-p(*-*') ± [N qp (±p) + 1] e^*-*') } 
G^ qp {p,t,t') = \{ ±AWTP) e^'-'') T [N qp (±p) + 1] e=F*-p(*-f) } 

G>^ qp {p,t,t') = ^{ N qp ( T p) e ±*"p0-O + [AV(±P) + 1] e=F*"p(*-f ) } , 

where for each momentum p the Green functions are freely oscillating in relative time t — t' 
with the on-shell energy u p . The time-dependent quasi-particle distribution functions are given 
with the energy variable fixed to the on-shell energy as N qp (p,i) = N(p,p n — u) p ,t), where 
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the on-shcll energies u> p might depend on time as well. Such a time variation e.g. might be due 
to an effective mass as generated by the time-dependent tadpole self-energy. In this case the 
on-shell energy reads 

Wp (t) = yW^ + ^U*)- (38) 

Vice versa we can define the quasi-particle distribution function by means of the quasi-particle 
Green functions at equal times t as 



N qp (p,t) 



[P,t,i) 



1 



2ajp{F) iG< Vtgp (p,i,t) 



(39) 



G< 01W (P,*,*) ~ G<^ qp {p,i,t) 



Using the equations of motions for the Green functions in diagonal time direction (19) (ex- 
ploiting G^ n (p,i,t) = — [G^(p, i, £)]*) the time evolution of the distribution function is given 

by 



d- t N qp (p, t) = -Re{ I< qp (p, i, t) } m Im { I< 2 . qp (p, i, t) } . 



(40) 



The time derivatives of the on-shcll energies cancel out since the quasiparticle Green functions 
obey 

G<Jp, t, t) = co 2 p (t) G< (p, t, t) (41) 

as seen from (37). Furthermore, it is remarkable that contributions containing the energy up - 
as present in the equation of motion for the Green functions (19) - no longer show up. The time 
evolution of the distribution function is entirely determined by (equal-time) collision integrals 
containing (time derivatives of the) Green functions and self-energies. 

I< qp (p, i, t) = jf dt' (e<(j>, t, t') G> , gp (p, t', t) - E>(p, i, t') G<^ qp (p, t', t)) , (42) 
Ii,2-, qP (P, i, t)= [ dt' (Z<(p, i, t') G>^ qp {p, t', t) - 2> (p, t, t') G<^ qp (p, t', t) 

J t y 

Since we are dealing with a system of on-shell quasi-particles within the Boltzmann approxima- 
tion, the Green functions in the collision integrals (42) are given by the respective quasi-particle 
quantities of (37). Moreover, the collisional self-energies (17) are obtained in accordance with 
the quasi-particle approximation as 

4<P.U')=-4 /^/(g,/^ (2^<%.-<,-r-s> ^-±^ (43, 



{ iV(Tq) N qp ( T r) N qp ( T s) e -H[*-f ][±a»,±« r ±<-.] 
+ 3 N qp (Tq) N qp ( T r) [JV w (±s) + l ] e+* I*"*' 1 [^'^ 1 
+ 3 N qp (T<l) [N qp (±r) + l}[N qp (±s) + l}e+^ t - t '^ ± ^^-^ 

+[N qp (±q) + l}[N qp (±r) + l}[N qp (±s) + l}e+^ t - t '^^^^ }. 
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For a free theory the distribution functions N qp (p) are obviously constant in time which, of 
course, is no longer valid for an interacting system out of equilibrium. Thus one has to specify the 
above expressions for the quasi-particle Green functions (37) to account for the time dependence 
of the distribution functions. 

The actual Boltzmann approximation is defined in the limit, that the distribution functions 
have to be taken always at the latest time argument of the two-time Green function [53] . Ac- 
cordingly, for the general non-equilibrium case, we introduce the ansatz for the Green functions 
in the collision term 

G L (p ' t ' t ' )= 2^ { ^ (TP ''- )e±MpM) + [^p(±P^^) + 1] e^(*-*')} (44) 

G^ qp {p,t,t') = X - { TN qp ( T p,t max ) e^*-*') ± [N qp (±p,t max ) + 1] e =F**p(*-*'> }. 

with the maximum time t max = max(t, t'). The same ansatz is employed for the time dependent 
on-shell energies which enter the representation of the quasi-particle two-time Green functions 
(44) with their value at t max , i.e. u> p — u) p (t max = max{t, t')). 

The collision term contains a time integration which extends from an initial time to to the 
current time i. All two-time Green functions and self-energies depend on the current time i as 
well as on the integration time t' < t. Thus only distribution functions at the current time, 
i.e. the maximum time of all appearing two-time functions, enter the collision integrals and the 
evolution equation for the distribution function becomes local in time. Since the distribution 
functions are given at fixed time t, they can be taken out of the time integral. When inserting the 
expressions for the self-energies and the Green functions in the collision integrals the evolution 
equation for the quasi-particle distribution function reads 

A 2 f d d q f d d r fd d s . . . . 1 



WP,*) - yJj^J^J^ (2.)^>(p-<,-r-s) («) 

| K,N-^-r,N-sr N P,^-^-r,N-s}jdt' C0S( [ t- 1' } [ LO p + C0 q + UJ r + UJ S } ) 
+ ^p.N-c l .N~r.iN s ^N p .^ cl .^ r ^N s .i fdt' COs([i-t'} [uJp+UJ^ + GJr-UJs}) 

+m p .^-c l jN r ^N s ^-N p ^^N r , i N s ,i l fdt' cos{[t-t'] [a;p+a; q -a; r -a; s ]) 

J t 

+ ^p^ q ,^ r ,^ s , t --^ p # q , t -A^ r , t -^ s , t -]^' COs([i-i'] [cJp-LU^-UJr-LUs}) J , 

where we have introduced the abbreviation N p j — N qp (p,t) for the distribution function at 
current time i and N p j — N qp (p, t) + 1 for the according Bose factor. Furthermore, a possible 
time dependence of the on-shell energies is suppressed in the above notation. 

The contributions in the collision term (45) for particles of momentum p are ordered as they 
describe different types of scattering processes where, however, we always find the typical gain 
and loss structure. The first line in (45) corresponds to the production and annihilation of four 
on-shell particles (0 — > 4, 4 — > 0), where a particle of momentum p is produced or destroyed 
simultaneous with three other particles with momenta q, r, s. The second line and the forth line 
describe (1 — > 3) and (3 — > 1) processes where the quasi-particle with momentum p is the single 
one or appears with two other particles. The relevant contribution in the Boltzmann limit is 
the third line which respresents (2 — > 2) scattering processes; quasi-particles with momentum p 
can be scattered out of their momentum cell by collisions with particles of momenta q (second 
term) or can be produced within a reaction of on-shell particles with momenta r, s (first term). 
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The time evolution of the quasi-particle distribution is given as an initial value problem 
for the function N qp (p) prepared at initial time to- For large system times t (compared to the 
initial time) the time integration over the trigonometric function results in an energy conserving 
(5-function: 

r* i - 

_ lim / dt' cos((t - t') us) = _ lim - sin((i - t ) us) = n 6(us) . (46) 

t — to— »oo J ta t— to— >oo US 



Here us = us p ± u> q ± us r ± us s represents the energy sum which is conserved in the limit i—to — > oo 
where the initial time to is considered as fixed. In this limit the time evolution of the distribution 
function amounts to 

A 2 f d d q fd d r fd d s . . . +1 1 



d t N qp ^t) = -j^j—j— (2.)^ 2ujp2 ^ 2u)r2uJs (47) 



|^ P;t - JV qif N rf JV Sjf - N Vi { N^ f N Ti { JV s>f ] 6 ( %> + q+r+s) (5(w p +w q +cj r +cj s ) 
+3[JVp,f JVq,f JVr,f JV.,f ~ Wp,f JVq,f JVp,f JVs,f] <5 (d) (p+q+r-s) <5(w p + Wq +w r -^ s ) 



+3[JV Pl f JV q ,f JVr,f JV.,f - N Pti N Tti JV 8>f ] ^ d \p + q-r-s) <S(cj p +w q -w r -w s ) 

+ [JVp,f JVq,f ^p,f JV.,f - JVp,t #q,t #r,t #s,t] 5 ( %> - q- r - s) <f(w p - W q - W r - W s ) | . 

In the energy conserving long-time limit (46) only the 2 — > 2 scattering processes are non- 
vanishing, because all other terms do not contribute since the energy <5-functions can not be 
fulfilled for on-shell quasi-particles. Furthermore, the system evolution is explicitly local in 
time because it depends only on the current configuration; there are no memory effects from 
the integration over past times as present in the full Kadanoff-Baym equation. 

In the following we will solve the energy conserving Boltzmann equation for on-shell parti- 
cles: 

o ^ , s X2 f dd< l f ddr f dd - s , n Nrf+i 1 



[NpjN^t N TtI N s j - N Pi{ JV q , t - N r ,i N Bti ] tf^p+q-r-s) 5{u> p +us^- u> r -ui s ) . 

The numerical algorithm employed for the solution of (48) is basically the same as for the 
solution of the Kadanoff-Baym equation (cf. Section 1.2). We explicitly calculate the time 
integral in (45). Energy conservation can be assured by a precalculation including a shift of the 
lower boundary to to earlier times. We note, that in contrast to the Kadanoff-Baym equation 
no correlation energy is generated in the Boltzmann limit! 

In addition to the procedure presented above we calculate the actual momentum-dependent 
on-shell energy for every momentum mode by a solution of the dispersion relation including 
contributions from the tadpole and the real part of the (retarded) sunset self energy. In this 
way one can guarantee that at every time t the particles are treated as quasi-particles with the 
correct energy-momentum relation. 

Before presenting the actual numerical results we comment on the derivation of the Boltz- 
mann equation within the conventional scheme that is different from the one presented above. 
Here, at first the Kadanoff-Baym equation (in coordinate space) is transformed to the Wigner 
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Kadanoff-Baym distribution d2 /Jm = J \& 

Fig. 12. Evolution of the Green function in momentum space. The equal time Green function is 
displayed for various times t ■ m — 0, 15, 30, 45, 75, 150. Starting from an initially non-isotropic shape 
it develops towards a rotational symmetric distribution in momentum space. 

representation by Fourier transformation with respect to the relative coordinates in space and 
time. The problem then is formulated in terms of energy and momentum variables together 
with a single system time. For non-homogeneous systems a mean spatial coordinate is nec- 
essary as well. As a next step the 'semiclassical approximation' is introduced, which consists 
of a gradient expansion of the convolution integrals in coordinate space within the Wigner 
transformation. For the time evolution only contributions up to first order in the gradients 
are kept. Finally, the quasi-particle assumption is introduced as follows: The Green functions 
appearing in the transport equation - explicitly or implicitly via the self-energies - are written 
in Wigner representation as a product of a distribution function N and the spectral function 
A. The quasi-particle assumption is then realized by employing a 5-like form for the spectral 
function which connects the energy variable to the momentum. By integrating the first order 
transport equation over all (positive) energies, furthermore, the Boltzmann equation for the 
time evolution of the on-shell distribution function (48) is obtained. 

Inspite of the fact, that the Bolzmann equation (48) can be obtained in different subsequent 
approximation schemes, it is of basic interest, how its actual solutions compare to those from 
the full Kadanoff-Baym dynamics. 



1.4.2 Boltzmann vs. Kadanoff-Baym dynamics 

In the following we will compare the solutions of the Boltzmann equation with the solution of 
the Kadanoff-Baym theory. We start with a presentation of the non-equilibrium time evolution 
of two colliding particle accumulations (tsunamis) within the full Kadanoff-Baym calculation 
(see Figure 12). 

During the time evolution the bumps at finite momenta (in p x direction) slowly disappear, 
while the one close to zero momentum - which initially stems from the vacuum contribution to 
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Fig. 13. Decrease of the quadrupole moment in time for different coupling constants A/m = 8 (2) 16 
for the full Kadanoff-Baym calculation and the Boltzmann approximation. 



the Green function - is decreased as seen for different snapshots at times t ■ m = 0, 15, 30, 45, 
75, 150 in Fig. 12. The system with initially apparent collision axis slowly merges - as expected 
- into an isotropic final distribution in momentum space. 

For the comparison between the full Kadanoff-Baym dynamics and the Boltzmann approx- 
imation we concentrate on equilibration times. To this aim we define a 'quadrupole' moment 
for a given momentum distribution n(p) at time T as 



(2ir) d VP ' ' 

which vanishes for the equilibrium state. For the Kadanoff-Baym case we employ the actual 
distribution function by the relation 

1 



n(p, t) = ^/G<^(p, t, t) G< T (p, t, f) - - . (50) 

Note that when constructing the distribution function by means of equal-time Green func- 
tions the energy variable has been effectively integrated out. This has the advantage that the 
distribution function is given independently of the actual on-shell energies. 

The relaxation of the quadrupole moment (50) has been studied for two different initial 
distributions: The evolution of distribution d2 is displayed in Fig. 12 while for distribution 
dl the position and the width of the two particle bumps have been modified. The calculated 
quadrupole moment (50) shows a nearly exponential decrease with time (cf. Fig. 13) and one 
can extract a relaxation rate 7q via the relation 

Q(t)~exp(-r g t). (51) 



Fig. 14 shows for both initializations that the relaxation in the full quantum calculation 
occurs faster for large coupling constants than in the quasi-classical approximation, whereas 
for small couplings the equilibration times of the full and the approximate evolutions are com- 
parable. We find that the scaled relaxation rate ^q/A 2 is nearly constant in the Boltzmann 
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Fig. 14. Relaxation rate (divided by the coupling A squared) for Kadanoff-Baym and Boltzmann 
calculations as a function of the interaction strength. For the two different initial configurations the 
full Kadanoff-Baym evolution leads to a faster equilibration. 



case, but increases with the coupling strength in the Kadanoff-Baym calculation (especially for 
initial distribution d2). 

These findings are readily explained: Since the free Green function - as used in the Boltz- 
mann calculation - has only support on the mass shell, only (2 <-> 2) scattering processes are 
described in the Boltzmann limit. All other processes with a different number of incoming and 
outgoing particles vanish (as noted before). Within the full Kadanoff-Baym calculation this is 
different, since here the spectral function - determined from the self-consistent Green function 
- aquires a finite width. Thus the Green function has support at all energies although it drops 
fast far off the mass shell. Especially for large coupling constants, where the spectral function 
is sufficiently broad, the three particle production process gives a significant contribution to 
the collision integral. Since the width of the spectral function increases with the interaction 
strength, such processes become more important in the high coupling regime. As a consequence 
the difference between both approaches is larger for stronger interactions as observed in Fig. 14. 
For small couplings A/m in both approaches basically the usual 2^2 scattering contributes 
and the results for the thermalization rate -Tq are quite similar. 

In summarizing this Subsection we point out that the full solution of the Kadanoff-Baym 
equations does include «-> 4, 1 <-> 3 and 2 «-> 2 off-shell collision processes which - in compari- 
son to the Bolzmann on-shell 2 2 collision limit - become important when the spectral width 
of the particles reaches ~ 1/3 of the particle mass. On the other hand, the simple Boltzmann 
limit works surprisingly well for smaller couplings and those cases, where the spectral function 
is sufficiently narrow. 
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2 Off-shell relativistic transport theory 

Formal derivations of off-shell transport equations have been presented about 40 years ago 
by Kadanoff and Baym [13] but actual solutions have barely been addressed. This Section 
is devoted to a transparent derivation of generalized transport equations including numerical 
illustrations as well as detailed comparisons to the full theory. 



2.1 Derivation of off-shell transport theory 

The derivation of generalized transport equations starts by rewriting the Kadanoff-Baym equa- 
tion for the Wightman functions in coordinate space (xi = (ti, xi), x 2 = {t 2 ,x 2 )) (8) as 

[d^+m 2 + S 5 ( Xl )} iG^(x u x 2 ) = ilf(x u x 2 ). (52) 

The collision terms on the r.h.s. of (52) are given in D = d + 1 space-time dimensions by 
convolution integrals over coordinate space self-energies and Green functions: 

lf(x u x 2 ) = -f 1 d D z [S>(x 1 ,z)-S<{x 1 ,z)] G^(z,x 2 ) (53) 

ft 



f 2 d D z E$(x u z) [G > (z,x 2 )-G < (z,x 2 )] . 



't 

In the general case of an arbitrary (scalar) quantum field theory S 5 is the local (non-dissipative) 

part of the path self-energy while resemble the non-local collisional self-energy contribu- 
tions. In the representation (53) the integration boundaries are exclusively given for the time 
coordinates, while the integration over the spatial coordinates extends over the whole spatial 
volume from — oo to +oo in d dimensions. 

Since transport theories are formulated in phase-space one changes to the Wigner repre- 
sentation via Fourier transformation with respect to the rapidly varying ('intrinsic') relative 
coordinate Ax = x\ — x 2 and treats the system evolution in terms of the ('macroscopic') mean 
space-time coordinate x = ( Xl + x 2 )/2 and the four-momentum p = (p ,p)- The functions in 
Wigner space are obtained as 

/oo 
d D Ax e +iAx » p " F( Xl = x + Ax/2, x 2 =x- Ax/2). (54) 
-oo 

For the formulation of transport theory in the Wigner representation we have to focus not only 
on the transformation properties of ordinary two-point functions as given in (54), but also of 
convolution integrals as appearing in Eq. (53). A convolution integral in D dimensions (for 
arbitrary functions F, G) , 

/oo 
d D z F{x u z) G{z,x 2 ) (55) 
-OO 

transforms as 

/oo 
d D Ax e +iAx » p " H( Xl ,x 2 ) (56) 
-OO 



/OO pOC 

d D Ax e +iAx » p " / d D z F( Xl ,z) G(z,x 2 ) 
-oo J —oo 

e 2 p M x * ! [F(p, x) G[p\ x )\ 



x'—x, p'—p 
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In accordance with the standard assumption of transport theory we assume that all func- 
tions only smoothly evolve in the mean space-time coordinates and thus restrict to first order 
derivatives. All terms proportional to second or higher order derivatives in the mean space- 
time coordinates (also mixed ones) will be dropped. Thus the Wigner transformed convolution 
integrals (55) are given in first order gradient approximation by, 

H(p,x) = F(p,x) G(p,x) + i±{F(p,x),G( P ,x)} + 0(d 2 x ) , (57) 

using the relativistic generalization of the Poisson bracket 

{ F{p, x) , G{ P , x) } = dl F(p, x) ■ d% G{p, x) - d» F{p, x) ■ G{p, x) . (58) 

In order to obtain the dynamics for the spectral functions within the approximate scheme 
we start with the Dyson-Schwinger equations for the retarded and advanced Green functions 
in coordinate space (14,15). - Note that the convolution integrals in (14) and (15) extend 
over the whole space and time range in contrast to the equations of motion for the Wightman 
functions given in (52) and (53)! - The further procedure consists in the following steps: First we 

i) transform the above equations into the Wigner representation and apply the first order 
gradient approximation. In this limit the convolution integrals yield the product terms and 
the general Poisson bracket of the self-energies and the Green functions { E R / A , G R I A }. We, 
furtheron, represent both equations in terms of real quantities by the decomposition of the 
retarded and advanced Green functions and self-energies as 

G R l A =Re G R ± ilm G R —Re G R =F iA/2 , A=T 2 Im G R ' A , 

(59) 

£ R / A =Re E R ± ilm Z R =Re S R =p iT/2 , r=T 2 Im E R / A . 

We find that in Wigner space the real parts of the retarded and advanced Green functions and 
self-energies are equal, while the imaginary parts have opposite sign and are proportional to 
the spectral function A and the width P, respectively. The next step consists in 

ii) the separation of the real part and the imaginary part of the two equations for the re- 
tarded and advanced Green functions, that have to be fulfilled independently. Thus we obtain 
four real- valued equations for the self-consistent retarded and advanced Green functions. In the 
last step 

iii) we get simple relations by linear combination of these equations, i.e. by adding/subtracting 
the relevant equations. 

This finally leads to two algebraic relations for the spectral function A and the real part of 
the retarded Green function Re G R in terms of the width r and the real part of the retarded 
self-energy ReS R as [64,65]: 

[pl-p 2 -m 2 - E s + ReS R ] Re G R = 1 + \ rA, (60) 



[Po ~ P 111 — U + Re S ] A = r Re G . (61) 

Note that all terms with first order gradients have disappeared in (60) and (61). A first 
consequence of (61) is a direct relation between the real and the imaginary parts of the re- 
tarded/advanced Green function, which reads (for t ^ 0): 

ReG R = ^- p2 - m2 -^- Re ^ R A. (62) 
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Inserting (62) in (60) we end up with the following result for the spectral function and the real 
part of the retarded Green function 

r r 

(63) 



[p§ - p 2 - to 2 - S s - ReS R ] 2 + 1 2 /4 M 2 + T 2 /4 



ReC R = lP 2 o-P 2 -™ 2 -Z 5 -ReZ R ] = M 
" ' '™ 2 -p 2 -m 2 - Z s - ReZ R } 2 +P 2 /4 M 2 + T 2 /4' 



where we have introduced the mass-function M(p,x) in Wigner space: 

M{p,x) =pl-p 2 -m 2 - S\x) - ReS R {p,x) . (65) 

The spectral function (63) shows a typical Breit- Wigner shape with energy- and momentum- 
dependent self-energy terms. Although the above equations are purely algebraic solutions and 
contain no derivative terms, they are valid up to the first order in the gradients! 

In addition, subtraction of the real parts and adding up the imaginary parts lead to the 
time evolution equations 

ffd?A=±{E s + ReZ: R ,A} + \{r , ReG R } , (66) 

p» dl Re G R = \ { S s + Re S R , Re G R } - i { f , A } . (67) 

The Poisson bracket containing the mass-function M leads to the well-known drift operator 
pfj, Qxp ^£ or an ar bitrary function F), i.e. 

{M,F} = {p 2 -p 2 -m 2 -E 5 -ReS R ,F} (68) 

= 2p»d*F - {E 5 + Re£ R , F} , (69) 

such that the first order equations (66) and (67) can be written in a more comprehensive form 
as 

{M,A} = {P,ReG R }, (70) 



{M,ReG R } = -\{r , A}. (71) 

When inserting (63) and (64) we find that these first order time evolution equations are solved 
by the algebraic expressions. In this case the following relations hold: 

{M,A} = { r,ReG R } = {M,r } ^-gjt? > ( ?2 ) 



1 m r ii 

{ M , i?e G R } = --{r,A} = {M,F} ^ 2+ ^ /4]2 - (73) 

Thus we have derived the proper structure of the spectral function (63) within the first-order 
gradient (or semiclassical) approximation. Together with the explicit form for the real part of 
the retarded Green function (64) we now have fixed the dynamics of the spectral properties, 
which is consistent up to first order in the gradients. 
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2.1.1 Kadanoff-Baym transport 

As a next step we rewrite the memory terms in the collision integrals (53) such that the time 
integrations extend from — oo to +oo. In this respect we consider the initial time to = — oo 
whereas the upper time boundaries t\, t 2 are taken into account by ©-functions, i.e. 

< < 
I>(x u x 2 ) = - d D x' 0{h-t') [E>( Xl , x')-E < {x 1 ,x')] G>(x',x 2 ) 

J — OO 

/CO < 
d D x' S>(x u x') 0(t 2 -t') [G > (x',x 2 )-G < (x',x 2 )] 
-OO 



/OO 
d 
-oo 



d D x' E R ( Xl ,x') G^{x',x 2 ) + S^( Xl ,x') G A {x',x 2 ) . (74) 

o 

We now perform the analogous steps as invoked before for the retarded and advanced Dyson- 
Schwingcr equations. We start with a first order gradient expansion of the Wigner transformed 
Kadanoff-Baym equation using (74) for the memory integrals. Again we separate the real and 
the imaginary parts in the resulting equation, which have to be satisfied independently At the 
end of this procedure we obtain a generalized transport equation [13,20,24,64,65,67,69,77]: 

2 p" d* iG^ - { S s + Re E R , iG^ } - { iS^ , ReG R } = iS< iG> - iS> iG< 



{M,iG^} - {iS^ , ReG R } = iS< iG> - iS> iG< (75) 



as well as a generalized mass-shell equation 

[p 2 -m 2 -E 5 -ReS R ] iG^ = iS^ ReG R + ^{iS>,iG<} - ^{iS< , iG> } (76) 



M 

with the mass- function M specified in (65). Since the Green function G > (x ll x 2 ) consists of an 
antisymmetric real part and a symmetric imaginary part with respect to the relative coordinate 
xi — x 2 , the Wigner transform of this function is purely imaginary. It is thus convenient to 

represent the Wightman functions in Wigner space by the real- valued quantities iG^ (p, x) . 
Since the collisional self-energies obey the same symmetry relations in coordinate space and in 

phase-space, they will be kept also as iS > (p,x) furthcron. 

In the transport equation (75) one recognizes on the l.h.s. the drift term p M c£f iG > , as well 
as the Vlasov term with the local self-energy S 5 and the real part of the retarded self-energy 
ReE R . On the other hand the r.h.s. represents the collision term with its typical 'gain and 
loss' structure. The loss term iE > «G < (proportional to the Green function itself) describes the 
scattering out of a respective phase-space cell whereas the gain term i£ < iG y takes into account 

scatterings into the actual cell. The last term on the l.h.s. { iS > , Re G R } is very peculiar since 
it does not contain directly the distribution function iG K . This second Poisson bracket vanishes 
in the quasiparticle approximation and thus does not appear in the on-shcll Boltzmann limit. As 

demonstrated in detail in Refs. [13,64,65,67,69,77] the second Poisson bracket { iZ 1 , Re G R } 
governs the evolution of the off-shell dynamics for noncquilibrium systems. 
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Although the generalized transport equation (75) and the generalized mass-shell equation 
(76) have been derived from the same Kadanoff-Baym equation in a first order gradient expan- 
sion, both equations are not exactly equivalent [24,64,77]. Instead, they deviate from each other 

by contributions of second gradient order, which are hidden in the term { iE > , Re G R } (see 
below or Refs. [24,77] for extended discussions). This raises the question: which one of these 
two equations has to be considered of higher priority? The question is answered in practical 
applications by the prescription of solving the generalized transport equation (75) for iG K in 
order to study the dynamics of the nonequilibrium system in phase-space. Since the dynamical 
evolution of the spectral properties is taken into account by the equations derived in first order 
gradient expansion from the retarded and advanced Dyson- Schwinger equations, one can ne- 
glect the generalized mass-shell equation (76). Thus for our actual numerical studies in Section 
2.2 we will use the generalized transport equation (75) supported by the algebraic relations 
(63) and (64). 



2.1.2 Transport in the Botermans-Malfliet scheme 

Furthermore, one recognizes by subtraction of the iG > and iG < mass-shell and transport equa- 
tions, that the dynamics of the spectral function A = iG > — iG< is determined in the same 
way as derived from the retarded and advanced Dyson-Schwingcr equations (63) and (70). 
The inconsistency between the two equations (75) and (76) vanishes since the differences are 
contained in the collisional contributions on the r.h.s. of (75). 

In order to evaluate the { i£ < 1 Re G R }-term on the l.h.s. of (75) and to explore the dif- 
ferences between the KB- and BM-form of the transport equations (see below) it is useful to 
introduce distribution functions for the Green functions and self-energies as 



iG<{p,x) = N{p,x) A(p,x) , iG>(p,x) = [1 + N{p,x)} A(p,x) , 



(77) 



iS<(p,x) = N s ( P ,x) P(p,x) , iS>(p,x) = [1 + N s (p,x)}P(p,x) . 



(78) 



In equilibrium the distribution function with respect to the Green functions N and the self- 
energies N s are given as Bose functions in the energy po at given temperature; they thus are 
equal in equilibrium but in general might differ out-of-equilibrium. Following the argumentation 
of Botermans and Malfliet [24] the distribution functions N and N E in (77) should be identical 
within the second term of the l.h.s. of (75) in order to obtain a consistent first order gradient 
expansion (without hidden higher order gradient terms) . In order to demonstrate their argument 
we write 



r n 2 



rN + k. 



(79) 



The 'correction' term 



K 



r(N s - N) = {iS < iG > - iS > iG < )A- 



(80) 



is proportional to the collision term of the generalized transport equation (75), which itself 
is already of first order in the gradients. Thus, whenever a distribution function N s appears 
within a Poisson bracket, the difference term (N s — N) becomes of second order in the gradients 
and should be omitted for consistency. As a consequence N s can be replaced by N and thus 
the self-energy S K by G < ■ t/A in the Poisson bracket term {S < ,ReG R }. The generalized 
transport equation (75) then can be written in short-hand notation 



1 



A r 



{M, iG<} - l{r,M-zG<} 



= iS 



< „• 



iG> - iS 



> „• 



G< 



(81) 
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with the mass-function M (65). The transport equation (81) within the Botermans-Malfiict 
(BM) form resolves the discrepancy between the generalized mass-shell equation (76) and the 
generalized transport equation in its original Kadanoff-Baym (KB) form (75). 

However, it is presently not clear in how far the generalized transport equations in KB-form 
(75) or in BM-form (81) reproduce the same dynamics as for the full Kadanoff-Baym theory 
(8). Moreover, the differences in the time evolution of nonequilibrium systems between the KB- 
form (75) and BM-form (81) are not known either. We will thus perform an explicit comparison 
between the different limits for 4> -theory in 2+1 dimensions in the following Section. 

2.2 Numerical studies on off-shell transport 

In this Section we will perform numerical studies of the dynamics inherent in the generalized 
transport equations in comparison to solutions of the full Kadanoff-Baym equations (8) for the 
</> 4 -theory in 2+1 space-time dimensions within the three- loop approximation for the effective 

action. This fixes the self-energies S° and il, in (75) and (76) to be the same as in the case 
of the full Kadanoff-Baym theory given by the tadpole (9) and the sunset (10) contributions, 
respectively. 

For the first investigation we concentrate on the dynamics of the generalized transport equa- 
tion (75). As in Section 1.2 we restrict ourselves to homogeneous systems in space. Consequently 
the derivatives with respect to the mean spatial coordinate x vanish, such that the generalized 
transport equation (75) reduces to 

2p d t iG< - { E 6 +Re E R , iG< } - { iS< , ReG R } = iS< iG> - iE> iG< (82) 

with the simplified Poisson brackets (for arbitrary functions F, G) 

{F(p,p ,t),G(p,p ,t)} (83) 

= d P0 F(p, Pa ,t) d t G(p,p ,t) - d t F(p,p , t) d Po G(p,p , t) . 

In order to obtain a numerical solution the Kadanoff-Baym equation (8) as well as the gener- 
alized transport equations (75) are transformed to momentum space. For the actual numerical 
integration of the generalized transport equation we will use the same grid in momentum space 
as employed for the full KB-theory in Section 1.2. 

2.2.1 Phases of Equilibration 

In order to specify the problem for the first order transport equation in time, the initial state has 
to be fixed, i.e., the initial Green function iG < (p,po, t = 0) has to be specified for all momenta 
p and all energies po- According to (77) iG < (p,po, t — 0)_can be written as a product of a distri- 
bution function N(p,po, t = 0) and a spectral function A(p,p , t = 0). In order to determine N 
and A we employ an iterative scheme: We first assume initial momentum-distribution functions 
n(p, t = 0) equivalent to those used in the investigation of the full Kadanoff-Baym theory in 
Section 1.2 in order to allow for a comparison of the two schemes. These initial momentum 
distributions n(p, t = 0) then fix iG < (p, t = 0, t = 0) on the time-diagonal axis by [72] 

2w p iG < (p,t = 0,t = 0) = n(p,f = 0) +n(-p,i = 0) + 1. (84) 

Next, the complete initial phase-space distribution function N(p,po, t = 0) in Wigner represen- 
tation as well as the spectral function A(p lPOl t — 0) is_specified as a function of the energy 
Po in an iterative way. As a starting point we assume iV(p,p_ ) to be constant in the energy 
variable po. In order to determine the initial Green function iG < we employ the self-consistent 
iteration procedure for the full spectral function as_ described in Appendix D of Ref. [72] using 
the specified nonequilibrium distribution function N. The iteration process then yields the fully 
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self-consistent spectral function A(p,po,t- 
thc initial Green function via 



: 0) for this initial distribution and thus determines 



iG<(p,po,t=0) - N(p,p ,t=0) A(p,po,t=0). 



(85) 



The advantage of the initialization prescription introduced above is that the actual spectral 
function - directly obtained by A = iG > — iG < from the Green functions - complies with the 
one determined from the self-energies (63) in accordance with the first order gradient expansion 
scheme. During the nonequilibrium time evolution this correspondence is maintained since the 
analytic expression for the spectral function already is a solution of the generalized transport 
equation itself. Furthermore, the real part of the retarded Green function, that enters the 
peculiar second Poisson bracket on the l.h.s. in (75), can be taken in the first order scheme (64) 
which simplifies the calculations considerably. 

Now we turn to the actual solutions of the generalized transport equation in the KB 
form (75). In Fig. 15 (left part) we show the time evolution of the equal-time Green func- 
tion iG < (\ p |, t, t) for the polar symmetric initial states Dl, D2 and D3. Displayed are several 
momentum modes | p |/m = 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 of the equal-time Green function on a 
logarithmic time scale. As in the full Kadanoff-Baym theory (right part) one finds that for all 
initializations the quantum system approaches a stationary state for t — > oo, i.e. all momentum 
modes approach a constant. 
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Fig. 15. Left part: Evolution of several momentum modes j p\/m — 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 of the 
equal-time Green function on a logarithmic time scale for the different initializations Dl, D2 and D3 
for the generalized transport equation (75). Right part: same as above but for the full Kadanoff-Baym 
equation (8). 



However, the respective momentum modes of the different initializations do not achieve 
identical values for t — > oo, as seen in particular for the low momenta | p |/m = 0.0, 0.8 in Fig. 
15 (left part). This is not surprising since the various initializations - obtained within the self- 
consistent scheme described above - do not correspond to exactly the same energy. This is why 
the respective long-time limits differ slightly. The small difference in energy is, of course, most 
prominently seen in the low momentum (energy) modes. Moreover, the dynamics within the 
generalized transport equation (75) is in general very similar to the full Kadanoff-Baym theory 
(right part). For all three initial states we find (apart from the very initial phase t ■ m < 5) the 
same structures during the equilibration process. In particular for the initializations Dl and 
D2 the characteristic overshooting for the low momentum modes is seen as in the full quantum 
evolution, which does not show up in solutions of the corresponding on-shell Boltzmann limit. 
Since in the Boltzmann limit a strictly monotonous evolution of the momentum modes is seen 
(cf. Ref. [72]) this overshooting has to be attributed to an off-shell quantum effect. Even the 
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positions of the maxima are in a comparable range: For the initialization Dl they are shifted to 
slightly larger times and are a little bit lower than in the full calculation; the same holds for the 
initial state D2. The initial distribution D3 yields a monotonous behaviour for all momentum 
modes within the generalized transport formulation which is again in a good agreement with 
the full dynamics. 

Some comments are worthwhile with respect to the comparison performed above: The spec- 
tral function in the Kadanoff-Baym calculation is completely undetermined in the initial state; 
it develops during the very early phase to an approximate form (which in the following still 
evolves in time). In contrast to this, the spectral function in the generalized transport formu- 
lation (75) has a well-defined structure already from the beginning. This principle difference 
results from the fact, that in the Kadanoff-Baym case we deal with a true initial value problem 
in the two time directions (ti, £2)- Additionally, the relative time integral in (t\ — £2) ~~ to obtain 
the spectral function in energy po by Wigner transformation - is very small. Consequently, the 
spectral shape in Wigner space is determined by the finite integration interval in time rather 
than by the interactions itself. On the other hand, we have used an infinite relative time range 
in deriving the generalized transport equation within the first order gradient expansion. Thus 
in this case we deal with a completely resolved spectral function already at the initial time! 
This demonstrates why both approaches can only be compared to a limited extent for the very 
early times! 

Finally, concentrating on the very early time behaviour, we find a significant difference 
between the full and the approximate dynamics in the gradient scheme (75). For the genera- 
lized transport equation we observe a monotonous evolution of the equal-time Green function 
momentum modes whereas strong oscillations are seen in the initial phase for the solution of the 
full Kadanoff-Baym theory. Thus - with respect to the early time behaviour - the generalized 
transport equation behaves much more like the Boltzmann approximation, which is a first 
order differential equation in time as well. However, the Kadanoff-Baym evolution is given by 
an intcgro-differential equation of second order in time. In this case the phase correlations 
between the Green functions G^, G^, G< n are kept and the instantaneous switching-on of 
the interaction results in an oscillatory behaviour of the single momentum modes. 



2.2.2 Evolution of the Spectral Function 

Since the Green functions develop in time also the spectral properties of the system change as 
well. In Fig. 16 the time evolution of the spectral functions for the initializations Dl, D2 and 
D3 within the gradient scheme are displayed for two particular momentum modes | p \/m = 0.0 
(l.h.s.) and \ p\/m= 2.0 (r.h.s.) for various system times t ■ m = 5, 20, 60, 120, 240, 360, 480, 
600 up to the long-time limit. This representation corresponds to Fig. 8, where the respective 
evolution of the spectral function is studied for the full Kadanoff-Baym theory. We find that the 
time evolution of the spectral functions obtained from the generalized transport equation (75) 
is very similar to the one from the full quantum calculation (see below). The zero-mode spectral 
function for the initial distribution Dl becomes sharper with time and is moving to slightly 
higher energies. The opposite characteristics is observed for the zero-mode spectral function for 
the initialization D3, which broadens with time (reducing the peak correspondingly) and slowly 
shifts to smaller energies. Together with the weak evolution for the distribution D2 (which only 
slightly broadens at intermediate times and returns to a narrower shape at smaller energies in 
the long-time limit) the evolution of all three initializations in the scmiclassical approximation 
is well comparable to the full Kadanoff-Baym dynamics (cf. Fig. 8). Furthermore, the maxima of 
the zero-mode spectral functions are located above the bare mass (as indicated by the on-shell 
arrow) for all initial states during the time evolution. 

The spectral functions for the momentum mode |p|/m = 2.0 are in a good agreement 
with the Kadanoff-Baym dynamics as well. Again we observe - for the initial distribution 
Dl - a narrowing of the spectral function, while for D3 the spectral function broadens with 
time. Moreover, the width of the spectral function starting from distribution D2 shows a non- 
monotonous behaviour with a maximum at intermediate times. 
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Fig. 16. Time evolution of the spectral function A(p,po,t) for the initial distributions Dl, D2 and 
D3 (from top to bottom) with coupling constant A/m = 18 for the two momentum modes | p \ /m — 
0.0 (l.h.s.) and |pj/m — 0.0 (r.h.s.). The spectral function from the transport eq. (75) is shown for 
times t ■ m = 5, 20, 60, 120, 240, 360, 480, 600 as indicated by the different line types. 



In order to study the dynamics of the spectral function in a more quantitative manner 
we display in Fig. 17 the time evolution of the on-shell energies (as derived from the maxima 
of the spectral function) for the momentum modes | p |/m = 0.0 (upper plot) and 2.0 (lower 
plot) for the initializations Dl, D2 and D3 with A/m = 18 (l.h.s.). By comparison with the 
corresponding results from the Kadanoff-Baym theory (r.h.s.) we observe a close similarity of 
the evolutions within the full and the semiclassical KB scheme. The effective mass of the zero 
momentum mode decreases for initialization D3, passes a minimum for D2 and increases for 
the initial state Dl. 

As familiar from the Kadanoff-Baym calculations in Section 1.2 the behaviour of the on-shell 
energies is different for higher momentum modes. We find for the momentum mode | p \/m = 2.0 
a monotonous decrease of the on-shell energy for the initializations Dl and D2 and an increase 
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Fig. 17. Time evolution of the on-shell energies uj(p,t) of the momentum modes | p \/m = 0.0 and 
| p \Jm = 2.0 for the different initializations Dl, D2 and D3 with X/m = 18 in the semiclassical KB limit 
(75) (l.h.s.). The on-shell self-energies are extracted from the maxima of the time-dependent spectral 
functions. The respective results from the full Kadanof-Baym theory are displayed on the r.h.s. 
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Fig. 18. Time evolution of the on-shell widths —Im S R (p, uj(p, t),t)/u}(p, t) of the momentum modes 
p |/m = 0.0 and | p \/m = 2.0 for the different initializations Dl, D2 and D3 with X/m — 18 in the 
semiclassical transport eq. (75) (l.h.s.). The respective results from the full Kadanof-Baym theory are 
displayed on the r.h.s. 



for distribution D3. Altogether, the evolution of the on-shell energies for the higher modes 
is rather moderate compared to the lower ones in accordance with the dominant momentum 
contribution and the weakening of the retarded self-energy for higher energy modes. 

Finally, the on-shell energies approach a stationary state for all modes and all initializations. 
However, the long-time limit of the equal momentum modes is not exactly the same for all initial 
distributions Dl, D2 and D3. As discussed above this small difference can be traced back to 
the specific initial state generation from the given momentum distribution. 
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Next we consider the time evolution of the on-shell width as determined by the imaginary 
part of the retarded self-energy at the maximum position of the spectral function. In Fig. 18 
(l.h.s.) the on-shell width is displayed for the two momentum modes | p \/m = 0.0 and | p|/m = 
2.0 for all three initial distributions Dl, D2 and D3 with \/m —18 as a function of time. For both 
momentum modes the on-shell width increases for the distribution D3, while it has a maximum 
at intermediate times (t-m « 40) for the initialization D2. Thus the results - together with the 
reduction of the on-shell width for both momentum modes for the initialization Dl - is in good 
agreement with the results obtained for the full Kadanoff-Baym theory (r.h.s.). However, the 
stationary values for the on-shell widths deviate again slightly in accordance with the different 
preparation of the initial state in the gradient scheme. 



In summarizing we find that the main characteristics of the full quantum evolution of the 
spectral function are maintained in the semiclassical transport equation (75) as well. This 
includes the evolution of the on-shell energies as well as the width of the spectral function. 
Since the generalized transport equation is formulated directly in Wigner space one has access 
to the spectral properties at all times, whereas the very early times in the Kadanoff-Baym case 
have to be excluded due to the very limited support in the relative time interval (ti — £2) for 
the Wigner transformation. 
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Fig. 19. Spectral function A(po) for various momentum modes as a function of energy po/m at the 
late time t ■ m = 600 (lower part) for initial distribution D2 with coupling constant A/m = 18 in the 
semiclassical transport limit (75) (l.h.s.). The corresponding distribution function N (at the same time 
for the same momentum modes) is shown in the upper part. All momentum modes can be fitted for 
all energies by a single Bose function of temperature T/m = 1.827 and a chemical potential fi/m close 
to zero. The respective results from the full Kadanof-Baym theory are displayed on the r.h.s. 
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2.2.3 Stationary State of the Semiclassical Evolution 

As we have observed in the previous Subsections the evolution within the generalized trans- 
port equation (75) leads to a stationary state for all three different initializations Df , D2 and 
D3. Thus we turn to the investigation of this long-time limit itself, here in particular for the 
initialization D2. In Fig. 19 (l.h.s.) we show the distribution function N of various momentum 
modes |p|/m = 0.0, 0.8, 1.6, 2.4, 3.2, 4.0 for large times (t ■ m_= 60_0) as derived from the 
Green function itself and the spectral function via the relation N = iG < j A. The distribution 
function for a given momentum mode is calculated for all energies po where the corresponding 
spectral function - as displayed in the lower part of Fig. 19 - exceeds a value of 0.5. Since the 
width of the late time spectral function decreases with increasing momentum, the energy range 
(for which the distribution function is shown), is smaller for larger momentum modes. We find, 
that all momentum modes of N can be fitted at all energies by a single Bose function with 
a temperature T/m = 1.827 and a very small chemical potential /i/m = 0.000817. Thus the 
generalized transport formulation (75) leads to a complete (off-shell) equilibration of the system 
very similar to the solution of the full Kadanoff-Baym equation (r.h.s. of Fig. 19). Furthermore, 
the long-time limit of the semiclassical time evolution exhibits a vanishing chemical potential 
n/m in accordance with the properties of the neutral 4 -theory. This might have been expected 
since in the generalized transport equation particle number non-conserving processes of the 
type 1 <-» 3 — which lead to the decrease of the chemical potential - are included by means of 
the dynamical spectral function. Thus the semiclassical approximation (75) solves the problems 
within the Boltzmann limit, which does not yield a relaxation of the chemical potential, since 
only on-shell 2^2 transitions of quasiparticles are taken into account as demonstrated in 
Section 1.4. 

After observing, that the chemical potential decreases to zero in the long-time limit, it is 
interesting to study the relaxation process itself. The relaxation of the chemical potential /i/m 
for the three different initializations Dl, D2 and D3 (with coupling constant A/to = 18) shows 
an approximately exponential decrease in time. The relaxation rates - as determined from the 
slope of the exponential decline - are also approximately the same for all distributions. They 
are given explicitely by r£ 1 w 0.98 • 10~ 2 for distribution Dl, I^f 2 w 1.01 • 10~ 2 for distribution 
D2 and 3 w 1.07 • 10~ 2 for distribution D3. Thus the relaxation rates are in the same range 
as those found within the full Kadanoff-Baym theory (Section 1.3). 

We conclude that - although there is a small relative shift of the different time scales of 
kinetic and chemical equilibration as a function of the coupling strength A with respect to the 
full KB solutions - the results of the generalized KB transport equations are very similar. The 
differences we attribute to higher-order multi-particle effects in off-shell transitions. While the 
kinetic equilibration proceeds approximately with the coupling constant squared (as indicated 
by the calculations for non-polar-symmetric systems), the chemical relaxation rate is a higher 
order process in the coupling constant A. 

2.2.4 Quadrupole Relaxation 

In this Subsection we no longer restrict to polar symmetric systems and discuss the time 
evolution of more general initial distributions within the generalized transport approximation 
(75). We start with conditions similar to those employed in Section 1.4 but combined with the 
initialization scheme for the semiclassical limit. Again the decrease of the quadrupole moment 
of the distribution 




(86) 



is approximately exponential in time (oc exp(— Jq • t)) and thus allows for the extraction of a 
quadrupole damping rate Tq. The scaled quadrupole damping rates - as obtained for the two 
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initial distributions Dl and D2 - are displayed in Fig. 20 as a function of the coupling strength 
X/m. 
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Fig. 20. Scaled Relaxation rate for the generalized transport equation as a function of the coupling 
strength (half-filled symbols) for the initial distributions dl=Dl (squares) and d2=D2 (triangles). Ad- 
ditionally the results obtained within the Kadanoff-Baym (full symbols) and the Boltzmann calculation 
(open symbols) are shown for comparison. 



The calculations show that the quadrupole relaxation rates within the semiclassical approxi- 
mation (75) for both initial distributions Dl and D2 is well within in the range of the full 
Kadanoff-Baym and the on-shell Boltzmann case. Additionally, the quadrupole relaxation rate 
is rather flat in the coupling A when divided by the coupling constant squared (A 2 ) as already 
observed for the other two evolution schemes in Section 1.4. Nevertheless, the relaxation in the 
full KB-theory (8) proceeds slightly faster than in the transport limits for large couplings. The 
latter effect is again attributed to higher order off-shell transition effects which are no longer 
incorporated in the generalized KB transport equation. 



2.2.5 Kadanoff-Baym versus Botermans-Malfliet transport 

In this Subection we will perform a comparison of the generalized transport equation in the 
original Kadanoff-Baym (KB) form (75) with the modified Botermans-Malfliet (BM) form (81). 
As discussed in detail in Section 2.1.2 the latter form results from the replacement of the 
collisional self-energy by 

iS< -> iG< ■ r/A (87) 

in the second Poisson bracket on the l.h.s. of the original kinetic equation (75). This replacement 
leads to a consistent first order equation in the gradients and achieves consistency of the re- 
sulting transport equation with the corresponding generalized mass-shell relation (76). As noted 
before, the replacement (87) might be questionable in the general case; thus explicit numerical 
investigations are needed. 

In Fig. 21 we compare the time evolution within the generalized transport equation in the 
KB form to the consistent equation in BM form. In this respect several momentum modes of 
the equal-time Green function are displayed evolving in time from an initial distribution D2 for 
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Fig. 21. Time evolution of various equal-time Green function momentum modes within the generalized 
transport equation (original KB form, solid black lines) and within the Botermans-Malfiet (BM) form 
(dashed red lines) for initial distribution D2 with coupling constant A/m = 16. 



a coupling constant A/m = 16. We find that the deviations between both approximations (KB 
and BM) are rather moderate. Only for very small momentum modes | p |/m < 1.6 deviations 
between both modes are visible. For the very low momentum modes the range of difference 
starts at t ■ m ss 10 and extends to t ■ m « 100 for the non-zero modes. For the zero momentum 
mode the deviation lasts even longer. In this region the semiclassical transport in the BM form 
is slightly 'slower' than in the original KB choice. Nevertheless, also the BM form exhibits the 
typical overshooting behaviour of the low momentum modes beyond the stationary limit as 
observed for the KB form. However, the maxima are shifted slightly to later times. Finally, 
both gradient approximations converge in the long-time limit to very similar configurations. 

In summarizing we point out that the approximation of the full Kadanoff-Baym dynamics 
by the generalized transport equations in Kadanoff-Baym (75) or Botermans-Malfliet form (81) 
holds very well for the different momentum modes of the Green function iG < itself. Slight 
deviations are only visible for the zero momentum mode at early to intermediate times (Figs. 
15 and 21) for a logarithmic representation of the time axis. Consequently, the characteristic 
features of quantum equilibration obtained for the full Kadanoff-Baym theory are retained in 
the generalized transport limits. 



2.3 Testparticle representation 



The generalized transport equation (81) allows to extend the traditional on-shell transport 
approaches for which efficient numerical recepies have been set up [29, 78-91] (and Refs. therein). 
In order to obtain a practical solution to the transport equation (81) we use a testparticle ansatz 
for the Green function G < , more specifically for the real and positive scmidefinite quantity 



"XP 



iG<(X,P) ~ * (3) ( X - M*)) ^ (3) (P - Pi(*)) S(P -et(t)) 



(88) 



i=i 



In the most general case (where the self energies depend on four-momentum P, time t and the 
spatial coordinates X) the equations of motion for the testparticles read [65] 



dt 



1 - C(i) 2ti 



2 Pi + V Pl ReE\ 



Ml 
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(89) 
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dt 



1-C, 
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V x , ReE? 
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(90) 
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(*) 2ei 
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9< 



r, 



at 



(91) 



where the notation i 7 ^) implies that the function is taken at the coordinates of the testparticle, 



i.e. F, 



F(t,Xi(t),Pi(t),ei(t)). 



In (89-91) a common multiplication factor (1 — C(i)) 1 appears, which contains the energy 
derivatives of the retarded self energy 



C (i ) - 



2e, : 



e? - P? - M 2 - Re£[$ q 



r, 



« 



(92) 



It yields a shift of the system time t to the 'eigentime' of particle i defined by t% = t/(l — Cu)). 
As the reader immediately verifies, the derivatives with respect to the 'eigentime', i.e. rfXj/rffj, 
dPi/dti and dei/dti then emerge without this renormalization factor for each testparticle i when 
neglecting higher order time derivatives in line with the semiclassical approximation scheme. 

Some limiting cases should be mentioned explicitly: In case of a momentum-independent 
'width' r(X) we take M 2 = P 2 — ReS ret as an independent variable instead of Po, which then 
fixes the energy (for given P and M 2 ) to 



p2 

o 



= P 2 + M 2 + ReE r x e *, M2 
Eq. (91) then turns to (AM 2 = M 2 - M 2 ) 
dAM 2 _ AM 2 dr {i) 



dt 



r, 



(i) 



dt 



d i 
— In 

dt 



(AM 2 



= 



(93) 



(94) 



for the time evolution of the test-particle i in the invariant mass squared. In case of r = const. 
the familiar equations of motion for testparticles in on-shell transport approaches are regained. 



2.3.1 Model simulations for the momentum-dependent transport equations in testparticle 
representation 

To demonstrate the physical content of the equations of motion for testparticles (89-91) we 
perform an exploratory study with a momentum-dependent trial potential. The potential is 
chosen of the type: 

2 l + exp{(|r| -i?)/a } ^ 1 + exp{(|r| - -R)/a } 2 / k ' 

with a constant (but finite) vacuum width Py . While the spatial extension of the potential is 
given by a Woods-Saxon shape (with parameters R = 5 fm and do = 0.6 fm) its momentum 
dependence for the real as well as for the imaginary part is introduced by 

V(Po,P) =Cv Ai _$_ p2) , W(P ,P) =CW ^_^_ pa) - (96) 

Here the constants Cy (C\y) give the 'strength' of the complex potential while Ay (A w ) play 
the role of cutoff-parameters. Due to the structure in the denominator of (96) the momentum- 
dependent part of this potential is explicitly Lorentz-covariant but should be considered for 
M 2 « A 2 N . 
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In the simulations the testparticles are propagated with different initial mass parameters 
Mi, which are shifted relative to each other by TV/20 [GeV] around a mean mass of 1.0 GeV. 
To each testparticle a momentum in positive ^-direction is attributed so that all of them have 
initially the same energy Pq — 2.0 GeV. All testparticles are initialized on the negative z-axis 
with (|Xj(£ = 0)| i=a 15 fm) and then evolved in time according to the equations of motion 
(89-91). 
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Fig. 22. L.h.s.: upper part: P;o, Piz and Mi as a function of z(t) for a momentum-dependent imaginary 
potential with Cw = 0.6 GeV 2 and Aw = 2.0 GeV (lower part). The vacuum width is chosen as 
IV = 0.02 GeV 2 and the initial separation in the mass parameter of the testparticles is AM — rv/20. 
R.h.s.: upper part: Pio, Piz and Mi as a function of z(t) for a momentum-dependent complex potential 
with C v = -0.3 GeV 2 , A v = 2.0 GeV, C w = 0.6 GeV 2 and A w = 2.0 GeV (lower part). 



In our first simulation we consider a purely imaginary potential with a strength of Cw = 0.6 
GeV 2 and a cutoff-parameter Aw = 2.0 GeV. The evolution in energy P oi , momentum P zi and 
in the mass parameter Mi for all testparticles is shown in Fig. 22 (l.h.s., upper part) as a 
function of z(t). When the testparticles enter the potential region, their momenta and mass 
parameters are modified. The imaginary potential leads to a spreading of the trajectories in the 
mass parameter Mj which in turn reflects a broadening of the spectral function. The relation 
between the imaginary self energy and the spreading in mass is fully determined by relation 
(94). Since we have chosen a potential with no explicit time dependence the energy of each 
testparticle is a constant of time. According to the explicit momentum dependence of our 
'trial' potential each single testparticle is affected with different strength. Since the imaginary 
potential is strongest for small momenta (which correspond to the highest lines in the lower 
graph of Fig. 22) (l.h.s.) the momentum and mass coordinates of those testparticles are changed 
predominantly that are initialized with the lowest momenta (i.e. with the largest masses). As 
a result one observes a rather asymmetric distribution in the mass parameters (and in the 
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momenta) in the potential zone. For z(t) 3> R the mass and momentum coordinates of the 
testparticles regain the proper asymptotic value. 

In the second example we allow for an addititonal real part of the self energy. The calculation 
is performed with the parameters Cy = —0.3 GeV 2 , Cw = 0.6 GeV 2 and Ay = Aw = 2.0 
GeV. The momentum-dependent real part (r.h.s., lower part of Fig. 22) causes an additional 
shift of the testparticle momenta. Since the real part of the potential is larger for small initial 
momenta, these testparticle momenta are shifted up somewhat more than for particles with 
larger momenta. This gives rise to a reduction of the asymmetry which was introduced by the 
momentum-dependent imaginary part of the self energy (r.h.s., upper part of Fig. 22). The mass 
parameters of the testparticles are only weakly influenced by the real part of the potential. 



> 

O 



2.0 
1.9 
1.8 
1.7 
1.6 
1.5 
1.4 
1.3 



-i 1 . 1 . 1 1 1 1 1 1- 




-15 



-10 




'correction' 
term 

on 

off 



_j , i i i i u 



10 15 

z[fm] 




'correction term' 



10 15 

z[fm] 

Fig. 23. Lower part: Correction term (1 — C"^)) -1 as a function of z(t) for the same imaginary potential 
as in Fig. 22. Upper part: P iz as a function of z(t) with and without including the correction term. 
Both curves are identical. 



We, furthermore, study the implications of the correction term (1 — C^))" 1 using the same 
imaginary potential as in Fig. 22. The time evolution of the correction factor for each testparticle 
i is displayed in the lower part of Fig 23. While it is > 1 for large initial mass parameters 
Mi > Mo, it is < 1 for mass parameters Mi < Mo. In the upper part of Fig. 23 the momenta 
of the testparticles are shown for two calculational limits: in the first one the correction term is 
taken into account (as in the previous calculations) , while in the second one the corrections due 
to the energy dependence of the retarded self energy are neglected. However, the calculations 
with and without the correction factor exhibit no difference in the testparticle momenta as a 
function of z. The same holds for the mass parameters Mi which are not displayed here since 
they provide no new information due to energy conservation. We thus conclude that the particle 
trajectory (in phase space) is independent of these correction factors, since the correction term 
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- when displayed as P z (z) in phase-space - leads only to a rescaling of the 'eigentime' of the 
testparticles as pointed out above. 



2.3.2 Collision terms 

The collision term of the Kadanoff-Baym equation can only be worked out in more detail by 
giving explicit approximations for S < and S > . A corresponding collision term can be formu- 
lated in full analogy to Refs. [79,92], e.g. from Dirac-Brucckner theory (cf. also Section 1.2), 
following detailed balance as 

I coll (X, P, M 2 ) = Tr 2 Tr 3 Tr 4 A(X, P, M 2 )A(X, P 2 ,M 2 )A(X, V 3 ,M 2 )A{X, P 4 , Mf ) 
|T((P, M 2 ) + (P 2 , Ml) (P 3 , Ml) + (P 4 , Ml))\\ s JW(P + P 2 - P 3 - P 4 ) 

NxP 3 M$ N XP 4 M% fxPM 2 fxP 2 M 2 - N X PM 2 N XP 2 M 2 fxP 3 M 2 /xP 4 M| (97) 

with 

Ixpm 2 = 1 + V Nxpmi (98) 

and 77 = ±1 for bosons/fcrmions, respectively. The indices ^4, <S stand for the antisymme- 
tric/symmetric matrix element of the in- medium scattering amplitude T in case of fcrmions/bo- 
sons. In Eq. (97) the trace over particles 2,3,4 reads explicitly for fermions 

where a 2 , t 2 denote the spin and isospin of particle 2. In case of bosons we have 

since here the spectral function is normalized as 

/" r/P 2 

J^-A B (X,P) = 1 (101) 



whereas for fermions we have 

^-A F (X.P) = I. (102) 



/ 



We mention that the spectral function Ap in case of fermions in (97) is obtained by considering 
only particles of positive energy and assuming the spectral function to be identical for spin 'up' 
and 'down' states. In general, the spectral function for fermions A a p(X, P) is a Dirac-tensor 
with af3 denoting the Dirac indices. It is normalized as 



/ 



dP 

-^A af3 (X,P) = ( 7 °) a0 , (103) 



which implies 

\j2[^° A (X,P)) aa = l. (104) 

a 

Now expanding A in terms of free spinors u s (P) (s=l,2) and v s (P) as 

2 

{A) afj = J2 u s {VM) (j AP S u r (P, M) a + v s {PM)jj A%% v r (P, M) a (105) 

r,s—l 
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one can separate particles and antiparticles. By neglection of the antiparticle contributions (i.e. 
Afv = 0) and within the assumption that the spectral function for the particles is diagonal in 
spin-space (i.e. AP S = 8 ra A p s ) as well as spin symmetric, one can define Ap as 

A F = A\ = A\ . (106) 

Neglecting the 'gain-term' in eq. (97) one recognizes that the collisional width of the particle 
in the rest frame is given by 

r coll (X, P, M 2 ) = Tr 2 Tr 3 Tr 4 |T((P, M 2 ) + (P 2 , M|) - (P 3 , M 3 2 ) + (P 4 , M 2 ))\ 2 A S (107) 
A(X, P 2 ,M 2 )A(X, P 3 ,M$)A(X, P 4 , Ml) 5 4 (P + P 2 - P 3 - P 4 ) N XP2M , f XVatq f XPtM , , 

where - as in eq. (97) - local on-shell scattering processes are assumed for the transitions 
P + P 2 — » P3 + Pa- We note that the extension of eq. (97) to inelastic scattering processes 
(e.g. NN — ► NA) or (irN — » Z\ etc.) is straightforward when exchanging the elastic transition 
amplitude T by the corresponding inelastic one and taking care of Pauli-blocking or Bose- 
enhanccment for the particles in the final state. 

For particles of infinite life time in vacuum - such as protons - the collisional width (107) has 
to be identified with twice the imaginary part of the self energy. Thus the transport approach 
determines the particle spectral function dynamically via (107) for all hadrons if the in-medium 
transition amplitudes T are known in their full off-shell dependence. Since this information is 
not available for configurations of hot and dense matter, which is the major subject of future 
development, a couple of assumptions and numerical approximation schemes have to be invoked 
in actual applications. 
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3 The dynamical quasiparticle approach to hot QCD 

This Section addresses the question how to achieve a suitable approximation for the spectral 
functions of partons (quarks and gluons) in hot QCD, i.e. at temperatures above the dcconfinc- 
ment transition T c . Since the spectral functions are fully determined by complex selfenergies 
U ret for the partons, it is sufficient to determine effective selfenergies, e.g., from lattice QCD. 



3.1 Application to parton dynamics 

The 'Big Bang' scenario implies that in the first micro-seconds of the universe the entire system 
has emerged from a partonic system of quarks, antiquarks and gluons - a quark-gluon plasma 
(QGP) - to color neutral hadronic matter consisting of interacting hadronic states (and reso- 
nances) in which the partonic degrees of freedom are confined. The nature of confinement and 
the dynamics of this phase transition has motivated a large community for several decades (cf. 
[93,94] and Refs. therein). Early concepts of the QGP were guided by the idea of a weakly 
interacting system of partons since the entropy density s and energy density e were found in 
lattice QCD to be close to the Stefan Boltzmann (SB) limit for a relativistic noninteracting 
system [95]. However, experimental observations at the Relativistic Heavy Ion Collider (RHIC) 
indicated that the new medium created in ultrarelativistic Au+Au collisions was interacting 
more strongly than hadronic matter and consequently this notion had to be given up. More- 
over, in line with earlier theoretical studies in Refs. [96-98] the medium showed phenomena of 
an almost perfect liquid of partons [99,100] as extracted from the strong radial expansion and 
elliptic flow of hadrons as well the scaling of the elliptic flow with parton number etc. All the 
latter collective observables have been severely underestimated in conventional string/hadron 
transport models [101-103] whereas hydrodynamical approaches did quite well in describing (at 
midrapidity) the collective properties of the medium generated during the early times for low 
and moderate transverse momenta [104,105]. The question about the constituents and prop- 
erties of this QGP liquid is discussed controversely in the literature (cf. Refs. [106-108]) and 
practically no dynamical concepts are available to describe the dynamical freezeout of partons 
to color neutral hadrons that are finally observed experimentally. Since the partonic system ap- 
pears to interact more strongly than even hadronic systems the notation strong QGP (sQGP) 
has been introduced in order to distinguish from the dynamics known from perturbative QCD 
(pQCD). 

Lattice QCD (1QCD) calculations provide some guidance to the thermodynamic properties 
of the partonic medium close to the transition at a critical temperature T c up to a few times 
T c , but 1QCD calculations for transport coefficients presently arc not accurate enough [109] 
to allow for firm conclusions. Furthermore, it is not clear whether the partonic system really 
reaches thermal and chemical equilibrium in ultrarelativistic nucleus-nucleus collisions [110] 
such that nonequilibrium models are needed to trace the entire collision history. The available 
string/hadron transport models [81,82,111] partly fail - as pointed out above - nor do partonic 
cascade simulations [112-115] (propagating massless partons) sufficiently describe the reaction 
dynamics when employing cross sections from perturbative QCD. Some models, e.g. the Multi- 
phase Transport Model AMPT [116], employ strong enhancement factors for the cross sections, 
however, use only on-shell massless partons in the partonic phase as in Ref. [113]. The same 
problem comes about in the parton cascade model of Ref. [117] where additional 2^-> 3 processes 
like gg <-» ggg are incorporated but massless partons are involved. 

On the other hand it is well known that strongly interacting quantum systems require 
descriptions in terms of propagators with sizeable selfenergies for the relevant degrees of freedom 
(cf. Sections 1 and 2). While the real part of the selfenergies can be related to mean-field 
potentials, the imaginary parts provide information about the lifetime and/or reaction rate of 
time-like 'particles' [118]. The studies of Peshier [119,120] indicate that the effective degrees of 
freedom in a partonic phase should have a width 7 in the order of the pole mass M already 
slightly above T c . This opens up the problem how to interpret/deal with the space- like part of 
the distribution functions and how to 'propagate' effective degrees in space-time in equilibrium 
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as well as out of equilibrium without violating 'microcausality'. These questions will now be 
addressed in the Dynamical QuasiParticle Model (DQPM) [121]. 



3.2 Off-shell elements in the DQPM 

The dynamical quasiparticle model goes back to Pcshier [119,120] and starts with the entropy 
density in the quasiparticle limit [120,122-124], 



where n B (u/T) = (exp(w/T) - 1) _1 and n F ((w - fJ, q )/T) = (exp((w - /J. q )/T) + 1) _1 denote the 
Bose and Fermi distribution functions, respectively, while A = (P 2 — 77) _1 , S q = (P 2 — Sq)^ 1 
and Sq = (P 2 — Sq)^ 1 stand for the full (scalar) quasiparticle propagators of gluons g, quarks 
q and antiquarks q. The degeneracy factor for gluons is d g — 2(N 2 — 1) = 16 while for quarks 
q and antiquarks q we get d q = dq = 2N c Nf = 18 for three flavors Nf. In Eq. (108) 77 and 
S = S q ;=y Sq denote the (retarded) quasiparticle selfenergies. In principle, 77 as well as A are 
Lorentz tensors and should be evaluated in a nonperturbative framework. The DQPM treats 
these degrees of freedom as independent scalar fields with scalar selfenergies. In case of the 
fermions S q ,Sq and S q ,Sq (for q and q) have Lorentz scalar and vector contributions but 
only scalar terms are kept in (108) for simplicity which are assumed to be identical for quarks 
and antiquarks. Note that one has to treat quarks and antiquarks separately in (108) as their 
abundance differs at finite quark chemical potential fj, q . 

Since the nonperturbative evaluation of the propagators and selfenergies in QCD is a 
formidable task (and addressed in Dyson-Schwingcr (DS) Bethe-Salpeter (BS) approaches 
[125-129]) an alternative and practical procedure is to use physically motivated Ansdtze with 
Lorentzian spectral functions for quarks 2 and gluons, 



and to fit the few parameters to results from 1QCD. With the convention E 2 (p) — p 2 + 
M 2 — 7 2 , the parameters M 2 and 7 are directly related to the real and imaginary parts of the 
corresponding (retarded) self-energy, e.g. 77 = M 2 — 2i^/uu in case of the 'scalar' gluons. 

Following [120] the quasiparticle mass (squared) for gluons is assumed to be given by the 
thermal mass in the asymptotic high-momentum regime, i.e. 




(108) 






(109) 




(110) 



and for quarks (assuming vanishing constituent masses here) as, 




(111) 



2 In the following the abbreviation is used that 'quarks' denote quarks and antiquarks if not specified 
explicitly. 
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with a running coupling (squared), 

2/ 

9 {T/Tc) = (UN C - 2N f ) ln(A 2 (T/T c - T s /T c f 1 {U2) 

which permits for an enhancement near T c [130-132]. Here N c = 3 stands for the number 
of colors while Nf denotes the number of flavors and [i q the quark chemical potentials. The 
parameters A = 2.42 and T s /T c = 0.46 have been fixed in Ref. [120]. As demonstrated in Fig. 
24 this functional form for the strong coupling a s = g 2 /(An) is in accordance with the 1QCD 
calculations of Ref. [133] for the long range part of the q — q potential. 




0.0 1 i , I 

1 2 3 456789 10 
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Fig. 24. The coupling a s (T) = g 2 (T)/(4n) (solid red line) as a function of T/T c in comparison to 
the long range part of the strong coupling as extracted from Ref. [133] from the free energy of a 
quark-antiquark pair in quenched 1QCD (for 7V T = 8). 



The width for gluons and quarks (for fj, q — 0) is adopted in the form [134] 

g 2 T 2c N 2 - 1 g 2 T 2c 

7 s CO = iV c — In—, ^ = — j ——. hl -. 

where c = 14.4 (from Ref. [97]) is related to a magnetic cut-off. 

The physical processes contributing to the width j g are both gg <-» gg, gq <-> gq scattering 
as well as splitting and fusion reactions gg <-> g, gg <-> ggg, ggg «-> gggg or g «-> qq etc. On 
the fermion side elastic fermion-fermion scattering pp <-> pp, where p stands for a quark q or 
antiquark q, fermion-gluon scattering pg <-> pg, gluon bremsstrahlung pp <^ pp + g ot quark- 
antiquark fusion qq ^ g etc. emerge. Note, however, that the explicit form of (113) is derived 
for hard two-body scatterings only. It is worth to point out that the ratio of the masses to their 
widths ~ gln(2c/g 2 ) approaches zero only asymptotically for T — * oo such that the width of 
the quasiparticles is comparable to the pole mass slightly above T c up to all terrestrial energy 
scales. 

Within the DQPM the real and imaginary parts of the propagators A and 5 now are 
fixed and the entropy density (108) can be evaluated numerically once the free parameters in 
(112) are determined. In the following we will assume 3 light quark flavors Nf = 3. Since the 
presently available unquenched 1QCD calculations (for three flavors) for the entropy density 
are still accompanied with rather large error bars the parameters of the DQPM are taken the 
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same as in the pure Yang-Mills sector: A = 2.42, T s /T c — 0.46 as determined in Refs. [97, 
120]. This is legitimate since an approximate scaling of thermodynamic quantities from 1QCD 
is observed when dividing by the number of degrees of freedom and scaling by the individual 
critical temperature T c which is a function of the different number of parton species [135]. 
However, these parameters will have to be refitted once more accurate 'lattice data' become 
available. 

The resulting values for the gluon and quark masses - multiplied by T c /T - are displayed in 
Fig. 25 (for \i q = 0) by the solid lines while the gluon and quark width (p/ g , j q ) - multiplied by 
T c /T - are displayed in terms of the dashed lines as a function of T/T c . The actual numbers 
imply very 'broad' quasiparticles already slightly above T c . For [i q = the ratio 7 9 /7 g =4/9 
is the same as for the ratio of the squared masses m 2 /M 2 = 4/9 and reflects the ratio of the 
Casimir eigenvalues in color space. Consequently the ratio of the width to the pole mass is 
smaller for quarks (antiquarks) than for gluons in the whole temperature range. 
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Fig. 25. The mass M (solid red line) and width 7 9 (dashed blue line) for gluons (l.h.s.) and the mass 
m (solid orange line) and width 7 q (dashed green line) for quarks (r.h.s.) as a function of T/T c in the 
DQPM for A = 2.42, T a /T c = 0.46, and c = 14.4 (for fi q = 0). All quantities have been multiplied by 
the dimensionless factor T c /T. 



In order to fix the scale T c , which is not specified so far, one may directly address unquenched 
1QCD calculations (for 3 light flavors). However, here the situation is presently controversal 
between different groups (cf. Refs. [136,137] and the discussion therein). An alternative way is 
to calculate the pressure P from the thermodynamical relation (at ji q — 0), 




(114) 



by integration of the entropy density s over T, where one may tacitly identify the 'full' entropy 
density s with the quasiparticle entropy density s dqp (108). Since for T <T C the DQPM entropy 
density drops to zero (with decreasing T) due to the high quasiparticle masses and the width 7 
vanishes as well (cf. Fig. 25) the integration constant may be assumed to be zero in the DQPM 
which focusses on the quasiparticle properties above T c . 



The energy density e then follows from the thermodynamical relation [130,138] 

e = Ts-P (115) 
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Fig. 26. The DQPM results for e(T/T c )(T c /T) 4 (solid violet line), P(T/T C )(T C /T) 4 (blue dash-dotted 
line) and the interaction measure W(T/T C )(T C /T) 4 (116) (green dotted line). Note the logarithmic 
scale in T/T c . The energy density e (115) practically coincides with the quasiparticle energy density 
T 00 from (123) (dashed red line). 



(for Hq = 0) and thus is also fixed by the entropy s(T) as well as the interaction measure 

W(T) :— e(T) ~ 3P(T) = Ts — 4P (116) 

that vanishes for massless and noninteracting degrees of freedom. 

The actual results for e ■ (T c /T) 4 are displayed in Fig. 26 (solid violet line), P ■ (T c /T) 4 
(blue dash-dotted line) as well as the interaction measure W ■ (T c /T) 4 (116) (green dotted 
line) and show the typical pattern from 1QCD calculations [135]. The scale T c may now be 
fixed (estimated) by requiring that the critical energy density e(T c ) is roughly the same for the 
pure Yang-Mills case as for the full theory with dynamical quarks in line with the approximate 
scaling of 1QCD [135]. Since e(T c ) in the Yang Mills sector is about 1 GeV/fm 3 (cf. [121]) the 
scaled energy density e(T/T c ) from Fig. 26 can be employed to fix the critical temperature T c w 
0.185 GeV. This leads to the 'thumb rule' e « 2 (T/T c ) 4 [GeV/fm 3 ] for T > 1.2T C which is 
roughly fulfilled according to Fig. 26. 



A first test of the DQPM in comparison to 1QCD calculations is given for the sound velocity 
(squared), 




(117) 



which does not depend on the absolute (uncertain) scale of T c . This comparison is shown in Fig. 
27 by the solid line where the 1QCD results have been taken from Ref. [139] and correspond 
to Nf = 2 + 1 with N t — 6. The DQPM is seen to reproduce the drop in c 2 s close to T c 
within errorbars and to reach the asymptotic value c 2 s = 1/3 approximately for T > 2 T c . This 
comparison, however, has to be taken with some care since the present DQPM assumes massless 
current quarks whereas the 1QCD calculations employ finite quark masses. 




3.2.1 Time-like and space-like quantities 

For the further analysis of the DQPM it is useful to introduce the shorthand notations 

TrJ ■■■ = d 9 J^-0f M n B^/T) 0{±P 2 ) • • • (118) 

Trf • ■ ■ = d q J^ ^ 2cu p q {u>) 9(w) n F ((u> - ^)/T) 0(±P 2 ) ■ ■ ■ 

*«£•■• = ^ p&) n ^ + ^)/ T ) °( ±p2 ) ■ ■ ■ 

with P 2 = uj 2 — p 2 denoting the invariant mass squared. The quark and antiquark degrees 
of freedom, i.e. 2d q = 36, are by roughly a factor of two more abundant than the gluonic 
degrees of freedom. The <9(±P 2 ) function in (118) separates time-like quantities from space- 
like quantities and can be inserted for any observable of interest. Note, however, that not all 
space-like quantities have a direct physical interpretation! 

We note in passing that the entropy density (108) is dominated by the time- like contributions 
for quarks and gluons and shows only minor space-like parts (cf. Refs. [120, 121]). Furthermore, 
the entropy density from the DQPM is only 10 - 15 % smaller than the Stefan Boltzmann 
entropy density ssb for T > 2 T c as in the case of 1QCD [95]. 

Further quantities of interest are the 'quasiparticle densities' 

JV±(T) - Tr* 1, W±(T) = fr± 1, iV±(T) = Tr± 1, (119) 

that correspond to the time-like (+) and space-like (-) parts of the integrated distribution func- 
tions. Note that only the time- like integrals over space have a particle number interpretation. 
In QED this corresponds e.g. to time-like photons (7*) which are virtual in intermediate pro- 
cesses but may also be seen asymptotically by dileptons (e.g. e + e _ pairs) due to the decay 
7* -> e+ + e (fi + + n~) [111,140] (cf. Section 4). 
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Scalar densities for quarks and gluons - only defined in the time-like sector - are given by 



N s (T)=fv: 



V s (T)=Tr c 



, JVJ(T) = fr+ 



(120) 



and have the virtue of being Lorentz invariant. 

Before coming to the actual results for the quantities (119) and (120) it is instructive to 
have a look at the integrand in the quark density (119) which reads as (in spherical momentum 
coordinates with angular degrees of freedom integrated out) 



I(v,p) = Y% P 2 u p q (uj,p 2 )n F ((uj - fJ, q )/T) . 



(121) 



Here the integration is to be taken over u> and p from to oo. The integrand I(u),p) is shown in 
Fig. 28 for T = 1.05T C (l.h.s.) and T = 3T C (r.h.s.) (fi q = 0) in terms of contour lines spanning 
both one order of magnitude. For the lower temperature the quark mass is about 0.55 GeV and 
the width 7 as 0.034 GeV such that the quasiparticle properties are close to an on-shell particle. 
In this case the integrand I{w,p) is essentially located in the time-like sector and the integral 
over the space-like sector is almost negligible. This situation changes for T = 3T C where the 
mass is about 0.7 GeV while the width increases to 7 ~ 0.25 GeV. This situation is close to the 
systems studied in Sections 1 and 2 in case of (/> 4 -theory and a strong coupling (A/to « 18). As 
one observes from the r.h.s. of Fig. 28 the maximum of the integrand is shifted towards the line 
lu = p and higher momentum due to the increase in temperature by about a factor of three; 
furthermore, the distribution reaches far out in the space-like sector due to the Fermi factor 
rip {uj/T) which favors small lu. Thus the relative importance of the time-like (+) part to the 
space-like (-) part is dominantly controlled by the width 7 - relative to the pole mass - which 
determines the fraction of N~ with negative invariant mass squared (P 2 < 0) relative to the 
time-like part N+ (P 2 > 0). 
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Fig. 28. The integrand I(iv,p) (121) as a function of u and p for quarks at temperatures T = 1.05 T c 
(l.h.s.) and T = 3 T c (r.h.s.). At T — 1.05 T c (l.h.s.) the quasiparticle mass amounts to m ~ 0.55 GeV 
and the width to j q w 0.034 GeV while at T = 3 T c (r.h.s.) m w 0.7 GeV and -y q w 0.25 GeV. The 
contour lines in both figures extend over one order of magnitude. Note that for a convergence of the 
integrals (119) the upper limits for lo and p have to be increased by roughly an order of magnitude 
compared to the area shown in the figure. 



The actual results for the different 'densities' (multiplied by (T c /T) 3 ) are displayed in Fig. 
29 for gluons (l.h.s.) and quarks (r.h.s.) including the antiquarks. The lower (magenta) lines 
represent the scalar densities N s , the red solid lines the time-like densities N + , the green lines 
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the quantities N~ while the thick solid blue lines are the sum N = N + + N~ as a function of 
T/T c (assuming T c — 0.185 GeV). It is seen that N + is substantially smaller than N~ in case 
of gluons in the temperature range 1.1 < T/T c < 10. The quantity N follows closely the Stefan 
Boltzmann limit N$b for a massless noninteracting system of bosons which is given in Fig. 
29 (l.h.s.) by the upper dash-dotted line. Though TV differs by less than 20% from the Stefan 
Boltzmann (SB) limit for T > 2Tc the physical interpretation is essentially different! Whereas 
in the SB limit all gluons move on the light cone without interactions only a small fraction of 
gluons can be attributed to quasiparticles with density N + within the DQPM that propagate 
within the lightcone. The space-like part N~ corresponds to 'gluons' exchanged in t-channel 
scattering processes and thus cannot be propagated explicitly in off-shell transport approaches 
without violating causality and/or Lorentz invariance. In case of quarks (or antiquarks) the 
results are qualitatively similar but now the time-like part N + comes closer to the space-like 
part N~ since the ratio of the width to the pole mass (jq/m) is smaller than the corresponding 
ratio for gluons as stated above. Furthermore, the quantity N is closer to the respective SB 
limit (for massless fermions) due to the lower effective mass of the quarks. In this respect the 
quarks and antiquarks are closer to (but still far from) the massless on-shell quasiparticle limit. 




1 2 345678 910 1 2 345678 910 

T/T c T/T c 

Fig. 29. The various 'densities' (119) for gluons (l.h.s.) and quarks (r.h.s.) including antiquarks (at 
fi q = 0). The lower magenta lines represent the scalar densities N" , the red solid lines the time- 
like densities AT , the green lines the quantities N~ while the thick solid blue lines are the sum 
iV = N + + N~ as a function of T/T c . The upper dash-dotted lines display the Stefan Boltzmann limits 
Nsb for reference. All densities are multiplied by the dimensionless factor (T c /T) 3 to divide out the 
leading scaling with temperature. 

The scalar densities N s (lower magenta lines) follow smoothly the time-like densities N + 
(for gluons as well as quarks+antiquarks) as a function of temperature and uniquely relate to 
the corresponding time-like densities N + or the temperature T in thermal equilibrium. 

The separation of N + and N~ so far has no direct dynamical implications except for the 
fact that only the fraction N + can explicitly be propagated in transport models as argued 
above. Following Ref. [121] we, furthermore, consider the energy densities, 

T± iB (T) = Tr± u> , (122) 

that specify time-like and space-like contributions to the quasiparticle energy densities (x = 

g,q,q)- 

The result for the quasiparticle energy densities and T 00 are displayed in Fig. 30 for 
gluons (l.h.s.) and quarks+antiquarks (r.h.s.) as a function of T/T c . All quantities have been 
multiplied by the dimensionless factor (T c /T) in order to divide out the leading temperature 
dependence. The lower red solid lines show the time-like components T^q while the dashed 
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Fig. 30. The time-like energy density Tqq (red solid line), the space- like energy density T^ (dashed 
green line) and their sum Too = 2qq + T^ (upper blue line) as a function of T/T c for gluons (l.h.s.) 
and quarks (+ antiquarks) (r.h.s.). All densities are multiplied by the dimensionless factor (T c /T) 4 . 



green lines display the space-like parts T 00 which dominate over the time-like parts except 
in the vicinity of T c . The general behaviour of the scaled energy densities TJ^j is similar as 
for the 'densities' N^ 1 given in Fig. 29 since the extra factor w in the integrand does not 
change significantly the time-like and space-like parts. As in Ref. [121] the space-like parts are 
interpreted as potential energy densities while the time-like fractions are the gluon and quark 
quasiparticle contributions which propagate within the lightcone. 

Summing up the time-like and space-like contributions for gluons, quarks and antiquarks 
we obtain the total energy density T 00 , 

yoo _ T QQ g + T QO g + T^ q + T 0Q q + Tqq - + T Q0 - , (123) 

which is displayed in Fig. 26 by the dashed red line (scaled by T c /T) 4 ). As in the case of the 
pure Yang-Mills system in Ref. [121] the quantity T 00 practically coincides with the energy 
density e (115) obtained from the thermodynamical relations. Small differences of less than 5% 
show up which indicates that the DQPM in its present formulation is not fully consistent in the 
thermodynamical sense. Since these differences are small on an absolute scale and significantly 
smaller than differences between present independent 1QCD calculations for 3 quark flavors one 
may consider T 00 (T) ss e(T) and separate the kinetic energy densities from the potential 
energy densities as a function of the temperature T or - in equilibrium - as a function of 
the scalar densities N s or time- like densities N + , respectively. 

It is instructive to show the potential energies per degree of freedom V gg /N+ = Tq g /N+ 
and V qq /N+ = Tqq^/N^ as a function of T/T c . The corresponding quantities are displayed in 
Fig. 31 (l.h.s.) multiplied by T c /T in terms of the solid red line and the dot-dashed blue line. It is 
seen that the potential energies per degree of freedom steeply rise in the vicinity of T c and then 
increase approximately linear with temperature T. As expected from the larger width of the 
gluons the latter also show a potential energy per degree of freedom which is roughly a factor of 
two larger than the corresponding quantity for quarks (antiquarks). Consequently rapid changes 
in the temperature (or density) - as in the expansion of the fireball in ultrarelativistic nucleus- 
nucleus collisions - are accompanied by a dramatic change in the potential energy density and 
thus to a violent acceleration of the quasiparticles. It is speculated here that the large collective 
flow of practically all hadrons seen at RHIC [99] might be attributed to the early strong partonic 
forces expected from the DQPM. 
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Fig. 31. l.h.s.: The potential energies per degree of freedom for gluons V gg /Nf = Tj^> g /N^ (dashed 
blue line) and for quarks (antiquarks) V qq /N q = T^ Q q /N q (solid red line) as a function of T/T c . All 
energies are multiplied by the dimensionless factor (T c /T). r.h.s.: The plasma parameter _T (124) for 
gluons (dot-dashed green line) and quarks (solid orange line) as a function of T/T c . Note that r ~ 1 — 2 
separates a gas phase from a liquid phase in case of Lenard- Jones type of interactions. 



Furthermore, the plasma parameter r defined by the ratio of the average potential energy 
per particle to the average kinetic energy per particle, 



r g = , r q = , (124) 



is displayed in the r.h.s. of Fig. 31 for gluons (dot-dashed green line) and quarks (solid orange 
line) as a function of T/T c . Here the kinetic energy densities are evaluated as [97] 

Tkin.g = Tr+ (w - y/F*), T kin>q = fr+(w - yfi"). (125) 

The present results clearly indicate that the plasma parameters r g ,r q are larger than unity 
for both quarks and gluons up to 10 T c (except for the vicinity of T c ) such that the system 
should be in a liquid phase provided that some attractive interaction between the constituents 
persists. Note that the present evaluation of the plasma parameter r is entirely carried out 
within the DQPM and no longer based on estimates for the potential energy as in Refs. [96,97]. 
The present results indicate that the sQGP should persist for a large range in temperature (or 
energy density) and thus also show up in nucleus-nucleus collisions at Large Hadron Collider 
(LHC) energies. Consequently, a partonic liquid is expected to be seen also at LHC energies 
and as a consequence the observed scaling of elliptic flow of practically all hadrons with the 
number of constituent quarks (as seen at RHIC) should persist also at LHC. 



3.2.2 Selfenergies and effective interactions of time-like quasiparticles 

Since in transport dynamical approaches there are no thermodynamical Lagrange parameters 
like the inverse temperature (3 — T^ 1 or the quark chemical potential which have to be 
introduced in thermodynamics in order to specify the average values of conserved quantities (or 
currents in the relativistic sense), derivatives of physical quantities with respect to the scalar 
densities iV® (or time-like densities AT+) (x = g,q,q) are considered in the following (cf. Refs. 
[121,141]). 

The independent potential energy densities V x :— T 00 x now can be expressed as functions of 
the scalar densities (or iV+) instead of the temperature T (and/or quark chemical potential 
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Fig. 32. The parton potential energy density V (126) - divided by the parton density p p - from the 
DQPM (l.h.s., blue line) as a function of the parton density p p (127). The functional form of V(p p ) is 
well reproduced by the expression (129) when divided by p v (green line) such that the lines cannot be 
separated by eye. The r.h.s. shows the gluon fraction a (127) from the DQPM (dark blue line) as a 
function of p v together with the fit (130) (green line). Again both lines cannot be distinguished by eye. 



fiq). For a determination of mean-field potentials for gluons and quarks (antiquarks) it is useful 
to consider the partonic potential energy density 

V ■■= T QQ . g + T - 0j9 + T - 0j9 - = V gg + Vqq + V qg (126) 

and to separate a pure gluonic interaction density V gg from a pure fermionic interaction density 
V qq as well as a gluon- fermion interaction density V qg . Correspondingly a parton density p p 
and gluon fraction a is defined via 

N+ 

Pp = N ^ + N ^ + N± , a = — c %r T . (127) 

Hp 9 q q ' N+ + N+ + N} y ' 

In the present DQPM (for T c = 0.185 GeV) the parton density p p (127) turns out to be a 
simple function of temperature, 

rj~i \ / rj~i \ 3.15 



* [tJ ^ {%) '' m ~ ' • ^> 

such that the average distance between the partons is given by d(T/T c ) — p p w (T c /T) 105 
[fm] w T c /T [fm]. These relations allow to convert temperatures scales to geometrical scales in 
a simple fashion. 

In Fig. 32 the parton potential energy density V (126) - divided by the parton density 
p p - is shown as a function of p p (l.h.s.). The functional dependence of V on p p can be well 
approximated by the expression (cf. l.h.s. of Fig. 32) 

V{p p ) m 0.975p^- 292 - 0.7lpp 21 [GeV/fm 3 ] , (129) 

where the numbers in front carry a dimension in order to match the units in GeV/fm 3 . The 
gluon fraction a (127) is shown on the r.h.s. of Fig. 32 and is well approximated by 

a(p p ) = 0.29 - 0.075pp°- 28 - 0.15 exp(-1.6p p ). (130) 

Adding half of the interaction density V qp to the gluon part and fermion part separately, 
we have T^ g = V gg + 0.5V qg and T m q + T m q = V qq + Q.5V qg such that V qq - V gg =: AV = 
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T 00 + T 00 g — Tq . The relative fraction of this quantity to the total potential energy density 
is evaluated as 

«(ft>) = 4?" = ^= ^ , m - ■ ( 131 ) 



Using the Ansatz: 
then gives 



V BB + V qg = (1 - 0V (132) 



V gg = 0.5(1 -£- k)V , V OT = 0.5(1 -£ + k)V , ^ M = £V, (133) 

with still unknown fraction £ for the interaction density V qg . 

In order to determine mean-field potentials U g (p p ) for gluons or U q (p p ) for quarks (in the 
rest frame of the system) one has to consider the derivatives 



<97V+ ' qKHpJ d(N+ + 



which by virtue of (132) can be computed as 

i d(i- K + Qv d Pp i d{i + K + j)v d Pp 
U 'M = 2 di P aN+> vM = ~ 2 9i p d{N + + N +)- (135) 

The fraction £ of the interaction density - in principle a function of p p but here taken to be a 
constant - now can be fixed in comparison to the gluon mean-field from Ref. [121] where the 
pure Yang-Mills sector has be investigated in the same way. This leads to £ « 0.3 and separates 
the total potential energy density V into 26% for the gluon-gluon interaction part, 30% 
for the quark-gluon interaction part (including the antiquarks) and w 44% for the fermionic 
interaction part. 




The corresponding results for U g (p p ) and U q (p p ) are displayed in the l.h.s. of Fig. 33 in 
terms of the solid red line and dashed blue line, respectively, and show distinct minima at p p w 
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2.2 fm 3 which corresponds to an average partonic distance of w 0.77 fm. The actual numerical 
results for the mean-fields can be fitted by the expressions, 

U g ( Pp ) w 70 e -ft=/°- 31 + 2.65 P ° p 21 + 0.45 p° p A [GeV] , (136) 

U q ( Pp ) « 32 e -P P /°- 31 + 1.1 pY 1 + 0.3 p° p A1 [GeV] , 

where p p is given in fm -3 and the actual numbers in front carry a dimension in order to match 
to the proper units of GcV for the mean-fields. The ratio U g /U q w 2.05 as can be seen in the 
r.h.s. of Fig. 33 for a very wide range of parton densities p p within 10% accuracy. 

Some comments on the Lorentz structure of the mean fields U g and U q appear appropriate. 
Note that by taking the derivatives with respect to the time-like densities one implicitly assumes 
that a 4- vector current is the physical source of the sclfcncrgies and that U g , U q are the 0'th 
components of vector fields Ug and U q (u = 0,1,2,3). The spatial components are assumed 
to vanish in the rest frame of the system and can be evaluated by a proper Lorentz boost 
to the frame of interest. This implies that the dynamical forces (as space-time derivatives 
of a Lorentz vector) are Lorentz tensors as in case of QED (or vector sclfenergies as in the 
nuclear physics context [142]). On the other hand one might consider the notion of purely 
scalar selfenergies where derivatives of the potential energy density with respect to the scalar 
density (e.g. dT^ Q x /dN*) define an effective mass M* [142]. In a selfconsistcnt framework then 
the quasiparticles masses (111) should be given by the derivatives with respect to the scalar 
densities. However, this relation is not fulfilled at all since a numerical evaluation of the scalar 
derivatives gives effective masses that are larger by more than an order of magnitude than the 
quasiparticle masses introduced in (111)! This result might have been anticipated since the 
effective forces for a gauge (vector) field theory should be dominated by Lorentz forces as in 
case of QED. 

Some information on the properties of the effective gluon-gluon, quark-gluon and quark- 
quark interaction may be extracted from the second derivatives of the potential energy density 
V, i.e. 

v ^-^^ ^-t-^ f^V (137) 
V «°W - 9N +2 ~ 2 dpi {dN+J ' {167 > 

d 2 V qq ^ l d\l-Z + K)V ( d Pp V 

u qq\Pp) ■ 



d(N++N^) 2 2 dpi \d(N++N±) 



d 2 V qg 



qg[Pp> ■ d(N+ + N+)dN+ ~ dpi [d(N+ + N+) ) \dN+J ■ 

The numerical results for the interactions (137) are displayed in Fig. 34 (l.h.s.) for the effective 
gluon-gluon (solid red line) , gluon-quark (solid blue line) and quark-quark interaction (dashed 
green line). All interactions show up to become strongly attractive at low parton density p p < 2.2 
fm~ 3 , change sign and become repulsive for all higher parton densities. Note that the change of 
quasiparticle momenta (apart from collisions) will be essentially driven by the (negative) space- 
derivatives —VUj(x) = —dUj(pp)/dp p Vpp(x) which implies that the partonic quasiparticles 
(at low parton density) will bind with decreasing density, i.e. form 'glueballs', mesons, baryons 
or antibaryons dynamically close to the phase boundary and repel each other for p p > 2.2 fm~ 3 . 
Note that color neutrality is imposed by color-current conservation and only acts as a boundary 
condition for the quantum numbers of the bound/resonant states in color space. 

The r.h.s. of Fig. 34 displays the ratios v gg /v qq (solid red line) and v qg /v qq (solid blue line) 
as a function of the parton density p p and demonstrates that v qg /v qq ~ 9/4, which is the ratio 
of the Casimir eigenvalues. 

A straight forward way to model the parton condensation or clustering to confined glueballs 
or hadrons dynamically (close to the phase transition) is to adopt screened Coulomb- like poten- 
tials v c (r, A) with the strength j d 3 r v c (r, A) fixed by the interactions v gg (p p ),v qg (p p ),v qq (p p ) 
from (137) and a screening length A from 1QCD studies. For the 'dilute parton regime' (p p < 
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Fig. 34. l.h.s. The effective gluon-gluon (solid red line), gluon-quark (solid blue line) and quark-quark 
interaction (dashed green line) from the DQPM for £ = 0.3 (see text). The r.h.s. displays the ratios 
Vgg/v qq (solid red line) and v qg /v qq (solid blue line). The straight light blue line is the ratio of the 
Casimir eigenvalues, i.e. C g /C q = 9/4. 



2.2 fm -3 ), where two-body interactions should dominate, one may solve a Schrodinger (Dirac 
or Klein-Gordon) equation for the bound and/or resonant states. This task is not addressed 
further since for the actual applications in the PHSD approach (cf. Section 4) the formation of 
glueballs is discarded and the formation of resonant hadronic states close to the phase boundary 
is described by density-dependent transition matrix elements (~ \v qq (p p )\ 2 ) between partons of 
'opposite' color with fixed flavor content. In this way the energy-momentum conservation, the 
flavor current conservation as well as 'color neutrality' are explicitly fulfilled in the (PHSD) 
transport calculations for interacting particles with spectral functions of finite width [143]. 



3.3 Finite quark chemical potential // f/ 

The extension of the DQPM to finite quark chemical potential fi q is more delicate since a 
guidance by 1QCD is presently very limited. In the simple quasiparticle model one may use the 
stationarity of the thermodynamic potential with respect to self-energies and (by employing 
Maxwell relations) derive a partial differential equation for the coupling g 2 (T, which may be 
solved with a suitable boundary condition for g 2 (T,p, q = 0) [130]. Once g 2 (T,fi q ) is known one 
can evaluate the changes in the quasiparticle masses with respect to T and fi q , i.e. dM 2 /dp q 
and dM 2 /dT (for x = g,q, q) and calculate the change in the 'bag pressure' AB (cf. Refs. [130, 
144] for details). However, such a strategy cannot be taken over directly since additionally the 
quasiparticle widths 7x(r, p q ) have to be known in the (T, p q ) plane in case of the DQPM. 

In hard-thermal-loop (HTL) approaches [145, 146] the damping of a hard quark (or gluon) 
does not depend on the quark chemical potential explicitly [147] and one might employ (113) 
also at finite /i. This, however, has to be considered with care since HTL approaches assume 
small couplings g 2 and should be applied at sufficiently high temperature, only. Present 1QCD 
calculations suggest that the ratio of pressure to energy density, P/e, is approximately indepen- 
dent of [i q as a function of the energy density e [148]. Accordingly, the functional dependence of 
the quasiparticle width 7 on fi q and T has to be modeled in line with 'lattice phenomenology' 
(see below). 



3.3.1 A scaling hypothesis 



Assuming three light flavors (q — u, d, s) and all chemical potentials to be equal — p<i — 
Us = fi) equations (110) and (111) demonstrate that the effective gluon and quark masses are 
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a function of 

T* 2 = T 2 + ^r. (138) 

7T 

Since the coupling (squared) (112) is a function of T/T c a straight forward extension of the 
DQPM to finite p is to consider the coupling as a function of T*/T c (p) with a ^-dependent 
critical temperature, 

T c (p) « T c (p = 0)(1 - ^2^2) « r c (0)(l - 0.05^^). (139) 

The coefficient in front of the /^-dependent part can be compared to 1QCD calculations at 
finite (but small) p which gives 0.07(3) [149] instead of 0.05 in (139). Consequently one has to 
expect an approximate scaling of the DQPM results if the partonic width is assumed to have 
the form (113), 

g 2 (T*/T c (p)) 2c 
7 ,(T, p) - N c T In g2(T « /Tc(p)) , (140) 

„ (T ,,) N 2 -lg 2 (T*/T c (p)) 2c 

2N C 4tt ln g 2 (T*/T c (p)) ' 

where g 2 (T/T c ) has been replaced by g 2 (T* /T c (p)). In fact, as will be demonstrated below, this 
choice leads to an approximate independence of the potential energies per degree of freedom 
as a function of p q . Nevertheless, the conjecture (140) should be explicitly controlled by 1QCD 
studies for Nf=3 at finite quark chemical potential. Unfortunately this task is presently out 
of reach and one has to live with the uncertainty in (140) which is assumed in the following 
investigations. 

Within the scaling hypothesis (139), (140) the results for the masses and widths in Section 
3.2 stay about the same as a function of T* /T c (p) when dividing by the temperature T. This 
also holds approximately when displaying the masses and widths as a function of the parton 
density p p for different chemical potentials p as demonstrated in Fig. 35. The latter quantities 
can well be fitted by the expressions 

M g ( Pp ) « 0.41 p° p 255 + 0.38 P - QJ [GeV] , (141) 

7s (Pp)« 0.235 p°- 245 - 0.14 p- 2 [GeV], 

2 4 
m g {p P ) w ^Mg(p p ), lq { Pp ) ss - j g (pp), 

with p p given in units of fm~ 3 . Note that according to the parametrization (141) the width 
might become negative for very small p p ; in actual transport applications it will be set to zero 
(cf. Section 4). 



3.3.2 Time-like and space-like quantities 

The more interesting question is how the energy density e (115) and the pressure P (from (114)) 
change with quark chemical potential p = p q in the DQPM. This information is provided in Fig. 
36 where the upper l.h.s. shows the energy density e (115) (scaled in terms of T c q = T c (p = 0) = 
0.185 GeV ) as a function of T*/T c (p). Here a scaling of the 'temperature' T* with T c (p) (139) 
is used since the phase boundary changes with the quark chemical potential p. The energy 
density e is seen to scale well with (T/T c o) 4 as a function of temperature for T*/T c (p) > 3, 
however, increases slightly with p close to the phase boundary where the scaling is violated 
on the level of 20%. This violation in the scaling (seen in the upper left part of the figure) is 
essentially due to an increase of the pressure P which is displayed in the lower left part of the 
figure as a function of T* /T c (p) for the same chemical potentials p from p = to 0.21 GeV in 
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Fig. 35. The gluon mass M and width j g (l.h.s.) as a function of the parton density p v for various 
chemical potentials \i q — /x from /j, = to /i = 0.21 GeV in steps of 0.03 GeV. The r.h.s. displays the 
mass m for quarks and width y q for the same quark chemical potentials. 



steps of 0.03 GeV. Note that a quark chemical potential of 0.21 GeV corresponds to a baryon 
chemical potential of fis = 3/i =0.63 GeV which is already substantial and the validity of (139) 
becomes questionable. 

Since the pressure P is obtained from an integration of the entropy density s over temper- 
ature (114) the increase in P with /i can directly be traced back to a corresponding increase 
in entropy density. The latter is dominated by the time-like quasiparticle contributions thus 
'counting' the effective degrees of freedom, 

Pp = K + N+ +q , N+ +q = N+ + N+ . (142) 

The upper r.h.s. of Fig. 36 shows p p versus T*/T c (/z) (multiplied by T 3 /T 3 ) for chemical 
potentials from /j, = to fj, = 0.21 GeV in steps of 0.03 GeV. Indeed the scaled parton density 
increases with (i; this increase is most pronounced for lower temperatures and becomes sub- 
stantial for /i = 0.21 GeV. One should recall, however, that a tri-critical endpoint (in the T, [i 
plane) is expected for fig = 3/j w 0.4 GeV [150] which corresponds to \i w 0.13 GeV. When 
restricting to the interval < \x < 0.13 GeV the explicit change in the parton density p p with 
fj, stays very moderate. This also holds for the energy density e and the pressure P. 

Quite remarkably the potential energy per time- like fermion T^~ q +q/^q+q changes very little 
with \x as can be seen from the lower right part of Fig. 36 where the latter quantity is displayed 
for the same chemical potentials as before (lower green lines) as a function of T*/T c (/x). This also 
holds for the potential energy per time-like gluon T ~^ g /N+ (upper magenta lines). Accordingly 
the potential energy per time-like degree of freedom is essentially a function of T* /T c (/x) alone. 

The time-like densities for the fermions Nq +q - (142) (upper l.h.s., dark blue lines) and the 
space-like quantities N~ +q - (upper r.h.s., light blue lines) are shown in Fig. 37 as a function of 
T* /T c (/j) for chemical potentials /i from [i = to /i = 0.21 GeV in steps of 0.03 GeV (scaled 
by T^/T 3 ). Both quantities increase with \i in a comparable fashion such that their ratio stays 
approximately constant for T* > 1.5 T c (/x). The space-like quantity N~ (lower l.h.s., dashed 
green lines) practically is independent from /x as well as the time-like gluon density N£ (lower 
l.h.s., red lines). Accordingly the gluon fraction 

N+ 

a ( 7 » = /v ++y ' (w3) 

which is displayed on the lower r.h.s. of Fig. 37 for the same set of quark chemical potentials as 
a function of T*/T c (/i), decreases with \i. It drops to zero below the phase boundary because 
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Fig. 36. The energy density e(T* /T c (p)) (115) for quark chemical potentials p from to 0.21 GeV 
in steps of 0.03 GeV as a function of the scaled temperature T*/T a (fi) for Nf —3 (upper l.h.s.). The 
pressure P(T* /T c (p)) from the thermodynamical relation (114) for quark chemical potentials p from 
to 0.21 GeV in steps of 0.03 GeV as a function of the scaled temperature T* /T c (p) for Nf —3 (lower 
l.h.s.). The parton density p p (for quark chemical potentials p from to 0.21 GeV) is shown in the 
upper r.h.s. The potential energy per time-like 'gluon' (upper magenta lines) and the potential energy 
per time-like 'quark+antiquark' (lower green lines) for the same quark chemical potentials are displayed 
in the lower r.h.s. Note that e and P are scaled by the dimensionless factor (T c o/T) 4 (with T c o = 0.185 
GeV) while the parton density p p is scaled by the factor (T c o/T) 3 . 



the difference between the gluon effective mass and the fermion effective mass becomes large 
below T c . 

We continue with the time-like and space-like energy densities for the fermions and gluons 
as a function of fj, and T* /T c (p) which are displayed on the l.h.s. of Fig. 38. The upper part 
shows that the space- like and time-like energy densities for the fermions increase with p roughly 
in a similar fashion as the fermion 'densities' such that their ratio is approximately independent 
on p,. The lower part of Fig. 38 (l.h.s.) indicates that the space-like energy density for gluons 
(green lines) as well as the time-like energy density for gluons (red lines) are approximately 
independent from p, within line width. When separating the time-like fermion energy density 
into contributions from quarks q and antiquarks q (r.h.s. of Fig. 38) we find an increase of 
with p, which is not fully compensated by a decrease of ^ with quark chemical potential. 
Similar dependences on p and T* are found for the space-like sector (lower part, r.h.s.). Thus 
when summing up all space-like and time-like energy densities from the fermions and the gluons 
a small net increase in the total energy density with p survives (see also upper left part of Fig. 
36). 

We note in passing that derivatives of the various energy densities w.r.t. the time-like gluon 
or fermion densities (as investigated in detail in Section 3.3) are approximately independent of 
p such that the effective potentials U g {p p ), U q (p p ) and U q (p p ) (135) stay practically the same. 
Since this result may be inferred already from the /x-(in)-dependence of the potential energy per 
time-like degree of freedom (displayed in the lower r.h.s. of Fig. 36) an explicit representation is 
discarded. This implies that the mean-fields (135) or parametrizations (136) may be employed 
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Fig. 37. Upper parts: The time-like fermion density N^ + „ (l.h.s., dark blue lines) and the space-like 
quantity N~ + = (r.h.s., light blue lines) for quark chemical potentials p, from to 0.21 GeV in steps of 
0.03 GeV as a function of the scaled temperature T* /T c (p). Lower parts: the time-like gluon density 
(l.h.s., lower red lines) and the space-like quantity Ng (l.h.s., dashed green lines) for the same 
quark chemical potentials as a function of T* /T c (p). Both, N~ and change only very little with 
p in the whole temperature range; the relative changes are of the order of the line width. The gluon 
fraction (143) is shown on the r.h.s. as a function of the scaled temperature T* /T c (p) for quark chemical 
potentials p from to 0.21 GeV in steps of 0.03 GeV. Note that all 'densities' have been multiplied by 
the dimensionless factor (T c o/T) 3 . 



also at finite (moderate) net quark density — which simplifies an implementation in 
parton transport models (as e.g. PHSD). 

Whereas for vanishing quark chemical potential p—0 the quark and antiquark densities are 
the same this no longer holds for finite p where the differences 

p± = N± - N± (144) 

are of separate interest since p q = /?+ + p~ is the zero'th component of a conserved flavor current 
(separately for each flavor u, d, s, ..). In order to obtain some idea about space-like and time-like 
contributions of the net quark density p q in the DQPM we first plot the time-like component 
p+ in the l.h.s. of Fig. 39 as a function of T* /T c (p) for different p from to 0.3 GeV (as 
before). In order to divide out the leading scaling with temperature the time-like densities have 
been multiplied here by T^ /T 2 on the r.h.s. of Fig. 39 (as known from Fermi systems in the 
nuclear physics context). The actual results show an approximately linear increase in p (l.h.s.) 
which suggests to study the scaled quantities • T c 2 /(/LtT 2 ). However, a rather good scaling is 
obtained when multiplying by T^ /(pT* 2 ). The latter quantities are displayed in the r.h.s. 
of Fig. 39 for the same set of chemical potentials as before. This approximate scaling allows to 
estimate the net quark density as 

Pq « 10 p T* 2 /T% [fm- 3 ] (145) 
with p given in units of GeV in case of Nf =3. 
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Fig. 38. l.h.s.: The time-like fermion energy density Tqq q+q - (blue lines) and the space-like fermion 
energy density Tj^ q+q (dashed red lines) for quark chemical potentials /j, from to 0.21 GeV in steps 
of 0.03 GeV as a function of the scaled temperature T* /T c (n) (upper l.h.s.). The lower l.h.s. shows 
the time-like (dashed red lines) and space-like gluon energy density (green lines) for the same chemical 
potentials as a function of the scaled temperature. Both, and are roughly independent from 

/i in the whole temperature range. The time-like (upper r.h.s) and space-like energy densities (lower 
r.h.s.) for quarks (upper red lines) and antiquarks (lower green lines) are separately displayed for quark 
chemical potentials /j, from to 0.21 GeV in steps of 0.03 GeV for Nf = 3. All quantities have been 
multiplied by the dimensionless factor T^/T 4 . 



3.3.3 Comparison to lattice QCD 



The approximate scaling depicted in Fig. 39 is a prediction of the DQPM in case of 3 quark 
flavors and should be controlled by 1QCD calculations. For two light flavors (Nf — 2) some 
comparison can be made in the low temperature (and low fj,) range. Unfortunately the simple 
scaling relations (145) do not hold at small T such that an explicit comparison has to be 
presented between the DQPM for Nj — 2 and the 1QCD calculations from Ref. [151]. The 
latter calculations have been carried out on a 16 3 x 4 lattice with two continuum flavors (of 
p4-improved staggered fermions) with mass m = 0.4T. Though the fermion masses are not 
really 'light' in the 1QCD calculations the actual 1QCD results may serve a test of the present 
DQPM in the 2-flavor sector (for low \i and T). The explicit 1QCD results for the net quark 
density p q (divided by T 3 ) are displayed in Fig. 40 as a function of T /T c $ in terms of the various 
symbols. The symbols of equal color correspond to ^/T c0 = 1.0, 0.8, 0.6, 0.4, and 0.2 from top 
to bottom and by eye show an approximately linear dependence on fj,. The explicit results from 
the DQPM for Nf = 2 are presented in terms of the dark green lines and approximately follow 
the 1QCD results (at least for smaller fj). Since the systematic errors of the 1QCD calculations 
are not known to the author an explicit refitting of the parameters A, T s /T c and c for Nf =2 is 
discarded here. Nevertheless, the qualitative (and partly quantitative) agreement between the 
DQPM and 1QCD provides a test for the basic concepts of the DQPM. 
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Fig. 39. l.h.s.: The time-like quark density (144) (scaled by the dimensionless factor T 2 /T 2 ) for 
quark chemical potentials p from to 0.21 GeV in steps of 0.03 GeV as a function of the scaled 
temperature T*/T c (fJ,). r.h.s.: The time-like quark density p+ (144) (red solid lines) and the space-like 
quantity p~ (144)(dashed green lines) for the same quark chemical potentials p. In this part of the 
figure all quantities are scaled by the factor T 2 /T* 2 / /i. 




Fig. 40. The quark density p q (144) (divided by the factor T 3 ) for quark chemical potentials fi/T c o 
from to 1.0 in steps of 0.2 as a function of the scaled temperature T/T c o from the DQPM (green 
lines) for Nf =2. The different symbols represent the 1QCD results from Ref. [151] for the same quark 
chemical potentials as a function of T/T c q (also for Nf =2). Note that the temperature axis here is 
given by T/T c0 = T/T c ((i = 0) and not by T*/T c (fj,)l 



3.4 Dilepton radiation from the sQGP 

The properties of the sQGP so far have been fixed in the DQPM by specifying the (vector) 
selfenergies/potentials as well as effective interactions for the time- like partons. As shown in Ref. 
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[97] this leads to a strongly interacting partonic system with a shear viscosity to entropy density 
ratio close to 77/s 0.2. However, the predictions from the DQPM should be controlled by 
independent 1QCD studies to get some idea about the reliability of the approach. As mentioned 
before transport coefficients like the shear viscosity 77 are available from 1QCD [109] but the 
present accuracy is not satisfactory. On the other hand some 1QCD information is available 
from the Bielefeld group on the electromagnetic correlator which is initimately related to the 
dilepton emission rate [140,152,153]. 
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Fig. 41. The 'back-to-back' dilepton emission rate (147) from the DQPM (solid red line) in comparison 
to the Born approximation (146) for massless partons (dashed blue line) and the results from the 1QCD 
analysis in Ref. [153] at T — 1.5 T c in case of Nf = 2 (full dots). 



In order to provide some further comparison to 1QCD results we consider the dilepton 
production rate in thermal equilibrium (at temperature T) in a two-flavor QGP as in Ref. [153]. 
For massless quarks and antiquarks the emission rate of 'back-to-back' leptons (qq — > e + e~) is 
given by [153,154] 

dW , 5a 2 lu uj 

^ (P = 0)= 36^^ ( 2T ) ^ ( 2T ) ' (W6) 

where u> is the invariant mass of the lepton pair and uf denotes the Fermi distribution function. 
In (146) a is the electromagnetic coupling constant. For massless partons, furthermore, the 
magnitude of the lepton momenta is given by |p| — lu/2 while their direction is opposite in the 
dilepton rest frame. Neglecting the rest mass of leptons their energy is uj/2 in the dilepton rest 
frame, too. The expression (146) changes in case of spectral functions with finite width to 

dw ha 2 r , r , f°° , uiw 2 



f f f LU UJ 

dudtp 30n* I J ^ I dp H P^P;T) Pq (^p;T) (147) 

with P? = uj 2 - — p 2 denoting the invariant mass of the annihilating partons j — 1,2. In (147) 
the factor vX/uj 2 gives a flux correction in case of massive quasiparticles with 

X(x,y,z) = (x-y- z) 2 -4yz. (148) 
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Since the spectral functions in the DQPM are fixed the lepton emission rate (147) can 
be evaluated without introducing any further assumption (or parameter). The results for the 
differential emission rate for Nf = 2 are shown in Fig. 41 for T = 1.5 T c (solid red line) in 
comparison to the limit (146) (dashed blue line) and the 1QCD results from [153] (full dots). 
The 1QCD dilepton rate has been obtained from the temporal correlators in 1QCD employing 
the 'maximum entropy method' which has a systematic error in the order of 30 to 50 % [153] 
depending on the energy scale considered. The results from Fig. 41 demonstrate a drastic 
suppression of low mass lepton pairs due to the finite mass of the partons. On the other hand 
the spectra from the DQPM are in qualitative agreement with the 1QCD results from Ref. [153] 
when including the systematic uncertainties of the latter approach. Thus the DQPM passes a 
further test in comparison to 1QCD. 

Some note of caution with respect to the present DQPM appears appropriate: the parameters 
in the effective coupling (112) and the width (113) have been fixed in the DQPM by the 
entropy (108) to 1QCD results for Nf=0 assuming the form (109) for the spectral function 
p(uj). Alternative assumptions for p(iv) will lead to slightly different results for the time-like and 
space-like densities, energy densities etc. but not to a qualitatively different picture. Also it is 
presently unclear if the three parameters (A, T s /T c , c) employed in Section 3.2 are approximately 
the same in case of two and three dynamical flavors. More precise calculations from 1QCD should 
allow to put further constraints on the form of the spectral function p(uj) and to fix the basic 
model parameters in the effective coupling. Also the 'modified' HTL expression (140) for the 
quasiparticle width has to be controlled by 1QCD at finite p q as well as transport coefficients 
like the shear viscosity r] or related correlators. 
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4 Dilepton production, parton propagation and hadronization 

The dynamics of strongly interacting systems are reflected in the spectral functions, i.e., in 
the imaginary part of the retarded/advanced Green functions. The question comes about how 
to measure spectral functions experimentally? A promising way to obtain information about 
the spectral properties of hadrons is to measure their electromagnetic decay to e+e~ or p + p~ 
pairs. This is particularly relevant for short lived particles in a hot and/or dense hadronic 
environment since this allows to look at 'in-medium' properties of time-like quanta because the 
decay products suffer only from a very small (electromagnetic) final-state- interaction. In fact, 
a variety of experimental collaborations have been set up to investigate this particular question 
(cf. the lecture by J. Stroth in this volume [71]). 



4.1 Dilepton production 

The theory of quantum-chromo-dynamics (QCD) describes hadrons as many-body bound or 
resonant states of partonic constituents. While the properties of hadrons are rather well known 
in free space (embedded in a nonperturbative QCD vacuum) the masses and lifetimes of hadrons 
in a baryonic and/or mesonic environment are subject of current research in order to achieve a 
better understanding of the strong interaction and the nature of confinement. In this context 
the modification of hadron properties (or spectral functions) in nuclear matter is of fundamental 
interest. Related research has focussed in the past on properties in the medium [86,155] as 
well as on vector mesons (p, u> and <p). QCD sum rules [156] as well as QCD inspired effective 
Lagrangian models [157] have predicted significant changes, e.g., of the p, lo and <j> with the 
nuclear density pm and/or tcmpcrture T. 

A direct evidence for the modification of the spectral properties of vector mesons has been 
obtained from the enhanced production of lepton pairs above known sources in nucleus- nucleus 
collisions at SPS energies [158, 159]. As proposed by Li, Ko, and Brown [160] and Ko et al. [161] 
the observed enhancement in the invariant mass range 0.3 < M < 0.7 GeV might be due to a 
shift of the p-meson mass following Brown/Rho scaling [157] or the Hatsuda and Lee sum rule 
prediction [156]. On the other hand also more conventional approaches that describe a melting 
of the p-mcson in the medium due to the strong hadronic coupling (along the lines of Refs. [162]) 
have been found to be compatible with the early CERES data [164]. This ambiguous situation 
has been clarified to some extent in 2005 by the NA60 Collaboration since the invariant mass 
spectra for p + pT pairs from In+In collisions at 158 A-GeV clearly favored the 'melting p' 
scenario [163]. Also the more recent data from the CERES Collaboration (with enhanced mass 
resolution) [165] show a preference for the 'melting p' scenario. 

Dileptons have also been measured in heavy-ion collisions at the BEVALAC by the DLS Col- 
laboration [166, 167] at incident energies that are two orders-of-magnitude lower than that at the 
SPS. The first published spectra at 1 A-GcV [166] (based on a limited data set) have been con- 
sistent with the results from transport model calculations [168] that include pn bremsstrahlung, 
7T°, r\ and A Dalitz decay and pion-pion annihilation. However, in 1997 the DLS Collaboration 
released a new set of data [167] based on the full data sample and an improved analysis, which 
showed a considerable increase in the dilepton yield: more than a factor of five above the early 
DLS data [166] and the corresponding theoretical results [168]. With an in- medium p spectral 
function, as that used in Ref. [164] for dilepton production from heavy-ion collisions at SPS 
energies, dileptons from the decay of both directly produced p's and pion-pion annihilation have 
been considered, and a factor of about two enhancement has been obtained in the theoretical 
studies compared to the case of a free p-spectral function. Furthermore, in Ref. [169] the alter- 
native scenario of a dropping p-meson mass and its influence on the properties of the iV(1520) 
resonance has been investigated. Indeed, an incorporation of such medium effects lead to an 
enhancement of the p-meson yield, however, was not sufficient to explain the DLS data. Since 
independent transport calculations by HSD (BUU) and UrQMD underestimated the DLS data 
for C+C and Ca+Ca at 1 A GeV by roughly the same amount these findings have lead to the 
denotation 'DLS-puzzlc' in 1999 [169] which persisted by about a decade. 
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In order to address the in- medium modifications of vector- mesons from the experimental side 
the HADES spectrometer has been built [170,71] which allows to study e + e~ pair production 
in a much wider acceptance region for elementary pp, pd reactions as well as it A, pA or even 
AA collisions up to about 8 A-GeV. Meanwhile the HADES Collaboration has presented first 
spectra [170 172] and the question is: what do these data tell us? The answer to the questions 
raised above is nontrivial due to the nonequilibrium nature of these reactions, and transport 
models have to be incorporated to disentangle the various sources that contribute to the final 
dilcpton spectra seen experimentally. This task requires an off-equilibrium transport approach 
which includes the dynamical evolution of spectral functions. In the following we briefly present 
results from an up-to-date relativistic transport model (HSD) that incorporates the relevant 
off-shell dynamics of the vector mesons (in line with Section 2) [173]. 

In the off-shell HSD approach the dilepton production by a (baryonic or mesonic) resonance 
R decay can be schematically presented in the following way: 



BB -> RX (149) 

mB -> RX (150) 

R-fe+e-X, (151) 

R^mX, m^e+e-X, (152) 

R -» R'X, R' -» e+erX, (153) 



i.e. in a first step a resonance R might be produced in baryon-baryon (BB) or meson-baryon 
(mB) collisions (149), (150). Then this resonance can couple to dileptons directly (151) (e.g., 
Dalitz decay of the A resonance: A — > e + e~N) or decays to a meson m (+ baryon) or in 
(152) produce dileptons via direct decays (p,uj) or Dalitz decays (n°,T),u). The resonance R 
might also decay into another resonance R' (153) which later produces dileptons via Dalitz 
decay. Note, that in the combined model the final particles - which couple to dileptons - can be 
produced also via non-resonant mechanisms, i.e. 'background' channels at low and intermediate 
energies or string decay at high energies [87]. 

The electromagnetic part of all conventional dilepton sources - 7r°, 77, u> Dalitz decays, direct 
decay of vector mesons p, u> and <fi - are treated as described in detail in Ref. [79] . Modifications 
- relative to Ref. [79] - are related to the Dalitz decay of baryonic resonances and especially 
the strength of the pp and pn bremsstrahlung since more recent calculations by Kaptari and 
Kampfcr [174] indicated that the latter channels might have been severely underestimated 
in previous studies on dilcpton production at SIS energies. In detail: For the bremsstrahlung 
channels in pp and pn reactions we adjust the previous expressions in order to match the recent 
results from the OBE model calculations in Ref. [174]. 



4.1.1 Vector-meson spectral functions 

In order to explore the influence of in-medium effects on the vector-meson spectral functions 
we incorporate the effect of collisional broadening, i.e. the vector meson width has been imple- 
mented as: 

r^(M, | P |, Pn ) = r v (M) + r coll (M, \ P \, PN ). (154) 

Here 7V(M) is the total width of the vector mesons (V — p, lo) in the vacuum. For the p meson 
we use 

n f \ 2 / ^ 3 
M \ q 



r p (M)~r p ^(M) = r Q ^) (j-j f(m) (155) 

(M 2 - Am 2 ,) 1 / 2 (M§ - Ami) 1 / 2 
9 = 7, 1 9o = 7. ■ 



Will be inserted by the editor 



71 



In (155) M is the vacuum pole mass of the vector meson spectral function, F(M) is a formfactor 
taken as 



2A 2 + M 2 

with a cut-off parameter A = 3.1 GeV. For the w meson a constant total vacuum width is used: 
r u = r u (Mo), since the a; is a narrow resonance in vacuum. The collisional width in (154) is 
approximated as 

r coll (M, \p\, PN ) = 7 p N (v a%) « a coll ^. (157) 

Po 

Here v = \p\/E; p, E are the velocity, 3-momentum and energy of the vector meson in the 
rest frame of the nucleon current, 7 is the Lorentz factor for the boost from the nucleon rest 
frame to the center-mass-system of heavy-ions; pn denotes the nuclear density and ay l N is the 
meson-nucleon total cross section. 

In order to simplify the actual calculations for dilepton production the coefficient a co u has 
been extracted in the HSD transport calculations from the vector-meson collision rate in C+C 
and Ca+Ca reactions (at 1 and 2 A- GeV) as a function of the density pjy. The numerical 
results for r co u(pif) then have been divided by pn/po to fix the coefficient a co u in (157). We 
obtain a co u « 150 MeV for the p and a co u w 70 MeV for ui mesons. In this way the average 
effects of collisional broadening are incorporated and allow for an explicit representation of the 
vector-meson spectral functions versus the nuclear density (see below). 

In order to explore the observable consequences of vector-meson mass shifts at finite nuclear 
density - as indicated by the CBELSA-TAPS data [175] for the u meson - the in-medium 
vector-meson pole masses are modeled (optionally) according to the Hatsuda and Lee [156] or 
Brown/Rho scaling [157] as 

M*(p N ) = -—^ (158) 
(1 + ap N /p ) 

where pn is the nuclear density at the resonance decay position r, p = 0.16 fm~ 3 is the normal 
nuclear density and a ~ 0.16 for the p and a ~ 0.12 for the u> meson. The parametrization 
(158) may be employed also at much higher collision energies (e.g. FAIR and SPS) and one 
does not have to introduce a cut-off density in order to avoid negative pole masses. Note that 
(158) is uniquely fixed by the 'customary' expression Mq(pn) ~ M (l — ap^/po) in the low 
density regime. 

The spectral function of the vector meson V with the mass M at baryon density pn is taken 
in the Breit-Wigner form: 

with the normalization condition for any p^: 

rM, im 

/ A v {M,p N )dM = 1, (160) 

JM mi „ 

where Mi im = 2 GeV is chosen as an upper limit for the numerical integration. The lower limit 
of the vacuum spectral function corresponds to the 27T decay M min = 2m x , whereas for the 
in-medium collisional broadening case M min — 2m e — > with m e denoting the electron mass. 
Mp is the pole mass of the vector meson spectral function which is Mq(pn = 0) = M in 
vacuum, however, shifted in the medium for the dropping mass scenario according to Eq. (158). 

The resulting spectral functions for the p and to meson are displayed in Fig. 42 for the case of 
'collisional broadening' (upper part) as well as for the 'dropping mass + collisional broadening' 
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Fig. 42. The spectral functions for the p and u) meson in the case of the 'collisional broadening' scenario 
(upper part) and the 'dropping mass + collisional broadening' scenario (lower part) for nuclear densities 
of 0,l,2,3,5xpo as employed in the transport calculations (see text for details). 



scenario (lower part) for densities of 0,1,2,3,5 xpo- Note that in vacuum the hadronic widths 
vanish for the p below the two-pion mass and for the w below the three-pion mass. With increas- 
ing nuclear density p^ elastic and inleastic interactions of the vector mesons shift strength to 
low invariant masses. In the 'collisional broadening' scenario we find a dominant enhancement of 
strength below the pole mass for the p-meson while the w-meson spectral function is drastically 
enhanced in the low- and high- mass region with density (on expense of the pole-mass regime). 
In the 'dropping mass + collisional broadening' scenario both vector mesons dominantly show 
a shift of strength to low invariant masses with increasing p^. Qualitatively similar pictures 
are obtained for the </>-meson but quantitatively smaller effects are seen due to the lower effect 
of mass shifts and a substantially reduced (j>N cross section which is a consequence of the ss 
substructure of the </)-meson. 

In order to demonstrate the importance of off-shell transport dynamics for dilepton produc- 
tion we present in Fig. 43 the time evolution of the mass distribution of p (upper part) and u> 
(lower part) mesons for central C+C collisions (b=l fm) at 2 A GeV for the dropping mass + 
collisional broadening scenario (as an example). The l.h.s. of Fig. 43 corresponds to the calcu- 
lations with on-shell propagation whereas the r.h.s. stand for the off-shell dynamics. As seen 
in Fig. 43 the initial p and u> mass distributions are quite broad even for a small system such 
as C+C where, however, the baryon density at 2 A GeV can reach (in some local cells) ~ 2po- 
The number of vector mesons decreases with time due to their decays and the absorption by 
baryons (pn — > ttN or 2irN). Most of the p mesons decay /disappear already inside the 'fireball' 
for density pn > 0. Due to the expansion of the 'fireball' the baryon density drops quite fast, 
so some amount of p mesons reach the very low density zone or even the 'vacuum'. Since for 
the off-shell case (r.h.s. of Fig. 43) the p spectral function changes dynamically by propagation 
in the dense medium, it regains the vacuum shape for p^r — > 0. This does not happen for the 
on-shell treatment (l.h.s. of Fig. 43) - the p spectral function does not change its shape by 
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Fig. 43. Time evolution of the mass distribution of p (upper part) and uj (lower part) mesons for 
central C+C collisions (b=l fm) at 2 A GeV for the dropping mass + collisional broadening scenario. 
The l.h.s. of Fig. 43 corresponds to the calculations with on-shell dynamics, whereas the r.h.s. shows 
the off-shell HSD results. 



propagation but only by explicit collisions with other particles. Indeed, there is a number of p's 
which survive the decay or absorption and leave the 'fireball' with masses below 2711^. Thus, 
the on-shell treatment leads to the appearance of the p mesons in the vacuum with M < 2™^, 
which can not decay to 2ir; thus they live practically 'forever' since the probability to decay to 
other channels is very small. Indeed, such p's will continuously shine low mass dileptons which 
leads to an unphysical 'enhancement/divergence' of the dilepton yield at low masses (note, that 
the dilepton yield is additionally enhanced by a factor 1/M 3 ). 

The same statements are valid for the uj mesons (cf. the lower part of Fig. 43): since the uj 
is a long living resonance, a larger amount of a/s survives with an in-medium spectral function 
shape in the vacuum (in case of on-shell dynamics). Such w's with M < Sm^ can decay only to 
7T7 or electromagnetically (if M < m n ). Since such unphysical phenomena appear in on-shcll 
transport descriptions including an explicit vector-meson propagation an off-shell treatment is 
mandatory! 



4.1.2 Dilepton production in C+C and Ca+Ca collisions at SIS energies 

A major aim of the HADES Collaboration was to clarify the 'DLS-puzzle', i.e. to verify /falsify 
the DLS data [167] from the experimental side [71]. Accordingly, the same systems have been 



74 



Will be inserted by the editor 



reinvestigated at 1 A- GeV with the HADES detector in order to clarify the issue. The situation, 
however, is not as easy since the DLS and HADES acceptances differ significantly [f 67, 170]. 
Recently the HADES Collaboration has performed a direct comparison with the DLS data for 
C+C at 1 A GeV [167] by filtering the HADES data with the DLS acceptance [172]. Both 
measurements were found to agree very well [71,172]! Thus, the 'DLS puzzle' has been solved 
from the experimental side; there is no obvious contradiction between the DLS and HADES 
data! 

The major goal for the transport models now was to verify/explain the solution of the 'DLS- 
puzzle' from the theoretical side. Note that transport calculations allow for a direct comparison 
between the DLS and HADES measurements by employing the different FILTER routines to 
the same set of calculated events. 



Comparison to the DLS data 



We start with a reinvestigation of the DLS data employing the additional (and enhanced) 
bremsstrahlung channels as mentioned above [173]. The results of the HSD transport calculation 
are displayed in Fig. 44 for 12 C+ 12 C (l.h.s.) at 1.04 A-GeV in case of 'free' vector-meson spectral 
functions (upper part) and in case of the 'dropping mass + collisional broadening' scenario 
(lower part) employing the DLS acceptance filter and mass resolution. 




0.3 0.4 0.5 

M [GeV/c 

Fig. 44. Results of the HSD transport calculation for the mass differential dilepton spectra in case 
of 12 C + 12 C (l.h.s.) and Ca + Ca (r.h.s.) at 1.04 A- GeV in comparison to the DLS data [167]. The 
upper parts show the case of 'free' vector-meson spectral functions while the lower parts give the result 
for the 'dropping mass + collisional broadening' scenario. In both scenarios the DLS acceptance filter 
and mass resolution have been incorporated. The different color lines display individual channels in the 
transport calculation (see legend). 
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Fig. 45. Results of the HSD transport calculation for the mass differential dilepton spectra - divided 
by the average number of 7r°'s - in case of 12 C + 12 C at 1.0 A-GeV in comparison to the HADES 
data [172]. The upper part shows the case of 'free' vector-meson spectral functions while the lower 
part gives the result for the 'dropping mass + collisional broadening' scenario. In both scenarios the 
HADES acceptance filter and mass resolution have been incorporated. The different color lines display 
individual channels in the transport calculation (see legend). 



In fact, the situation has substantially improved compared to the early studies, and the 
missing yield in the 'free' scenario is now reduced to a factor of about 1.5 in the mass region 
from 0.25 to 0.5 GeV. This is due to slightly higher contributions from A and r\ Dalitz de- 
cays and a significantly larger yield from bremsstrahlung channels. As seen from Fig. 44, the 
bremsstrahlung yield now is similar in shape and magnitude as the Z\-Dalitz decay contribu- 
tion and enhanced by factor of up to 5 due to the novel NN bremsstrahlung cross section from 
Kaptari [174] and accounting for the contribution from pp bremsstrahlung additionally to pn. 
The contribution of the pion-nucleon bremsstrahlung is quite small in the C+C system due to 
the limited energy available in meson-baryon collisions and the moderate rescattering rate of 
pions. 

Medium modifications for the vector mesons turn out to yield an enhancement in the region 
of 0.4 < M < 0.5 GeV but are very moderate due to the light system C+C. Though a slightly 
better description of the DLS data is achieved it would be premature to claim the presence of 
in-medium effects on the vector mesons from these data. 

Additional information is provided by the 40 Ca+ 40 Ca system at 1.04 A-GeV where the 
measured DLS spectra are reproduced in a qualitatively and quantitatively similar manner by 
the transport calculations (cf. r.h.s. Fig. 44) as in case of the lighter system. Again the combined 
bremsstrahlung channels provide a contribution in the same order as the A Dalitz decay; the 
r\ Dalitz decay remains subdominant but the ttN bremsstrahlung increases compared to the 
C + C system due to more frequent pion rescattering on nucleons. In- medium effects for the 
vector mesons can be indentified in the transport calculations but are hard to see in the total 
dilepton spectra. 
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Is the 'DLS puzzle' solved? This question can only be answered in a convincing manner by 
comparison with the recent HADES data. 

Comparison to HADES data 

The same (HSD) multi-differential dilepton spectra - used for comparison with the DLS data 
- are now filtered by the HADES acceptance routines and smeared with the HADES mass 
resolution [170]. A comparison of the HSD calculations for the 12 C + 12 C system at 1.0 A-GeV 
with the HADES data from [172] is presented in Fig. 45 and demonstrates that the agreement 
between data and transport calculations is reasonable for the HADES data, too. In this case 
the differential dilepton spectrum is divided by the average number of 7r°'s which in experiment 
is determined by half of the average number of (7r + + 7r~). The spectra look slightly different 
in shape due to the much higher acceptance at low invariant mass where the 7r° Dalitz decay 
practically exhausts the dilepton spectrum. As in case of the DLS data the r\ Dalitz decay turns 
out to be subdominant and the A Dalitz decay to be of similar magnitude as the combined 
bremsstrahlung contribution. Again the effect of in-medium vector-meson spectral functions is 
hard to see in the total spectra. 

In summary: The 'DLS puzzle' appears to be solved from the theoretical side, too, if en- 
hanced Bremsstrahlung channels are accounted for [173]. Note, however, that the problem of 
radiative corrections especially in inelastic collisions of charged hadrons is presently barely 
understood from the theoretical side! 



4.1.3 Dilepton production at SPS energies 

Let's continue with the actual results from HSD for In+In collisions at 160 A-GeV [176] which 
are displayed in Fig. 46 for peripheral, semi-peripheral, semi-central and central collisions fol- 
lowing the centrality definition of the NA60 Collaboration [163]. The dashed (blue) lines give the 
dilepton mass spectra when incorporating only the vacuum p- spectral function in all hadronic 
reaction processes. These reference spectra overestimate the data in the region of the p-meson 
pole mass and underestimate the experimental spectra in the region below (and above) the 
pole mass such that in-medium modifications can clearly be identified. The solid red lines show 
the result from HSD in the 'collisional broadening' scenario where no shift of the p pole mass 
is incorporated but an increase of the p-meson width due to hadronic collisions proportional 
to the baryon density pb (calculated as p%(x) = j fJ '(x)j ll (x) with j^(x) denoting the baryon 
four-current). As can be seen from Fig. 46 the dilepton mass spectrum is rather well described 
up to invariant masses of 0.9 GeV. For higher invariant masses the data signal additional con- 
tributions. This result is practically identical to the calculations of van Hees and Rapp [177] in 
the expanding fireball model - when incorporating the spectral function from Ref. [178] - and 
demonstrates that the dominant in-medium effect seen in the p, + p~ spectra from NA60 is a 
broadening of the p meson. 

In Ref. [179] a very detailed analysis of the NA60 data has been performed and shown that 
this additional yield above 0.9 GeV partly is due to open charm decays, four pion collisions 
or 'quark-antiquark' annihilation. Wc mention that our HSD calculations give only a small 
contribution from tt + a\ collisions in this invariant mass range but a preliminary study within 
the Parton-Hadron-String-Dynamics (PHSD) model [180] suggests that - apart from open charm 
decays - the extra yield seen experimentally by NA60 should be due to massive 'quark-antiquark' 
annihilations. Since this question is presently open and discussed controversely we concentrate 
on low mass dilepton pairs in the following. 

The next step in our study is related to an update of the HSD calculations in comparison 
to the recent data from the CERES Collaboration [181] (with enhanced mass resolution). In 
Fig. 47 we present the HSD results for the mass-differential dilepton spectra in 30% central Pb 
+ Au collisions at 40 A GeV (upper part) and 7% central 158 A GeV (lower part) in comparison 
to the data from CERES [181]. The dashed lines show the results in case of vacuum spectral 
functions (for p, uj, <f>) whereas the solid lines correspond to the 'collisional broadening' scenario. 
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Fig. 46. The HSD results for the mass differential dilepton spectra in case of In + In at 158 A-GeV for 
peripheral, semi-peripheral, semi-central and central collisions in comparison to the data from NA60 
[163]. The actual NA60 acceptance filter and mass resolution have been incorporated. The solid red 
lines show the HSD results for a scenario including the collisional bradening of the p-meson whereas 
the dashed blue lines correspond to calculations with 'free' p spectral functions for reference. 

Similar to the In+In case the experimental data agree better with the calculations employing 
the 'collisional broadening' scenario. This is also in line with Ref. [179]. On the other hand, the 
HSD model underpredicts the yield between the u> and cj> peaks which might again be attributed 
to possible contributions from 'quark-antiquark' annihilations. 

4.1.4 Dilepton production at RHIC energies 

We step on to RHIC energies and first compare the HSD results [176] for the dilepton invariant- 
mass spectrum from pp collisions at y^s = 200 GeV with the data from PHENIX (incorporating 
the acceptance cuts and mass resolution from PHENIX) [182] in Fig. 48. Actually the calcu- 
lations well reproduce the experimental spectrum which can entirely be described by meson 
Dalitz and direct decays as well as some contribution from open charm decays (light blue solid 
line as calculated by PYTHIA). This comparison demonstrates that the hadron production 
channels in HSD for elementary pp collisions are well under control. 

We recall that HSD also provides a reasonable description of hadron production in Au+Au 
collisions at yfs = 200 GeV [183] such that we can directly continue with the results for e + e~ 
pairs which are shown in Fig. 49 in case of inclusive Au + Au collisions in comparison to the 
data from PHENIX [184]. In the upper part of Fig. 49 the results are shown for vacuum spectral 
functions (for p,uj,<fi) including the channel decompositions (see legend for the different color 
coding of individual channels). Whereas the total yield (upper blue line) is quite well described in 
the region of the pion Dalitz decay as well as the u>- and 0-mass regime we clearly underestimate 
the measured spectra in the regime from 0.2 to 0.6 GeV by an average factor of 3. 

When including the 'collisional broadening' scenario for the vector mesons we achieve the 
sum spectrum shown by the green line in the lower part of Fig. 49 which is only slightly 
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Fig. 47. The HSD results for the mass differential dilepton spectra in 30% central Pb + Au collisions 
at 40 A-GeV (upper part) and 7% central 158 A GeV (lower part) in comparison to the data from 
CERES [181]. The dashed lines show the results for vacuum spectral functions (for p, u), <f)) whereas the 
solid lines correspond to the 'collisional broadening' scenario. 



enhanced compared to the 'free' scenario (blue line). Thus the question emerges if the PHENIX 
data might signal dropping vector meson masses contrary to the NA60 data (in Fig. 46)? 
To answer this question we have performed additional calculations in the 'dropping mass + 
collisional broadening' model where the p and u) masses have been dropped with baryon density 
in accordance with Eq. (158). The respective HSD results are displayed in the lower part of 
Fig. 49 by the red line and indeed show a further enhancement of the dilepton yield which, 
however, is only small in the mass range 0.2 GeV < M < 0.4 GeV such that also this possibility 
has to be excluded in comparison to the PHENIX data. 

In summary the presently available dilepton data from In+In collisions at 158 A-GeV (from 
the NA60 Collaboration) as well as the low mass dilepton spectra for Pb+Au collisions at 40 and 
158 A-GeV (from the CERES Collaboration) are well described in the 'collisional broadening 
scenario'. However, the low mass dilepton spectra from Au+Au collisions at RHIC (from the 
PHENIX Collaboration) are clearly underestimated in the invariant mass range from 0.2 to 0.6 
GeV. This also holds for the 'dropping mass + collisional broadening' scenario, i.e., when ten- 
tatively assuming a shift of the vector meson mass poles with the baryon density. We attribute 
this additional low mass enhancement seen by PHENIX to non-hadronic sources, possibly to 
virtual gluon-Compton scattering. 



4.2 Off-shell parton dynamics 

The Parton-Hadron-String Dynamics (PHSD) approach is an extension of the more familiar 
HSD model [79,91] with respect to the off-shell dynamics of partons (gluons, quarks, anti- 
quarks) . This extension is demanded by experiment since the HSD approach - without explicit 
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partonic degrees of freedom - failed in a couple of aspects to describe the rich phenomena ob- 
served experimentally at Relativistic-Heavy-Ion-Collider (RHIC) energies of y/s — 200 GeV (cf. 
Refs. [101-103]). 

The off-shell dynamics of the partons is entirely described by the generalized transport 
equations in BM form within the testparticle representation for iG < (Section 2) incorporating 
the effective selfenergies from the DQPM (Section 3). In PHSD the following elastic and inelastic 
interactions are included qq <-> qq, qq <-> qq, gg <-> gg, gg <-> g, qq <-> g etc. on the partonic 
side exploiting 'detailed-balance' with cross sections extracted in Refs. [97,121]. In this way 
the partonic evolution is fixed, however, the effective parton quasi-particlcs are ill defined as 
asymptotic states in vacuum since the 'masses' diverge for T/T c — 0.46 (cf. equations (3.3) 
- (3.5)). Consequently the partons (in PHSD) live only at sufficiently high temperature (or 
energy density) and have to recombine to color neutral hadrons in the expansion of the system. 

In PHSD the transition from partonic to hadronic degrees of freedom is described by local 
covariant transition rates e.g. for q + q fusion to a meson m of four-momentum p = [u>, p) at 
space time point x = (i,x) [143]: 



dN m (x,p) 4 . . 4 / x q + x q 

- Tr q TTq S (p - p q - pq) 6 



q -x 



d*xd 4 p v r« W y 2 

X ^g PqiPq) Uq Pq(Pq) \v q q\ 2 W{x q - Xq,p q ~Pq) 

xNq(x qi p q ) Nq(xq,pq) <5(flavor, color). (161) 



In (161) we have introduced the shorthand notation 

^ (277 



^ = E / ^ (162) 



where ■ denotes a summation over discrete quantum numbers (spin, flavor, color); Nj(x,p) 
is the phase-space density of parton j at space-time position x and four-momentum p. 5( 
flavor, color) stands symbolically for the conservation of flavor quantum numbers as well as 
color neutrality of the formed hadron m. Furthermore, v q q(p p ) is the effective quark- antiquark 
interaction as extracted from the DQPM in Section 3.3 as a function of the local parton (q+q+g) 
density p p (or energy density e) while W(x,p) is the phase-space distribution of the formed 
hadron. It is taken as a Gaussian in coordinate and momentum space with width \J (r 2 ) = 
0.66 fm for a meson and 1 fin for a baryon/antibaryon (in the rest frame). The width in 
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Fig. 49. The HSD results for the mass differential dilepton spectra in case of inclusive Au + Au 
collisions at y^s = 200 GeV in comparison to the data from PHENIX [184]. In the upper part the 
results are shown for vacuum spectral functions (for p,U),<ft) including the channel decompositions (see 
legend for the different color coding of the individual channels). The lower part shows a comparison 
for the total e + e~ mass spectrum in case of the 'free' scenario (blue line), the 'collisional broadening' 
picture (green line) as well as the 'dropping mass + collisional broadening' model (red line). 



momentum space then is fixed by the minimal uncertainty principle. We note that the final 
hadron formation rates are approximately independent on these parameters (within reasonable 
variations). Related transition rates are defined for the fusion of three off-shell quarks (q+q+q «-> 
B) to color neutral baryonic (B) resonances (or B) of finite width (or strings) with the help 
of Jacobi coordinates. The transition rates (161) fulfill energy and momentum conservation as 
well as flavor current conservation since the parton flavors fix the flavor content of the formed 
meson m or baryon B. 

According to the DQPM the effective interaction turns to be strongly attractive below 
p c « 2.2 fm -3 (cf. Section 3.3) which implies a dynamical binding of partons to composite 
(color neutral) hadrons once the parton density drops below p c . Thus the hadronization starts 
in PHSD when the local parton density falls below 2.2 fm -3 , i.e., the transition rate (161) 
becomes nonzero. For comparison we note that the average parton density in a baryon is about 
1.5 fm~ 3 while it is about 1.25 fm -3 in case of a meson. Gluons in PHSD decay to color octet 
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qq pairs below p c ; their reformation rate is strongly suppressed at low parton densities due to 
the large gluon masses from the DQPM (cf. Fig. 25). 

On the hadronic side PHSD includes explicitly the baryon octet and decouplct, the 0~ 
and 1 _ meson nonets as well as selected higher resonances. Hadrons of higher masses (> 1.5 
GcV in case of baryons and > 1.3 GeV in case of mesons) are treated as 'strings' that reflect 
the continuum excitation spectrum of mesons or baryons and decay to the known (low mass) 
hadrons within HSD [79] according to the JETSET algorithm [185]. 



4.2.1 The dynamics of an expanding partonic fireball 

Due to the complexity of the PHSD approach it is illustrative to explore the parton and 
hadronization dynamics in a transparent model case [143] rather than directly step to rela- 
tivistic nucleus- nucleus collisions at SPS or RHIC energies. Accordingly the following study 
addresses the expansion dynamics of a partonic fireball at initial temperature T = 1.7 T c (T c = 
0.185 GeV) with quasiparticle properties and four- momentum distributions determined by the 
DQPM at temperature T/T c = 1.7 for quark chemical potential /x 9 =0. The initial distribution 
in coordinate space is taken as a Gaussian ellipsoid with a spatial eccentricity 



in order to allow for the build-up of elliptic flow (as in semi-central nucleus- nucleus collisions 
at relativistic energies). In order to match the initial off-equilibrium strange-quark content 
in relativistic pp collisions the number of s (and s quarks) is assumed to be suppressed by 
a factor of 3 relative to the abundance of u, d quarks and antiquarks. In this way we will 
be able to investigate additionally the question of strangeness equilibration. Recall that the 
dynamical evolution of the system now is entirely described by the transport dynamics in 
PHSD incorporating the off-shell propagation of the partonic quasiparticles according to Ref. 
[65] (cf. Section 2) as well as the transition to resonant hadronic states (or strings) in line 
with (161) which is solved by Monte Carlo sampling on a time-depending expanding grid with 
local cells of volume dV(t)= 0.25 (1+0.7 t) 3 fm 3 . This choice approximately corresponds to a 
comoving grid for the expanding system. 

In Fig. 50 (upper part) the energy balance for the expanding system at initial temperature 
T = 1.7T C , ji q = and eccentricity e =0 is shown. The total energy E to t (upper line) is conserved 
within 3% throughout the partonic expansion and hadronization phase such that for t > 10 fm/c 
it is given essentially by the energy contribution from mesons and baryons (+antibaryons). The 
initial energy splits into the partonic interaction energy V p and the energy of the time-like 
partons 



with fractions determined by the DQPM (Section 3). The hadronization mainly proceeds during 
the time interval 1 fm/c < t < 7 fm/c as can be extracted also from the lower part of Fig. 50 
where the time evolution of the quark (+antiquark), gluon, meson and baryon (+antibaryon) 
number is displayed. 

The hadronic decomposition is dominantly determined via parton fusion to mesonic or bary- 
onic resonances and strings and their hadronic decays. Since on average the number of hadrons 
from the resonance or string decays is larger than the initial number of fusing partons, the 
hadronization process leads to an increase of the total entropy and not to a decrease as in case 
of coalescence models [186,187]. This is a direct consequence of the finite quark and antiquark 
masses which - by energy conservation - lead to hadron masses above 0.8 GeV (1.3 GeV) for 
meson and baryonic states, respectively; this solves the entropy problem in hadronization in a 
natural way\ 



e =<2/ 2 -* 2 >/<2/ 2 +- 2 } 



(163) 




(164) 
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time [fm/c] 

Fig. 50. Upper part: Time evolution of the total energy E to t (upper line), the partonic contributions 
from the interaction energy V p and the energy of time-like partons T v in comparison to the energy 
contribution from formed mesons Em and baryons (+ antibaryons) E B+ g. Lower part: Time evolution 
in the parton, meson and baryon number for an exploding partonic fireball at initial temperature 
T = 1.7 T c with initial eccentricity e = 0. 

4.2.2 Hadronization 

It is, furthermore, interesting to have a look at the final particle ratios K + /tt + , p/ir + , A/K + 
etc. (after hadronic decays) which are shown in Table 1. The latter ratios may be compared to 
the grandcanonical statistical hadronization model (SM) at baryon chemical potential /is = 
[188-190]. In this particular case the particle ratios only depend on temperature T and one may 
fix a freeze-out temperature, e.g., by the proton to 7r + ratio. A respective comparison is given 
in Table 1 for T = 160 MeV and 170 MeV for the SM which demonstrates that the results from 
PHSD are close to those from the SM for T w 170 MeV. This also holds roughly for the A/K+ 
ratio. On the other hand the K + /n + ratio only smoothly depends on the temperature T and 
measures the amount of strangeness equilibration. Recall that the system has been initialized 
with a relative strangeness suppression factor of 1/3. The deviation from the SM ratio by about 
13% indicates that strangeness equilibration is not fully achieved in the calculations. This is 
expected since the partons in the surface of the fireball hadronize before chemical equilibration 
may occur. 

The agreement between the PHSD results for the baryon to meson ratio in the strangeness 
S=0 and S—l sector may be explained as follows: Since the final hadron formation dominantly 
proceeds via resonance and string formation and decay - which is approximately a micro- 
canonical statistical process [191] - the average over many hadronization events with different 





p/ir + 


A/K+ 


K+/TT+ 


PHSD 


0.086 


0.28 


0.157 


SM T = 160 MeV 


0.073 


0.22 


0.179 


SM T = 170 MeV 


0.086 


0.26 


0.180 



Table 1. Comparison of particle ratios from PHSD with the statistical model (SM) [190] for T= 160 
MeV and 170 MeV. 
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Fig. 51. Time evolution of the elliptic flow V2 for partons and hadrons for the initial spatial eccentricity 
e = 0.33 for an expanding partonic fireball at initial temperature T — 1.7 T c . 

energy /mass and particle number (in the initial and final state) leads to a grandcanonical en- 
semble. The latter (for [i^ — 0) is only characterized by the average energy or an associated 
Lagrange parameter (3 = 1/T as well known from quantum statistical mechanics. 

4.2.3 Elliptic flow 

Of additional interest are the collective properties of the partonic system during the early time 
evolution. In order to demonstrate the build-up of elliptic flow we show in Fig. 51 the time 
evolution of 

V2 = ((pI-pI)/( P 2 x +P 2 v )) ( 165 ) 

for partons (solid line), mesons (long dashed line) and baryons (dot-dashed line) for an initial 
eccentricity e = 0.33. As seen from Fig. 51 the partonic flow develops within 2 fm/c and the 
hadrons produced in time essentially pick up the collective flow from the accelerated partons. 
In total the hadron v-i is smaller than the maximal parton i>2 since by parton fusion the average 
V2 reduces and a fraction of hadrons is formed early at the surface of the fireball without a 
strong acceleration before hadronization. Is is worth to point out that in PHSD the elliptic flow 
of partons predominantly stems from the gradients of the repulsive parton mean-fields at high 
parton (energy) density. This is quite analoguous to A + A collisions at SIS energies where the 
collective flow is dominated by the nucleon mean-field potentials [192]. To demonstrate this 
statement we display in Fig. 51 the result for a simulation without elastic partonic rescattering 
processes by the short dashed line. 

Fig. 52 shows the final hadron vi versus the initial eccentricity e and indicates that the 
ratio vijt is practically constant (ss 0.2) as in ideal hydrodynamics (cf. Fig. 3 in Ref. [193]). 
Accordingly the parton dynamics in PHSD are close to ideal hydrodynamics. This result is 
expected since the ratio of the shear viscosity r) to the entropy density s in the DQPM is on 
the level of rj/s 0.2 [97] and thus rather close to the suggested lower bound oirj/s — l/(47r) 
[194]. 

In summarizing this section we like to point out that the expansion dynamics of an anisotropic 
partonic fireball within the PHSD approach - including dynamical local transition rates from 
partons to hadrons (and vice versa) - shows collective features as expected from ideal hydro- 
dynamics in case of strongly interacting systems. The hadronization process conserves four- 
momentum and all flavor currents and slightly increases the total entropy (by about 15% in 
the model case investigated here) since the 'fusion' of massive partons dominantly leads to the 
formation of color neutral strings or resonanaces that decay microcanonically to lower mass 
hadrons. This solves the entropy problem associated with the simple coalescence model in case 
of massless partons! 

Furthermore, the hadron abundancies and baryon to meson ratios are found to be compati- 
ble with those from the statistical hadronization model [188, 189] - which describes well particle 
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Fig. 52. The elliptic flow vi versus the initial spatial eccentricity e for an expanding partonic fireball 
at initial temperature T = 1.7 T c . 

ratios from AGS to RHIC energies - at a freeze-out temperature of about 170 MeV. Strangeness 
equilibration is approximately achieved in the dynamical expansion and driven by the processes 
qq «-> g ss, which is a resonant process in the DQPM. Nevertheless, some note of caution 
has to be added here: Although the final hadron ratios are compatible with a fixed freeze-out 
temperature (~ 170 MeV) we find that the actual hadronization occurs at very different energy 
densities (or temperatures) (cf. also Ref. [195]) such that the microscopic studies do not support 
the sudden freeze-out picture. 

In closing these lectures the author likes to point out that the field of quantum physics out- 
of-equilibrium and especially for strongly interacting systems is far from being closed. Though 
substantial progress has been achieved in the last decade many fundamental questions still lack 
a satisfactory answer. This subject is challenging but also exciting! 
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